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ABSTRACT 


•4 

The  pickup  and  delivery  problem  (PUDP)  represents  a  class  of 
sequencing  or  routing  problems  where  the  key  facet  of  the  routing  is 
that  a  pickup  must  precede  the  corresponding,  subsequent  delivery. 

Other  considerations  such  as  service  time  windows,  quality  of  service 
parameters  or  operational  constraints  on  either  the  driver  o.-  the 
vehicle  are  possible.  As  such,  the  PUliP  is  a  constrained  version  of 
the  ubiquitous  travelling  salesman  problem  (TSP),  which  seeks  a 
minimum  cost  route  that  from  an  initial  point  visits  each  city  or 
stop  once  and  only  once,  ending  at  the  initial  stop.  There  arc  also 
similarities  between  the  PUliP  and  the  much  studied  vehicle  routing 
problem  (VRP),  although  the  two  problems  arc  distinctly  different 
because  of  the  origin  preceding  destination  requirement . 

The  TSP  and  VRP  literature  is  extensive,  offering  both  theory  and 
algorithms  for  the  solution  of  these  problems.  Given  the  similarities 
of  these  problems  to  the  PUDP,  those  algorithms  that  performed  well  on 
TSP’s  and  VRP's  arc  discussed  in  detail  and  served  as  the  basis  for 
developing  both  exact  and  heuristic  algorithms  to  solve  the  PUDP. 

■Assuming  that  all  problem  constraints  arc  cxprossablc  in  terms  of 
stop  numbers  along  a  vehicle's  route,  dynamic  programing  can  be  used 
to  optimally  solve  the  problem.  The  algorithm  developed  is  signifi¬ 
cantly  more  powerful  on  heavily  constrained  problem  instances  than  any 
other  known  technique.  Solutions  to  problems  in  excess  of  45  customers 
(equivalent  to  a  91  city  TSP)  arc  solved  on  an  IBM  3031  computer  in  a 


matter  of  seconds.  The  efficiency  is  achieved  by  only  generating 


feasible  state  spa qo,  vectors,  thus  greatly  reducing  the  storage  and 
execution  requirements.  The  same  dynamic  programming  algorithm  is  also 
used  to  solve  the  multiple  vehicle  PUDP  but  with  less  impressive 
results. ^  Other  exact  techniques  could  not  be  effectively  used  on  the 
PUDR-dorf^ primarily  to  the  precedence  requirement. 

-  Heuristic  algorithms  were  also  developed  and  tested.  Most  of  the 
algorithms  commonly  used  to  solve  tho  related  TSP  and  PUDP  perform 
poorly  on  the  PUDP,  often  producing  solutions  as  much  as  50»  above 
optimal.  v  An  interchange  (3-optimal)  heuristic  consistently  produced 
superior  results  for  the  single  vehicle  PUDP.  Solutions  averaging 
within  1*  of  optimal  were  obtained  for  heavily  constrained  problem 
instances. 

.  The  multiple  vehicle  problem  is  significantly  more  complex  than 
is  the  single  vehicle  problem.  Results  for  the  multiple  vehicle 
problem  were  acceptable  but  inconclusive.  Consequently,  the  multiple 
vehicle  area  is  judged  to  be  the  most  promising  area  for  future 


research. 
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ABSTRACT 


The  pickup  and  delivery  problem  (I'UDI'J  represents  a  class  of 
sequencing  or  routing  problems  where  the  key  facet  of  the  routing  is 
that  a  pickup  must  precede  the  corresponding,  subsequent  delivery. 

Other  considerations  such  as  service  time  windows,  quality  of  service 
parameters  or  operational  constraints  on  either  the  driver  or  the 
vehicle  are  possible.  As  such,  the  rillil’  is  a  constrained  version  of 
the  ubiquitous  travelling  salesman  problem  (ISP),  which  seeks  a 
minimum  cost  route  that  from  an  initial  |ioint  visits  each  city  or 
stop  once  and  only  once,  ending  at  the  initial  stop,  (here  are  also 
similarities  between  the  1*111*1*  and  the  much  studied  vehicle  routing 
problem  (VKP),  althou.-.h  the  two  problems  are  distmctl>  different 
because  of  the  origin  preceding  destination  requirement. 

The  TSP  and  VHP  literature  is  extensive,  ottering  both  theory  and 
algorithms  for  the  solution  of  these  problems.  Given  the  similarities 
of  these  problems  to  the  I1IUP,  those  algoiithms  that  performed  well  oil 
ISP's  and  VHP's  are  discussed  m  detail  and  served  as  the  basis  for 
developing  both  exact  and  heuristic  algorithms  to  solve  the  PIIPP. 

Assuming  that  all  problem  constraints  are  cxprcssahlc  in  terms  of 
stop  numbers  along  a  vehicle's  route,  dynamic  programming  can  be  used 
to  optimally  solve  the  problem.  Jhc  algorithm  developed  is  signifi¬ 
cantly  more  powerful  on  heavily  constrained  problem  instances  titan  any 
other  known  technique.  Solutions  to  problems  in  excess  of  4S  customers 
(equivalent  to  a  91  city  TSP)  are  solved  on  an  IBM  3051  computer  in  a 

i  i 


4 


i  i  i 

matter  of  seconds.  The  efficiency  is  achieved  by  only  generating 
feasible  state  space  vectors,  thus  greatly  reducing  the  storage  and 
execution  requirements.  The  same  dynamic  programming  algorithm  is  also 
used  to  solve  the  multiple  vehicle  I'UDI’  hut  with  less  impressive 
results.  Other  exact  techniques  could  not  he  effectively  used  on  t lie 
PUDI*  die  primarily  ta  the  precedence  requirement. 

Ihuristic  algorithms  were  also  developed  and  tested.  Most  of  the 
algornhms  commonly  used  to  solve  the  related  t SI*  and  perform 

poorly  «n  the  I'UDI’,  often  producing  solutions  as  much  as  .%(>'.  above 
optimal  An  interchange  (.t-optim.il)  heuristic  consistently  produced 
super i>  r  results  for  the  single  vehicle  HUH'.  Solutions  averaging 
within  i'  of  optimal  were  obtained  for  heavily  constrained  problem 
instances. 

Tlie  multiple  vehicle  problem  is  significantly  more  complex  than 
is  the  single  vehicle  problem.  Results  for  the  multiple  vehicle 
problem  sere  acceptable  but  inconclusive.  Consequently,  the  multiple- 
vehicle  area  is  judged  to  be  the  most  promising  area  for  future 


research . 
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OVliRVlI.K  OR  HU:  PiCkUP  AND  DCLIVhRY  PROBLEM 


1  .  PROBLEM  CUSS i  1- 1 CAT  1  ON 

The  pickup  and  delivery  problem  (PUOP)  represents  a  class  of 
sequencing  or  routing  problems  where  the  key  facet  of  the  routing  is 
that  a  pickup  must  precede  the  corresponding,  subsequent  delivery. 

Many  other  constraints  arc  possible  based  on  the  particular  applica¬ 
tion.  As  such,  the  pickup  and  delivery  problem  is  a  constrained 
version  of  the  basic,  much  studied  travelling  salesman  problem  (TSP). 
The  travelling  salesman  problem  seeks  to  find  a  minimum  cost  path 
that  from  an  initial  point,  visits  each  city  or  stop  once  and  only 
once,  ending  at  the  initial  stop.  More  rigorous  definitions  arc 
provided  below  and  in  Chapter  11. 

The  vehicle  routing  problem  (VRP)  is  also  a  constrained  version 
of  the  TSP.  In  the  vehicle  routing  problem  the  key  consideration  is 
vehicle  capacity,  although,  as  with  the  PUOP,  other  constraints  may  be 
applicable.  It  is  not  correct  to  classify  the  PUOP  as  a  further 
constrained  vehicle  routing  problem.  Vehicle  capacity  need  not  be  a 
factor  in  the  PUOP.  Given  that  capacity  is  a  consideration,  the 
manner  in  which  it  affects  the  problem  is  signi ficant ly  different 
between  the  two  problems.  Consequently,  it  appears  the  proper  classi¬ 
fication  is  to  treat  both  problems  as  constrained  versions  of  the  TSP. 


The  PUDP  is  representative  of  several  practical  routing 
situations.  Examples  include  dial-a-ride  services  and  courier  services. 
Notwithstanding,  few  articles  devoted  to  the  PUDP  have  appeared  in  the 
published  literature.  This  is  in  sharp  contrast  to  the  literature 
devoted  to  the  TSP  and  the  VRP  which  is  voluminous  and  extends  over  the 
last  quarter  century.  All  evidence  suggests  that  the  PUDP  is  a 
relatively  unexplored  subject  area. 

11.  PROBLEM  DESCRIPTION 

The  dial-a-ride  service  (DARS)  offers  a  convenient  means  of 
conceptualizing  the  PUDP.  Suppose  that  an  organization  provides 
transportation  services  for  the  handicapped.  Vehicles  must  pick  up 
these  people  at  their  individual  origins  and  take  them  to  their  desti¬ 
nations.  Return  trips  sometime  later  in  the  day  are  also  possibilities. 
The  objective  is  to  satisfy  all  requests  for  service  as  economically 
as  possible.  In  addition  to  the  origin-destination  (0/D)  constraint, 
other  intuitively  appealing  constraints  may  include: 

1.  A  limitation  on  the  number  of  passengers  who  can  occupy  the 
vehicle  at  any  one  time  (capacity  constraint); 

2.  A  limitation  on  the  amount  of  time  that  any  one  passenger  must 
remain  in  the  vehicle  (quality  of  service  constraint); 

3.  A  range  of  times  in  which  pickup  and/or  delivery  must  he  made 
(time  window  constraint); 

4.  A  requirement  that  the  same  driver  provide  both  legs  of  a 
person's  round  trip  (stop-back  constraint); 


5.  Limits  on  the  total  distance  that  a  vehicle  may  he  driven 
(operational  constraint). 

III.  PROBI.LM  S ICN 1 F I CANCL 

As  mentioned  above,  the  PUMP  is  representative  of  real  world 
routing  situations.  Fisher  and  Jaikumar  (10)  estimated  that  urban 
delivery  vehicles  in  1975  travelled  approximately  70  billion  miles  at 
a  fuel  cost  in  excess  of  $5.5  billion.  Fuel  costs  have  more  than 
doubled  in  the  last  five  years.  Consequently,  annual  fuel  costs  for 
urban  delivery  vehicles  in  excess  of  $10  billion  is  clearly  probable. 
The  percentage  of  urban  vehicles  engaged  in  services  that  could  be 
classified  as  fitting  the  pickup  and  delivery  model  is  not  known. 
However,  were  the  figure  as  low  as  5°,  a  figure  in  excess  of  $500 
million  is  obtained. 

Dial-a-ride  services  are  playing  an  increasingly  important  role 
in  urban  public  transportation.  Such  services  for  the  handicapped, 
or  transportationally  disadvantaged,  are  available  in  nearly  every 
American  city,  cither  provided  by  the  public,  by  charitable  organi¬ 
zations,  or  by  volunteers.  Unquestionably,  a  5.  improvement  in  route 
efficiency  would  produce  a  tremendous  savings  in  both  dollars  and 
barrels  of  oil . 

experience  with  the  vehicle  routing  problem  suggests  savings  of 
from  5"o  to  10°.  are  possible  bv  applying  fairly  simple  computer  assisted 
algorithms  to  route  the  vehicles.  Because  the  PUMP  is  inherently 
more  complex,  it  appears  less  likely  that  a  dispatcher,  acting  without 


benefit  of  some  algorithm,  could  produce  a  good  route.  Consequently, 
even  a  greater  savings  percentage  appears  possible  for  the  PLJUP. 

IV.  RELATIONSHIP  TO  THE  TSP  AND  THE  VRP 

There  are  similarities  as  well  as  differences  among  the  PUDP,  the 
TSP,  and  the  VRP.  Since  the  TSP  and  VRP  have  been  extensively  investi¬ 
gated,  exploitation  of  the  similarities  was  a  logical  course  of  action. 

Let  G  =  {N,A,C}  be  a  complete  network  with  N  representing  the  set 
of  nodes,  A  the  set  of  arcs,  and  C  =  |c..|  a  matrix  of  costs  represent¬ 
ing  the  cost  of  going  from  node  i  to  node  j.  A  Hamiltonian  cycle  is  a 
cycle  that  passes  through  each  node  it-N  exactly  once.  The  TSP  is  the 
problem  of  finding  a  least  cost  Hamiltonian  cycle  on  G.  The  multiple 
travelling  salesman  problem  (MTSP)  requires  that  for  m  salesmen  one 
find  m  cycles  on  G  such  that  every  itN  is  visited  exactly  once  and 
the  total  cost  of  the  m  cycles  is  minimal. 

If  we  further  constrain  the  MTSP  by  requiring  that  for  any  cycle 
or  tour  the  capacity  (in  weight,  volume...)  of  the  corresponding 
vehicle  cannot  be  exceeded,  the  resulting  problem  is  the  VRP.  Although 
other  constraints,  such  as  time-windows  or  stop-backs,  are  possible, 
very  few  articles  address  them  even  in  passing.  One  notable  exception 
is  the  paper  by  Fisher  and  Jaikumar  (10),  which  explicitly  considers 
the  time-window  constraints. 

In  none  of  these  problem  definitions  is  there  any  mention  of  a 
requirement  that  one  stop  be  visited  before  another.  Incorporating 
this  procedence  relationship,  required  to  define  the  PUDP,  complicates 
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the  TSP  or  MTSP  much  more  than  docs  capacity.  The  PUDP  is  considered 
more  complicated  and  at  least  as  difficult  as  any  of  these  related 
problems . 


V.  PR0BL1M  D1ITIOJLTY 

The  TSP  has  been  shown  to  be  NP-complcte  (12).  One  consequence 
of  this  classification  is  that  there  is  no  known  polynomial  time 
algorithm  that  optimally  solves  it,  despite  hundreds  of  man-years 
devoted  to  finding  one.  If  one  could  find  a  polynomial  time  algorithm 
for  any  one  of  the  more  than  300  NP-complcte  problems,  one  could  solve 
all  problems  in  NP  in  polynomial  time.  As  the  following  lemma  shows, 
the  PUDP  is  at  least  as  hard  as  the  TSP.  If  a  polynomial  time 
algorithm  exists  that  solves  the  PUDP,  it  could  be  used  to  solve  the 
TSP  in  polynomial  time.  This  in  turn  implies  that  it  could  be  used  to 
solve  all  of  the  other  NP-completc  problems.  Consequently,  the  likeli¬ 
hood  of  anyone  finding  a  polynomial  time  algorithm  for  the  PUDP  is  not 
considered  high. 

Lemma  1 

Unless  P  =  NP,  there  docs  not  exist  a  polynomial  time  algorithm 
that  optimally  solves  the  general  PUDP. 

Proof 

It  will  be  shown  that  one  instance  of  the  PUDP  can  be  reduced  to 
two  travelling  salesmen  problems.  Since  the  TSP  is  NP-complcte  and 
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cannot  be  solved  in  polynomial  time  unless  1’  =  NP,  this  implies  that  the 
PUDP  can  not  be  solved  in  polynomial  time  unless  P  =  NP. 

Consider  the  instance  where  all  origins  are  to  be  visited  before 
lunch,  all  destinations  after  lunch,  and  the  driver  must  return  to  the 
depot  for  lunch.  Clearly  the  optimal  solution  for  the  PUDP  is  the 
optimal  sequencing  over  the  set  of  origins  coupled  to  the  optimal 
sequencing  over  the  set  of  destinations.  But  these  optimal  solutions 
are  the  TSP  solutions  taken  over  their  respective  sets.  Hence,  unless 
P  =  NP,  a  polynomial  time  algorithm  for  the  general  PUDP  does  not 
exist . 


VI.  Rl.ShArCi:  JULSTIONS 

The  TSP  and  VRP  literature  offers  an  abundance  of  theory  and 
algorithms  for  the  solution  of  *  'O'r  respective  problems.  Given 
the  similarities  between  tl.  1  UUP  and  the  VRP  and  TSP,  one  would  suspect 
that  many  of  the  same  algorit!  <s  could  be  used  for  the  PUDP.  To 
what  extent  and  how  well  remained  to  be  determined.  Consequently,  this 
research  was  guided  by  the  following  three  general  questions: 

1.  Using  what  algorithms,  and  under  what  conditions,  can  the 
PUDP  be  efficiently  solved  optimally? 

2.  How  well  do  heuristics  commonly  applied  to  related  problems 
perform? 

3.  Which  heuristics  (s)  provide  the  best  PlIDP  solutions? 


i 


VI 1.  ORGAN  HAT  I  ON 


Chapter  11  provides  a  review  of  the  literature  pertinent  to  the 
pickup  and  delivery,  travelling  salesman,  and  vehicle  routing  problems. 
Those  algorithms,  both  exact  and  heuristic  which  have  been  successfully 
applied  to  the  TSP  and  the  VRP  are  developed  in  detail. 

In  Chapter  111,  a  detailed  mathematical  formulation  for  t lie  PUDP 
is  presented.  The  formulation  provides  the  insight  necessary  to 
explain  why  some  of  the  algorithms  developed  in  Chapter  II  do  not  work 
efficiently  when  applied  to  the  PUDP,  while  others  do. 

Chapter  IV  deals  with  the  optimal  solution  to  the  single  vehicle 
problem.  A  dynamic  programming  solution  is  developed  which  is  cxtrcmcl 
powerful  when  applied  to  highly  constrained  PUDP  instances.  The  same 
dynamic  programming  model  is  also  used  to  solve  the  multiple  vehicle 
PUDP  discussed  in  Chapter  VI.  Lxact  solutions  other  than  by  the 
dynamic  programming  approach  arc  shown  to  be  much  more  difficult  to 
obtain  for  the  PUDP  than  for  the  TSP. 

Single  vehicle  heuristics  are  discussed  in  Chapter  V.  Many  of  the 
heuristics  that  are  widely  applied  to  the  VRP  are  shown  to  perform 
poorly  on  the  PUDP,  especially  when  side  constraints  become  more 
binding.  Special  attention  is  given  to  those  instances  of  the  PUDP 
for  which  the  dynamic  programming  technique  provides  an  optimal 
solution.  This  allows  for  a  precise  evaluation  of  how  well  a  given 
heuristic  performs.  Only  relative  performance  has  previously  been 
possible  for  all  but  extremely  small  vehicle  routing  problems. 


The  multiple  vehicle  PUDP,  discussed  in  Chapter  VI,  is  seen  to  be 
a  most  difficult  problem.  The  precedence  requirements  render  ineffi¬ 
cient  the  heuristic  determined  to  be  the  most  powerful  for  the  single 
vehicle  problem.  Further,  the  additional  alternatives  available  due  to 
more  than  one  vehicle  being  available  for  customer  assignment  destroy 
the  efficiency  of  the  dynamic  programming  technique.  Heuristics  that 
proved  somewhat  successful  as  well  as  those  that  failed  are  detailed. 
When  failure  is  encountered,  an  explanation  is  offered. 

During  the  period  of  the  research,  many  areas  ripe  for  research 
were  encountered.  Practical  limitations  on  available  time  precluded 
investigation  of  these  opportunities  as  a  part  of  this  effort. 

Chapter  VII,  therefore,  details  several  of  the  more  interesting  areas 
remaining  to  be  explored. 
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L1TLRATURL  RIA1LK 

The  pickup  and  delivery  problem  (I’UUP)  has  been  classified  as  a 
constrained  version  of  the  travelling  salesman  problem  (ISP).  The 
vehicle  routing  problem  (VRP)  is  also  a  constrained  TSP.  Any  attempt 
to  solve  the  I’UDP  must,  therefore,  logically  begin  with  an  examination 
of  those  algorithms  successfully  used  in  solving  these  two  related 
problems.  This  chapter  summarises  the  results  of  such  an  examination. 
The  few  articles  directly  relating  to  the  I’UDP  are  also  discussed. 

The  travelling  salesman  problem  and  the  related  vehicle  routing 
problem  arc  two  of  the  most  studied  problems  in  operations  research. 
Literally  hundreds  of  algorithms,  many  representing  minor  modifications 
of  others,  have  been  proposed  for  their  solution.  It  would  neither  be 
practical  nor  useful  to  attempt  to  address  all  of  these  algorithms. 
Rather,  a  more  useful  approach  suggests  discussing  the  basic  approaches 
and  underlying  concepts  of  those  techniques  that  have  shown  the 
greatest  success  in  solving  the  TSP. 

1.  MATIILMATICAL  FORMULATION 

Travelling  Salesman  Problem 

A  word  description  of  the  ISP  was  presented  in  Chapter  I. 
Specifically,  the  TSP  seeks  to  minimize 
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where, 


c-j  =  cost  of  Roinj;  directly  from  city  i  to  i,  and 


xu  =  V 


1,  if  the  salesman  ^oes  directly  from  city  i  to  city  j 
jO,  otherwise. 

Expressions  (2.2)  and  |2.o)  insure  that  the  '-alesman  visits  and 
departs  from  each  city  exactly  once.  Express  ions  (2.1)  and  (2.o)  are 
the  subtour-el iminat ion  conditions  derived  originally  by  Miller, 
Tucker  and  Zcmlin  (2b)  and  often  quoted  by  others. 

Another  common  expression  for  el iminat inj;  subtours  is 
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for  every  nonempty  proper  subset  (j  of  N,  where  is  the  complement  of 
ij.  Bellmore  and  Malone  (2j  showed  that  this  formulation,  which  is  of 
order  2n  leads  to  an  effective  solution  algorithm  by  branch  and  bound. 
The  branch  and  bound  algorithm,  as  well  as  other  solution  approaches, 
is  discussed  later  in  this  chapter. 

Although  the  former  formulation  is  more  compact  than  the  latter, 
both  express  the  problem  precisely  in  mathematical  notation.  The 
compact  formulation  is  perhaps  more  elegant  and  may  be  of  value  compu¬ 
tationally  if  the  constraints  must  be  explicitly  considered.  The  less 
compact  formulation  is  sometimes  more  advantageous  if  the  constraints 
are  implicitly  handled,  lor  the  purpose  of  defining  the  1ST ,  either 
is  satisfactory,  lor  the  research  discussed  later,  there  is  no  compu¬ 
tational  preference  for  one  notation  over  the  other.  However,  the 
more  compact  formulation  is  used  throughout  simply  because  it  is 
not  at  iorial  ly  more  elegant. 

Multiple  Ira veil  mg  Salesman  Problem 

the  multiple  travelling  salesman  problem  (Ml  SI' )  can  be  easily 
transformed  into  a  standard  iSP.  Suppose  there  are  M  salesmen. 

M  copies  of  the  origin  are  included  in  the  cost  matrix  C.  =  {e  .  .  j ,  with 
each  of  the  copies  representing  a  unique  stop  but  with  the  same  costs 
relative  to  the  other  nodes  as  the  origin  lias.  Ilio  M  copies  are 
connected  with  arcs  of  infinite  lor  extremely  large  I  cost  so  that  it  is 
never  profitable  to  include  one  of  these  arcs  in  i  solution.  flic 


resulting  solution  to  the  ISP  taken  over  this  expanded  network  provides 


the  solution  to  the  MTSP.  Consequently,  subsequent  discussion  on  the 
TSP  is  equally  applicable  to  the  more  general  SUSP. 
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Vehicle  Routing  Problem 

The  VKP  formulation  is  similar  to  the  TSP  formulation  except  that 
it  is  necessary  to  include  the  number  of  vehicles,  k,  into  the  formula¬ 
tion  primarily  to  define  the  capacity  constraints,  <j^,  for  each  vehicle. 
Each  customer  has  a  requirement,  r^  An  integer  programming  formulation 
for  the  VHP  is: 
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)  ^  arbitrary  real  number,  i  =  l,2,...,ii  (2,15) 

where 

I 

,  1,  1 1‘  vehicle  k  visits  eusto::.er  ,i  immediately  after 

* j j  =  '  visiting  customer  i 
0,  otherwise. 

Many  of  these  expressions  are  logical  extensions  of  those  for  the 
TSI’.  i.xpression  (2.'.*)  requires  that  if  a  vehicle  visits  a  customer,  it 
must  depart  from  this  same  location.  Expressions  (’.10)  and  (2.11) 
represent  the  capacity  and  operational  limitations  of  the  kth  vehicle, 
while  (2.12)  insures  that  a  vehicle  must  he  used  once  and  only  once. 

Ihc  formulation  given  above  is  refered  to  hereafter  as  the  standard  VKI'. 

Other  constraints,  such  as  delivery  windows,  are  seldom  addressed 
at  all.  One  notable  exception  is  the  paper  by  1 i slier  and  Jaikumnr 
(10).  which  explicitly  considers  the  delivery  time  window  constraints. 
Their  formulation  is  much  more  complex  and  is,  therefore,  not  included 
here . 

II.  EXACT  AEliORI  llt'ES 

finding  exact  solutions  to  ail  but  relatively  small  (less  than  50 
cities  l  travelling  salesman  problems  has  proven  to  be  a  difficult  task. 
Exact  solutions  to  the  vehicle  routing  problem  are  much  more  difficult 
to  come  by.  Cbristof ides  was  credited  with  claiming  that  the  largest 
vehicle  routing  problem  of  any  complexity  that  had  been  solved  exactly 
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involved  only  23  customers  (18).  Christofides  subsequently  reports 
optimal  solutions  of  VRP's  of  up  to  about  30  customers  (4).  The  reason 
for  these  difficulties  is  the  exponential  growth  in  computations 
required  to  guarantee  optimality. 

lour  approaches  for  finding  exact  TSP  solutions  will  be  discussed. 
These  techniques  are  integer  programming  by  means  of  cutting  planes, 
branch  and  bound  based  on  subtour  elimination,  dynamic  programming, 
and  Lagrangcan  relaxation  using  minimum  1-trccs.  None  of  these 
approaches  arc  new,  the  most  recent  introduced  to  the  literature  in 
1970  (22).  In  addition,  a  Benders  decomposition  approach  for  solving 
the  vehicle  routing  problem  will  be  outlined. 

Integer  Programming 

Other  than  the  integrality  conditions  on  x..  in  expressions  (2.5) 
and  (2.14)  the  above  formulations  allow  for  solution  by  ordinary  linear 
programming  (LP) .  Suppose  these  conditions  are  relaxed  to 

0  x.j  <  1  for  every  i,j  .  (2.16) 

The  I.P  solution  will  not  generally  be  integral.  However,  it  is 
possible  to  add  additional  constraints  to  the  final  LP  tableau  to 
eventually  obtain  integrality.  These  additional  constraints  are  called 
cutting  planes,  Garfinkel  and  Nemhauser  (13)  provide  a  treatment  of  the 
theory  of  cutting  planes.  The  original  theory  of  cutting  planes  is  due 
ma inly  to  Gomo r y  (19). 


The  basic  idea  of  a  cutting  plane  algorithm  is  to  "cut  away," 
using  hyperplanes,  the  nonintegcr  portions  of  the  feasible  convex  hull 
of  the  relaxed  linear  program.  These  hyperplanes  are  constraints  that 
can  be  generated  at  each  step  from  the  current  LI'1  tableau,  and  taken  in 
such  a  manner  that  no  feasible,  integer  solutions  are  ignored. 

finite  cutting  plane  algorithms  have  been  proposed  and  used  to 
solve  the  ISP,  as  well  as  other  integer  programming  problems.  For  the 
most  part  they  have  not  performed  well,  the  exception  being  the  recent 
work  of  Miliotis  (25).  Fisher  and  Jaikomar  (10)  in  their  VRP  algorithm 
use  the  Miliotis  cutting  plane  algorithm  to  solve  TSP  subproblems  with 
impressive  results. 

Branch  and  Bound 

Branch  and  bound  techniques  are  by  far  the  most  common  type  applied 
to  the  TSP,  especially  those  employing  a  subtour  elimination  scheme. 
Expressions  (2.1),  (2.2),  (2.3)  and  (2.5)  taken  separately  define  the 
assignment  problem  (AP)  -  the  objective  being  the  most  efficient  assign¬ 
ment  of  n  men  to  n  distinct  jobs.  The  assignment  problem,  a  relaxation 
of  the  TSP,  is  easily  solved  in  polynomial  time.  Since  the  AP  is  a 
relaxation  of  the  TSP,  the  optimal  solution  to  the  AP,  which  is  most 
generally  not  feasible  to  the  TSP,  provides  a  lower  bound  on  the 
optimal  value  of  the  TSP  solution.  Any  feasible  TSP  solution  provides 
an  upper  bound.  If  the  AP  solution  is  not  feasible  to  the  TSP,  because 
of  subtours,  one  branches  into  k  subproblems,  where  k  is  the  number  of 
arcs  in  one  of  the  subtours.  In  each  subproblem,  one  of  the  k  original 
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subtour  arcs  is  assigned  an  infinite  cost,  this  breaking  or  eliminating 
the  subtour.  The  k  new  AP's  are  then  solved,  bounds  computed,  and  the 
process  repeated  if  subtours  are  present  and  the  lower  bound  is  less 
than  the  best  feasible  tour  as  yet  found.  Lventually,  one  is  assured 
of  finding  the  optimal  solution,  although  the  size  of  the  branching 
tree  may  become  enormous  for  large  problems. 

The  approach  outlined  above  represents  one  method  of  using  branch 
and  bound  techniques  to  solve  the  TSP.  The  key  to  the  success  of  a 
branch  and  bound  algorithm  rests  with  obtaining  very  tight  bounds,  thus 
greatly  reducing  the  size  of  the  branching  tree,  and  with  branching 
rules,  which  also  minimize  the  resulting  tree  size.  Branch  and  bound 
algorithms  often  find  optimal  or  near  optimal  solutions  toward  the 
beginning  of  the  enumeration.  Thus,  much  of  the  time  required  by  the 
algorithm  is  spent  in  verifying  the  optimality  of  a  tentative  solution. 
Consequently,  branch  and  bound  techniques  can  often  be  terminated  early, 
producing  a  very  good  solution.  In  this  sense,  the  technique  is  used 
as  a  heuristic.  Most  cutting  plane  algorithms  do  not  have  this 
characteristic . 

Dynamic  Programming 

Dynamic  programming  has  been  used  less  on  the  TSP  than  some  other 
techniques.  It  does  not  appear  that  dynamic  programming  has  been  used 
with  any  success  on  the  VHP.  The  computer  storage  requirements  arc  the 
primary  limiting  factor,  liven  though  the  dynamic  programming  recursions 
allow  for  treating  combinations  rather  than  permutations  of  the 
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operations  required,  the  combinations  increase  exponentially  with  the 
size  of  the  problem.  Consequently,  core  storage  requirements  to  solve 
a  20  city  TSP  exceed  900,000  words. 

The  basic  dynamic  programming  recursion  to  determine  a  shortest 
partial  TSt’  tour  from  the  origin  (node  1)  to  node  j  that  passes 
through  i  x ,  i2 .  ifc  ,  is 

fkO  1 1 1 .  i  . .  ik.,>  -»in  [<wyi,.i2 . i..!,!.., . ik. 

m=l, . . . ,k-l 

jl  (2.17) 

nr 


By  letting 


sk  =  •••*lk-i 


(2.18) 


we  can  simplify  (2.17)  to  obtain 


fk(j|Sk)  =  minltk_1(ini|Sk  . 
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(2.19) 


The  initial  recursive  equation  thus  becomes 


f,  (j  I  i  l ) 


ci.ii  +  c.  .  for  all  i i  ,  j  f  1  and  i  /  j,  (2.20) 

A  1  »  J  1 


while  the  final,  which  terminates  with  an  optimal  tour  value,  is 
obtained  by  solving 


f  (1  |S  j 
n  1  n  ’ 
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The  number  of  f^  values  is  given  by 

gln.k)  =  (n-lj!  /  (k-1 j ! (n-k-1 ) !  (2.22) 

which  reaches  a  maximum  halfway  through  the  computations. 

The  procedure  for  determining  the  optimal  tour  is  a  two  phase  one. 

First,  fj,.  k  =  2,3 . n  are  computed  recursiv  ly  by  (2.10).  Then 

the  optimum  ordering  (i i , i2 , . . . , i^)  is  obtained  by  picking  the  i  such 
that  (2.10)  holds  in  decreasing  order  of  k,  k  =  n,  n-1,  n-2,...,2. 

Lagrangean  Technique 

One  of  the  most  powerful  techniques  for  solving  the  symmetric  TSP 
is  the  1-tree  approach  of  Held  and  Karp  (22,23).  A  TSP  is  symmetric  if 

c.j  =  c..  for  every  i , j  .  (2.23) 

A  1-tree  is  a  tree  taken  over  vertices  2,3 . n,  connected  to  vertex 

1  with  two  edges.  A  TSP  tour  is  a  1-tree  for  which  each  vertex  lias 
degree  2.  Figure  1  depicts  a  1-tree  over  8  vertices.  A  minimum  weight 
1-tree  can  be  found  by  first  finding  a  minimum  spanning  tree  over  tire 
vertex  set  {2,3,...,n},  and  then  adding  the  two  least  cost  edges  at 
vertex  1.  The  minimum  spanning  tree  is  easily  solved  in  order  n  time. 
Changing  the  cost  matrix  for  intercity  distances  by 

c,.=c..+n.+tl.  (l 

ij  ij  i  j 

does  not  alter  the  optimal  solution  to  the  TSP  but  does  alter  the 
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solution  of  the  minimum  1-trce.  An  iterative  method  for  approaching 
the  optimal  solution  from  below  is  based  on  altering  the  and  tt^'s  so 
that  each  vertex  is  forced  toward  degree  2.  Normally,  the  algorithm 
must  employ  branch  and  bound  to  obtain  the  final  TS11  tour.  However, 
as  Held  and  Karp  note  the  bounds  computed  by  the  final  minimum  1-trees 
.  .are  so  sharp  that  the  search  trees  are  miniscule  compared  to  those 
normally  encountered.  .  .  ."  Consequently,  the  1-tree  approach  could  be 
considered  a  branch  and  bound  algorithm  that  uses  a  Lagrangean 
relaxation  to  compute  tight  lower  bounds.  It  is  important  to  note 
that  the  technique  is  only  valid  for  the  symmetric  travelling  salesman 
problem. 

Christofides ,  Mingozzi,  and  Toth  have  extended  Held  and  Karp's 
Lagrangean  relaxation  concept  to  the  VRP  (4).  It  is  this  work  for  which 
solutions  up  to  30  customers  are  reported.  The  key  factor  allowing  for 
success  is  again  tight  bounds  computed  by  similar  Lagrangean  penalty 
or  relaxation  procedures. 

Benders  Decomposition 

One  solution  technique  (10)  for  the  VRP,  which  specifically 
addresses  delivery  windows,  employs  Benders  decomposition.  The 
algorithm  iterates  between  a  generalized  assignment  problem  (GAP)  to 
assign  customers  to  vehicles  such  that  vehicle  capacity  is  not  exceeded, 
and  a  TSP  subproblcm  for  sequencing  the  stops.  As  mentioned  above,  a 
cutting  plane  algorithm  is  used  to  solve  the  TSP.  Infeasible  sub- 
problems  produce  the  Benders  cuts  needed  to  obtain  delivery  window 
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feasibility  on  subsequent  iterations.  Computational  results  remain 
incomplete,  but  appear  promising. 

fisher  and  Jaikumar  comment  that  although  their  algorithm  will 
produce  an  optimal  solution  in  a  finite  number  of  iterations,  they 
expect  that  for  practical  sice  problems  tiie  algorithm  will  be  used  as 
a  heuristic.  This  is  accomplished  by  terminating  the  iterations  prior 
to  achieving  optimality  and  using  the  best  feasible  solution  found  thus 
far . 

The  upper  limit  on  the  sice  of  a  general  TSi’  that  can  be  optimally 
solved,  within  practical  limits  on  time  and  computer  storage,  by  am- 
known  technique  is  questionable,  but  certainly  is  not  greater  than  100 
cities  and  is  probably  less.  Tor  the  VRP,  this  limit  is  s igni f icant ly 
lower.  Tor  problem  instances  of  such  size  that  obtaining  an  optimal 
solution  is  impractical,  heuristic  solutions  provide  the  only  alterna¬ 
tive  . 


III.  HEURISTIC  ALGORITHMS 

A  multitude  of  heuristic  algorithms  have  been  proposed  for  both 
the  TSP  and  the  VRP.  Many  heuristics  for  the  VRP  are  just  slight  modi¬ 
fications  of  ones  used  for  the  TSP.  The  heuristics  to  be  discussed 
fall  into  one  of  three  broad  classes:  tour  construction  approaches, 
tour  improvement  approaches,  and  composite  approaches.  The  latter  is 
the  logical  combination  of  the  first  two.  Each  class,  and  the  common 
representatives  of  each,  will  be  examined. 
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Tour  Construction  Heuristics 

Tour  construction  algorithms  generate  a  tour,  one  stop  at  a  time, 
from  the  given  distance  or  cost  matrix.  When  constraints  are  present, 
as  in  the  VRP,  the  next  stop  is  added  to  t lie  partially  constructed  tour 
only  if  it  can  lead  to  a  feasible  final  solution.  With  the  VRP,  this 
would  normally  be  a  check  to  insure  that  vehicle  capacity  had  not  been 
exceeded,  a  relatively  simple  procedure.  Insuring  feasibility  with 
delivery  windows  is  a  much  more  complicated  procedure  and  may,  to 
some  extent,  explain  why  such  constraints  are  not  treated  in  the 
literature. 

By  far  the  most  common  of  the  tour  construction  procedures  is  the 
Clarke-tVright  savings  algorithm  which  dominates  VRP  solution  attempts. 
Others  include  the  nearest  neighbor  or  greedy  approach  and  various 
insertion  approaches.  Since  the  literature  for  the  VRP  indicate  the 
general  superiority  of  the  Clarkc-W'right  models  over  other  tour  con¬ 
struction  heuristics,  they  will  be  discussed  first. 

Clarke-Wr i glit  savings.  The  basic  concept  of  the  travel  time  saved 
heuristic  was  developed  by  Clarke  and  Wri,  t  (?),  who  credited  the 
earlier  work  of  Hanttig  and  Rainscr  (8).  In  both  cases,  the  procedure 
was  developed  to  solve  the  VRP.  Numerous  modifications  have  been 
suggested,  but  the  underlying  concept  remains  invariant  III,  28,  48, 

42,  4.4)  . 

The  procedure  begins  by  selecting  one  node  as  the  origin.  With 
the  VRP,  no  choice  is  required.  Initially  one  assumes  that  everv  stop 
is  visited  directly  from  the  origin.  linn  the  savings  that  can  be 


achieved  by  combining  two  subtours  into  one,  by  link’ng  stops  i  and  j, 
is  computed  by 
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ij 


li 


c  •  • 


ij 


(2.25) 


Beginning  with  the  largest  of  these  savings  values,  routes  are 
assembled  such  that  the  next  stop  added  has  the  largest  remaining 
savings  —  provided  that  a  constraint  is  not  violated.  Customers  that 
have  been  linked  are  treated  as  a  single  macro  customer,  figure  2 
demonstrates  the  algorithm  for  two  stops  i  and  j.  Once  a  pair  of  stops 
has  been  linked,  they  remain  linked.  The  algorithm  continues  until  all 
stops  have  been  assigned. 

Most  commercially  available  algorithms  for  routing  vehicles  arc 
based  on  the  Clarke-Wr ight  savings  concept.  Survey  results  indicate 
a  percentage  in  excess  of  80°u ,  including  the  IBM  code,  YSPX,  which  is 
the  most  used  package  in  America  (18).  empirical  evidence  suggests 
that,  on  average,  the  Clarke-Wr ight  savings  method  produces  vehicle 
routes  that  are  as  good  as,  and  often  better  than,  routes  produced  by 
other  heuristics  for  the  standard  VRP.  Christofides  and  Li  Ion  (0) 
obtained  an  average  deviation  of  3.20„  from  optimal  on  10  small  vehicle 
routing  problems. 

One  of  the  attempts  to  improve  the  Clarkc-Wright  model  is  by 
modification  of  the  savings  equation.  Yellow  (13)  suggested  a  route 
shape  parameter  0  such  that  the  equation  is 
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(a)  Before  linking 


( b J  After  linking 


Figure  2.  Concept  of  the  Clarke-Kr ight 
savings  heuristic. 


Yarwtig  the  parameter  •:  alters  the  emphasis  placed  on  the  cost  between 
stops  i  and  i  in  relation  to  their  costs  relative  to  the  central  depot 
whin  =1,  Yellow's  model  reduces  to  the  l.’l  arke-Kr  i  ght  model, 
lypically,  several  values  of’  ••  are  tried,  including  v  =  1.  1  liere lore , 

Yellow's  method  always  produces  at  least  as  good  result  as  the 
i  1  arke  -  l\r  i  ght  technique. 

1  l  1  lman  and  Cochran  I  A'.>  I  and  'li  liman  and  Cain  138)  surest  ed 
another  approach  based  on  extending  the  savings  calculations  to  more 
than  one  stop  in  the  future.  Although  their  combined  reported  works 
are  limited  to  two  choices  in  sequence,  t lie  concept  is  extendable  to 
three  or  more  in  sequence.  However,  the  possible  combinations  are  of 
order  n‘  .  where  is  the  number  of  positions  examined.  Therefore,  thi 
procedure  rapidly  becomes  expensive  in  tents  of  computational  effort . 

Nearest  nejjlhhor.  The  nearest  neighbor  procedure  is  the  simplest 
and  perhaps  most  naive  of  all  heuristics  suggested.  One  always 
selects  the  nearest  unvisited,  feasible  stop  until  the  tour  is  com¬ 
pleted.  Because  of  the  way  in  which  the  algorithm  works,  it  is  often 
referred  to  as  a  greedy  algorithm.  I  he  primary  appeal  of  the  greed)' 
approach  is  its  simplicity.  It  is  easily  understood  and  large  problem 
can  easily  be  solved  by  hand.  Unfortunately,  the  greedy  approach  does 
not  normally  produce  "good"  tours  for  either  the  1ST  or  the  VKP. 

Insertion  procedures.  An  insertion  procedure  uses  a  specific 


selection  rule  to  pick  a  stop  not  yet  in  the  partial  solution  and  then 
determines  where  to  insert  this  stop.  Iluce  common  types  of  selection 


criteria  tor  the  next  stop  to  be  included  are  nearest,  farthest, 
and  random,  with  reference  to  any  one  of  the  nodes  included  in  the 
partial  tour.  Insertion  is  accomplished  for  node  k  by  finding  arc 
(i,j)  in  the  subtour  which  minimizes 

c.,  ♦  c,  -  c. . 

1  k  kj  i  j 

subject  to  the  problem  constraints,  if  any.  Computational  experience 
indicates  results  within  x-5".  of  optimal  or  best  known  solution  are 
attainable  with  insertion  algorithms  (1~J. 

lour  Improvement  lleurist ics 

The  basic  idea  behind  tour  improvement  procedures  is  to  take  a 
given  feasible  tour  and  systematically  modify  the  tour  to  obtain  a 
better  one.  The  procedure  is  one  of  arc  interchange.  'I wo  related 
heuristics  have  worked  well  on  the  TSP  and  the  VHP:  the  r-optimal 
heuristic  and  the  '-optimal  heuristic  of  bin  and  Kcrnighan  iJl). 

r-opt imaj .  Ihe  concept  of  r-opt ima 1 i t v  is  an  outgrowth  of  research 
on  the  ISP.  Christofidcs  and  i  t  ion  (<’l  appear  to  have  been  the  first 
to  apply  the  concept  to  the  vehicle  routing  problem,  with  results  at 
least  .is  good  as  those  obtained  using  t  lie  Cl  arke-Kr  i  ght  approach. 

Ihe  term  r-opt inal  implies  that  no  improvement  in  a  given  feasible 
solution  is  possible  by  eliminating  any  r  links  and  replacing  them  by 
r  new  links.  In  other  words,  the  procedure  terminates  at  a  local 
optimum.  The  r-optimal  procedure  is  an  O(n')  algorithm,  Consci|uent ly. 


only  2-optimal  and  3-optimal  interchanges  have  been  used  for  the  TSI’ 
and  Vltl'.  Larger  values  of  r  could  be  used  and  would  result  in  at  least 
as  good  a  solution  since  r  optimality  implies  r-1  optimality.  However, 
the  increased  computational  cost  could  apparently  not  be  justified. 

Christofidcs  and  Lilon  suggest  that  a  3-optimal  algorithm  first 
find  a  2-optimal  tour  and  then  use  this  tour  as  the  input  to  the 
3-optimal  algorithm.  Figure  3  shows  the  only  legitimate  reconnection 
pattern  for  the  2-optimal  algorithm  and  the  five  reconnection  patterns 
for  the  3-optimal  algorithm  which  exclude  duplication  of  a  2-optimal 
pattern.  In  Figure  3,  the  letters  represent  the  arcs  that  are  removed 
and  the  numbers  represent  the  specific  stops  that  are  to  be  reconnected. 
Initially,  the  stops  arc  visited  in  numerical  sequence,  with  inter¬ 
vening  stops  unnumbered.  The  heavy  lines  represent  all  intervening 
Stops  between  each  of  the  two  endpoints;  the  dotted  lines  the  single 
arc  removed;  and  the  remaining  lines  the  new  arcs.  Ihe  patterns  do  not 
depend  on  whether  the  network  is  directed  or  undirected.  However,  in 
the  directed  case,  some  of  the  patterns  may  he  infeasible  due  to  the 
nonexistence  of  an  arc  in  the  reverse  direction. 

* -opt  iron  1 .  Closely  related  to  the  r-optimal  procedure  is  the 
k -optimal  technique  of  I, in  and  kernighan.  Ihe  .‘-optimal  heuristic 
was  developed  to  solve  t lie  ISP  and  does  so  with  remarkable  results, 
ihe  probability  of  optimally  solving  a  given  IP  city  ISP  is  reportedly 
close  to  one,  while  an  optimal  solution  to  ,i  loo  city  ISP  is 
reported  at  close  to  23"o. 
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The  X-optimal  approach  works  essentially  in  the  same  manner  as  does 
the  r-optimal  approach  except  that  X  is  not  fixed  at  any  step  of  the 
iteration  process.  Rather,  for  each  iteration,  X  starts  at  a  value  of 
2,  i.e.,  two  arcs  are  initially  removed  and  a  2-optimal  solution  found. 
Then  a  third  arc  is  removed  and  a  3-optimal  solution  found.  Additional 
arcs  are  incrementally  removed  until  no  further  feasible  improvement 
in  total  tour  cost  can  be  made.  Consequently,  four,  five  or  even  more 
arcs  may  be  removed  at  one  time.  With  five  arcs  removed,  the  algorithm 
is  searching  for  a  better  solution  than  the  current  4-optimal  one. 

The  algorithm  insures  that  if  an  arc  is  removed,  a  feasible  reconnection 
pattern  does  exist,  thus  precluding  investigation  of  profitable,  but 
infeasible  reconnection  patterns.  The  creation  of  separate  subtours 
is  an  example  of  an  infeasible  reconnection  pattern. 

Composite  Heuristics 

A  composite  procedure  constructs  an  initial  tour  using  one  of  the 
tour  construction  techniques  and  then  attempts  to  improve  on  this 
solution  using  a  tour  improvement  technique.  Since  the  composite  will 
always  do  at  least  as  well  as  the  tour  construction  algorithm  without 
improvement,  the  question  of  whether  to  use  a  composite  is  one  of 
accuracy  desired  and  resources  available.  Several  unique  composite 
approaches  have  been  suggested,  specifically  for  application  to  the 
VRP.  Two  of  these  arc  of  special  interest  to  the  material  developed 
in  later  chapters. 
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Sweep  algorithm.  The  sweep  algorithm,  developed  by  Gillett  and 
Miller  (15),  employs  the  Lin-Kernighan  algorithm  to  sequence  stops  on 
cluster  generated  routes.  Other  models  of  the  vehicle  routing  problem 
can  be  characterized  as  solving  one  big  problem.  Gillett  and  Miller 
use  a  technique  which  divides  the  big  problem  into  k  subproblems, 
where  k  is  the  number  of  vehicles. 

The  algorithm  orders  all  delivery  points  sequentially  by  their 
polar  coordinate  angle  and  then  selects  routes  by  sweeping  through  the 
angles.  Points  are  not  added  if  one  of  the  constraints  would  be 
violated.  An  exchange  routine  is  applied  to  the  resultant  route 
structure  to  obtain  any  additional  route  improvement.  The  sweep  and 
improvements  are  accomplished  in  both  a  forward  and  a  backward  direction 
with  the  best  result  taken  as  the  final  solution.  Gillett  and  Miller 
recommended  and  used  the  Lin-Kernighan  algorithm  for  the  improvement 
step.  They  claim  that  problems  well  in  excess  of  100  customers  are 
computationally  feasible  with  their  sweep  algorithm,  and  solutions  are 
very  competative  with  other  solution  approaches. 

MTOUR.  MTOUR  is  the  name  given  to  an  algorithm  developed  by 
Russell  (36)  to  solve  constrained  multiple  travelling  salesman  problems 
or  constrained  VRP's.  The  additional  constraint  is  a  delivery  window 
or  a  due  date  requirement.  Russell  solves  the  multiple  problem  as  a 
single  problem  on  an  extended  network,  as  discussed  above,  using  the 
Lin-Kernighan  approach. 
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As  with  all  of  the  tour  improvement  heuristics,  an  initial 
feasible  solution  is  required.  In  the  more  heavily  constrained 
environment,  such  a  starting  solution  may  not  be  obvious.  In  the  more 
constrained  problem,  the  feasibility  of  arc  interchanges  greatly 
complicates  the  A-optimality  procedures.  Notwithstanding,  Russell  is 
one  of  the  few  who  does  address  the  time  window  constraint  and  he 
reports  good  success  in  using  his  MTOUR  on  his  test  problems. 

The  heuristics  discussed  here  are  not  intended  to  be  exhausive. 
However,  they  do  represent  the  basic  concepts  and  approaches  that 
have  been  successfully  applied  to  the  TSP  and  to  the  VRP.  Art.cles 
devoted  specifically  to  the  class  of  problems  described  by  the 
pickup  and  delivery  problem  are  far  less  numerous. 

IV.  PUDP  Rli  LATHI)  LITLRATURli 

Two  articles  and  one  working  paper  discuss  subject  areas  that 
fall  within  the  framework  of  the  PUDP  model.  Only  the  article  by 
Psaraftis  (34)  is  considered  significant  to  the  research  effort.  A 
brief  summary  of  the  other  two  is  presented  before  Psaraftis's  article 
is  discussed  in  more  detail. 

Single  Vehicle  PUDP 

In  their  1979  unpublished  paper,  Driscoll  and  Lmmons  (9)  discuss 
routing  of  a  single  vehicle  to  meet  pickup  and  delivery  requirements. 


The  only  problem  constraint  is  that  a  pickup  must  precede  a  subsequent 
delivery.  The  vehicle  is  taken  to  have  unlimited  capacity.  Driscoll 


and  Emmons  propose  a  heuristic  that  employs  a  2-optimal  algorithm.  All 
origins  are  2-optimal  sequenced,  and  then  all  destinations  are 
2-optimal  sequenced.  Finally,  the  two  lists  are  joined,  and  a  final 
2-optimal  sequence  is  formed. 

Dial-a-Ride  Systems 

Stein  (37)  presents  an  investigation  of  the  quality  of  any 
proposed  heuristic  in  a  theoretical,  asymptotic,  probabilistic  sense. 

No  specific  solution  is  offered.  All  of  the  results  cited  implicitly 
assume  vehicles  of  unlimited  capacity. 

Dynamic  Programming  and  Dial-a-Ride 

The  May  1980  edition  of  Transportation  Science  includes  the  article 
by  Psaraftis  that  outlines  the  use  of  dynamic  programming  to  optimally 
solve  the  single  vehicle,  many-to-many,  immediate-request  dial-a-ride 
problem.  This  article  appeared  subsequent  to  the  research  reported  in 
latter  chapters,  and  in  no  way  influenced  that  work. 

The  term  many-to-many  is  used  to  imply  that  the  pickup  points  and 
the  delivery  points  of  the  individual  customers  are  all  distinct  points. 
Immediate-request  is  used  to  imply  that  every  customer  wishes  to  be 
serviced  as  soon  as  possible.  The  article  addresses  both  the  static 
and  the  dynamic  cases.  Static  implies  that  no  additional  requests  for 
service  can  be  satisfied  by  the  single  vehicle  once  it  begins  servicing 


customer  requests. 


Objective. 


Psaraftis  uses  an  objective  function  that  is  a  weighted 


combination  of  total  service  time  for  all  customers  and  of  the  total 
service  time  for  all  customers  and  of  the  total  dissatisfaction  experi¬ 
enced  by  customers  in  waiting  to  be  picked  up  and  awaiting  delivery  at 
their  respective  destinations.  A  more  conventional  objective  function 
could  be  substituted  without  altering  the  algorithm.  A  maximum 
position  shift  (MPS)  limits  the  absolute  difference  between  the  order 
in  which  a  customer  is  served  on  the  route  developed  and  the  position 
a  customer  held  on  a  first  come  first  serve  list.  This  "position” 
orientation  to  the  problem  is  well  suited  for  solution  by  dynamic 
programming.  The  precedence,  MPS  and  capacity  constraints  are  all 
taken  care  of  as  the  recursive  procedure  develops. 


State  vector.  To  define  the  state  space,  a  state  vector 


(L, k i ,k2 , . . . , k.,)  is  used  where  L  represents  the  point  or  stop  the  vehicle 
is  currently  visiting,  N  is  the  number  of  customers  and  kj  is  the  status 
of  the  j^1  customers  such  that 

customer  j  has  yet  to  be  picked  up, 
customer  j  is  in  the  vehicle,  and 
customer  j  has  been  delivered. 

A  series  of  screening  steps  are  used  to  determine  if  a  particular  state 
vector  is  feasible  or  not.  One  of  the  screening  procedures  requires, 
for  each  state,  examining  subsequent  or  next  states  reachable  from  the 
current  state  vector  in  order  to  determine  if  the  present  vector  can  be 
feasibly  extended.  1'he  current  state  vector,  tentatively  classified  as 
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J 


k 


rnmtm 


34 


feasible,  is  reclassified  as  infeasible  if  all  subsequent  states  are 
determined  to  be  infeasible.  Consequently,  much  of  the  computational 
requirement  involves  these  screening  procedures. 

Storage  requirements.  Psaraftis'  dynamic  programming  algorithm 

N- 1 

requires  a  minimum  of  2  N  3  storage  locations,  where  N  is  the  number 
of  customers.  For  example,  for  IS  customers  a  total  of  over  430  million 
storage  locations  would  be  required.  As  actually  implemented,  Psaraftis' 
algorithm  would  require  in  excess  of  1  billion  storage  locations  to 
solve  a  15  customer  problem. 

Computational  results.  The  algorithm  runs  as  an  exponential 
function  of  the  size  of  the  problem  (number  of  customers),  but  is  shown 
to  be  asymptotically  more  efficient  than  dynamic  programming  applied 
to  the  travelling  salesman  problem.  The  largest  problem  reported  on 
involved  nine  customers  or  18  specific  stops  and  required  nearly  600 
seconds  to  solve.  Slight  improvements  in  running  times  were  noted  as 
the  MPS  and  capacity  values  were  decreased,  thus  making  the  problem 
more  heavily  constrained.  The  amount  of  improvement  between  an  uncon¬ 
strained  problem  and  a  fully  constrained  problem  (pickup  followed  by 
immediate  delivery  of  each  customer  in  first  come  first  serve  order) 
was  less  than  75",,. 

The  algorithm  as  constructed  can  only  be  used  to  solve  a  single 
vehicle  dial -a- ride  problem.  It  appears  to  be  practically  limited  to 


problems  of  nine  customers  or  less  due  to  storage  and  execution  require¬ 
ments.  In  Chapter  IV,  a  different  dynamic  programming  algorithm  is 
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developed.  This  dynamic  programming  algorithm  is  much  more  powerful 
than  Psaraftis'  is  on  the  more  heavily  constrained  problem.  Solutions 
to  problems  of  well  over  50  customers  are  readily  attainable.  The 
algorithm  also  can  be  used  to  solve  the  multiple  vehicle  problem. 


CHAPTER  III 

MATHEMATICAL  FORMULATION 

1.  MEMORY  REQUIREMENT 

In  Chapter  I  a  word  picture  of  the  pickup  and  delivery  problem 
(PUDP)  was  presented.  The  key  problem  facet  was  seen  to  be  the  prece¬ 
dence  relationship,  requiring  an  origin  be  visited  before  its  corre¬ 
sponding  destination.  Chapter  II  includes  the  mathematical  formulations 
of  the  related  travelling  salesman  and  vehicle  routing  problems.  These 
formulations  are  not  directly  extendable  to  describe  the  PUDP.  The 
precedence  relationship  and  any  capacity  restrictions  preclude  such  an 
extension . 

In  order  to  define  the  precedence  relationship  as  well  as  many  of 
the  other  suggested,  possible  constraints  for  the  PIJUP,  it  is  necessary 
to  include  in  the  model  a  memory  to  keep  track  of  the  sequence  or  timing 
of  the  stops.  Identifying  the  sequence  of  stops  is  required  to  insure 
that  a  destination  is  not  sequenced  before  its  origin.  Guaranteeing 
that  time  windows,  quality  of  service  standards  and  especially  the 
vehicle  capacity  limitations  are  met,  are  examples  of  t he  other 
constraints  that  are  time  or  sequence  dependent.  The  memory  provision 
can  be  achieved  by  incorporating  integral  time  periods  into  the  model. 
The  time  period  approach  is  mandatory  if  an  integer  linear  programming 
formulation  to  include  t lie  capacity  constraint  is  desired. 
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1 1 .  NOTAT I  ON 


The  following  notation  will  be  employed  through  the  remainder  of 
this  chapter: 


L 


n 

n  =  1./2 
N  =  n+1 
0 
I) 


K 

Mk 

R  =  I  r] 


number  of  nodes.  L  equals  the  number  of  origins  plus 
number  of  destinations  plus  one  for  each  vehicle's 
starting  and  ending  point,  bach  item's  or  person's 
origin  and  destination  is  assigned  a  unique  node. 
Consequently,  if  one  specific  location  were  to  serve  as 
a  destination  and  a  multiple  origin  of  multiplicity  tn, 
with  each  origin  having  a  different  destination,  m+1 
distinct  nodes  would  be  required.  There  would  be  zero 
cost  between  each  of  these  m+1  nodes.  The  depot  from 
which  all  vehicles  arc  dispatched  is  assigned  stop 
number  1 . 

number  of  origin/destination  pairs. 

#  vehicles. 

the  total  number  of  distinct  points. 

set  of  origins. 

set  of  destinations. 

number  of  vehicles. 

capacity  of  the  kth  vehicle. 

a  vector  whose  ith  component  represents  the  amount  or 
quantity  to  be  moved  from  origin  i  to  destination  1+n. 


We  assume  r. 
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1  if.O  for  the  PAHS. 
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-  cost  of  direct  travel  from  stop  i  to  stop  j.  The  cost  can 
he  distance,  time  or  other  suitable  measure. 

-  travel  time  from  stop  i  to  stop  j.  Since  travel  time  may 
be  a  valid  measure  of  cost,  it  is  possible  that  c.  .  =  i  .  . 
for  all  i  and  j.  However,  it  is  not  necessary  that  any 
specific  relationship  exist  between  these  values. 

-  earliest  time  that  stop  i  (either  a  customer's  origin  or 
destination)  can  be  serviced.  The  value  can  either  be 
expressed  in  terms  of  clock  time  or  by  a  period  number. 

-  latest  time  that  stop  i  can  be  serviced,  exprcssable  in  same 
manner  as  e.. 

i 

-  maximum  allowable  time  between  pick  up  and  delivery  for  the 
jth  customer  requirement.  can  either  be  expressed  in 
clock  time  or  by  the  number  of  intervening  periods. 

-  upper  limit  on  the  distance  or  other  appropriate  measure 
of  operational  limits  for  t he  kth  vehicle. 

j  1,  if  vehicle  k  proceeds  directly  from  stop  i  to  stop  \ , 

j 

)  where  j  is  the  vehicle's  tth  stop. 

i  0,  otherwise. 

-  clock  time  at  which  customer  stop  i  is  visited.  It  is 
assumed  that  the  starting  time  is  assigned  a  zero  time. 

-  a  scalar  taken  to  be  larger  than  the  length  of  any  feasible 
tour . 


-  set  of  time  periods. 


III.  TOUR  COM'  1  Nil  I  IV  CONS  I  RA I  MS 


As  with  the  multiple  travelling  salesman  and  the  vehicle  routing 
problems,  the  total  of  all  individual  tours  must  incorporate  all  N 
stops.  Because  of  the  parameter  t  in  the  formulation,  these  constraints 
can  be  expressed  more  compactly  than  they  can  be  by  using  the  more 
conventional  constraints.  As  before,  compactness  is  preferred  only  for 
mathematical  elegance. 

Convei itional  Constraints 

expressions  similar  to  those  for  the  travelling  salesman  and  vehicle 
routing  problems  can  be  used  to  force  tour  continuity  for  the  HUM'.  I  he 
problem  is 

minimice  c .  .  x ^  (Mil 

ijtk 


subject  to 


y .  arbitrary  real  number. 


(3.7) 


Expression  (3.2)  requires  that  every  stop  be  visited  exactly  once, 
while  (3.3)  states  that  if  vehicle  k  visits  a  given  stop,  it  must  also 
depart  from  it.  Expression  (3.-1)  insures  that  no  vehicle  is  used  more 
than  once.  The  subtour  elimination  equations  (3.5)  are  an  extension 
of  those  proposed  by  Miller,  'tucker  and  5emlin  for  the  the  travelling 
salesman  problem  (2(>) .  The  number  of  decision  variables  is  of  order 
KN3.  It  should  be  noted  that  with  the  exception  of  the  addition  of  the 
period  subscript  on  the  decision  variable  \ ,  (3.1)  to  (3.7)  are  identi¬ 
cal  to  those  defining  the  vehicle  routing  problem.  In  fact,  these  tom 
continuity  expressions  could  be  expressed  without  including  periodicity 
Periodicity  is  required  in  defining  other  constraints,  however,  as  will 
be  seen  later. 

Oomptic  t  ( amt  i  nu  i  t  y  Con s  t  r a  i  n t  s 

A  paper  by  fox,  (lavish  and  Craves  (11)  provides  a  compact 
formulation  of  order  N  for  the  time  dependent  travelling  salesman 
problem  (TUTSI’).  The  T 1  <  I S I  *  i a  variation  of  the  ISP  in  which 
C  =  [c.  |.  The  cost  of  going  from  i  to  i  depends  on  the  time  period 
t.  It  is  assumed  that  travel  time  between  any  two  cities  is  one  tune 
period.  Clearly  if  C  is  invariant  over  t,  the  1 1  >  1  SI'  reduces  to  the  1ST 
The  TUTSI’  formulation  of  (11)  can  he  expanded  to  define  the  tour 
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continuity  constraints  for  t he  I’llDP.  The  problem  is 
V  k 

minimize  )  c .  .  x .  . .  ( 3 .  8  | 
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X.  ;,  <  »  ,  k  =  1.2 . K 
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(3.12) 


x .  .  <.  f  0 , 1  }  for  a  1 1  i  ,  j ,  t  and  k 

i  jt 


(3.1 3 ) 


these  equations  arc  of  order  Nk,  while  the  conventional  expressions 
require  equations  of  order  N2 .  As  before,  the  number  of  decision 
variables  is  of  order  KN?.  repressions  (3.11)  and  (3.12)  arc  analogous 
to  (3.3)  and  (3.4).  The  contribution  of  (ox,  (lavish  and  braves  is  that 
expressions  (3.1>)  and  (3.10)  effect  suhtour  elimination. 
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IV.  ORIGIN/ DKST I  NAT I  ON  CONSTRAINTS 


The  key  element  that  distinguishes  the  pickup  and  delivery  problem 
from  the  TSP  and  the  VRP  is  the  precedence  requirement  of  origin  before 
destination.  This  requirement  can  be  expressed  as 
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for  ever>'  jeO  .  (3.14) 

Were  one  willing  to  include  the  decision  variable  t  in  the  formulation, 
then 


t.  <  x  for  every  jeO 

J  J+n 


(3.15) 


would  express  the  same  requirement.  The  latter  expression  does  not 
require  the  periodicity  parameter  t,  but  does  increase  the  number  of 
decision  variables. 

Unless  K  =  1,  it  is  also  necessary  to  insure  that  the  vehicle 
visiting  the  origin  is  the  same  vehicle  that  visits  the  correspondi ng 
destination.  This  can  be  accomplished  by 
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(3.10) 


V.  VFHICLli  CAPACITY  CONSTRAINTS 
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The  impact  of  vehicle  capacity  on  the  PUDP  is  totally  different 
than  it  is  on  the  VRP.  For  the  VRP,  a  given  set  of  customers  either  can 
or  can't  be  serviced  by  a  given  vehicle,  depending  on  the  sum  of  the 
customer's  requirements.  IVith  the  PUL)P(  it  is  the  sequence  of  stops 
which  determines  the  quantity  on  the  vehicle  at  any  time.  This  fact  is 
what  necessitates  the  use  of  the  period  parameter  t.  The  capacity  of 
a  vehicle  is  not  exceeded  at  any  time  so  long  as 
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VI.  T I  Ml:  W I NDOW  CONSTRA I  NTS 


There  are  two  approaches  to  defining  time  window  constraints.  One 
approach  employs  clock  time,  while  the  other  uses  the  periodicity 
parameter.  For  the  former,  the  expressions  are 


and 


h  i  H.  ,  i  =  2,3, . . . ,N  .  |  A .  IP) 

For  many  instances  of  the  problem,  it  is  possible  to  approximate 


clock  time  by  a  stop  number.  For  these  cases 
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l. 

J  N  K  , 

l  l  l  =  1  ,  J  =  2,3,... ,N  (3.20) 

t=e.  i=l  k=l  J 
J 

insures  that  the  time  window  constraints  are  included. 

VII.  QUALITY  or  si.rvici: 

As  above,  either  clock  time  or  the  number  of  periods  may  be  used 
to  represent  the  desired  standards  specified  by  Q..  Hence,  one  could 
use  either 
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or 
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to  insure  that  quality  of  service  is  maintained. 


VIII.  OPI: RATIONAL  CONSTRAINTS 


The  operational  constraint  may  be  expressed  by 


I  c .  .  x .  .  <  l),  ,  k  =  1 ,2 . 1 
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Altliough  the  above  formulation  does  in  fact  represent  the  i'UDP, 
the  necessary  quadruply  scripted  variables  should  be  noted.  The 
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The  implication  of  this  fact  will  become  apparent  in  the  next  chapt 
where  discussion  of  optimal  solution  techniques  will  be  discussed. 


CHAPTER  IV 


EXACT  SOLUTION  ALGORITHMS 

It  was  noted  in  Chapter  II  that  optimal  solutions  to  vehicle 
routing  problems  (VRR's)  were  difficult  to  come  by.  Because  the  pickup 
and  delivery  problem  (I’UOP)  is  more  complex  than  the  VRP,  one  might 
suspect  that  exact  solutions  to  it  would  be  even  more  elusive.  In  many 
cases  this  suspicion  is  correct.  However,  for  one  subclass  of  the 
PUDP,  large  problems  (over  100  stops)  are  readily  solved  optimally  using 
a  dynamic  programming  algorithm. 

Other  than  dynamic  programming,  the  exact  solution  techniques 
detailed  earlier  for  the  travelling  salesman  problem  (TSP)  and  the  VRP 
are  not  effective  in  solving  the  PUDP.  Not  all  PIJDP  instances  can  be 
solved  by  dynamic  programming.  Consequently,  this  chapter  is  devoted 
to  explaining  why  the  traditionally  successful  algorithms  fail,  or  are 
not  effective  when  applied  to  the  PUDP,  and  to  developing  the  conditions 
and  algorithms  that  allow  for  optimal  solution  to  relatively  large 
PUDP's.  A  detailed  discussion  of  the  multiple  vehicle  PUDP  is  deferred 
until  Chapter  VI. 


I.  INLITLCTJVE  1ECHN1QUES 

Integer  programming  by  means  of  cutting  planes,  branch  and  bound 
approaches,  and  Lagrangean  techniques,  has  been  used  with  varying 
degrees  of  success  for  solving  the  TSP  and  the  VRP.  In  addition,  a 
Benders  decomposition  approach  has  been  successfully  applied  to  the  VRP 

4(» 
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as  both  a  heuristic  and  an  exact  solution  technique.  The  complexity 
of  the  constraints  renders  each  of  these  techniques  ineffective  when 
applied  to  the  PUDP. 


4? 


Integer  Programming 

The  integer  programming  (IP)  formulations  of  the  PUDP  given  in 
Chapter  III  are  considered  as  compact  as  possible.  Although  the  compact 
formulation  is  only  of  order  KN,  the  number  of  decisions  variables 
presents  the  primary  difficulty.  At  least  K.N 3  (0,1)  decision  variables 
are  required. 

Consider  a  hypothetical  problem  with  15  customers,  or  31  stops. 
Nearly  30,000  decision  variables  are  required  for  the  single  vehicle 
version  of  the  problem  and  more  than  100,000  variables  are  needed  if 
the  customers  are  to  be  serviced  by  four  delivery  vehicles.  For  the 
TSP ,  integer  programming  techniques  such  as  cutting  planes  and  implicit 
enumeration  generally  have  performed  less  efficiently  than  have  others. 
Given  the  large  number  of  decision  variables  and  tiic  suspect  efficiency 
of  cutting  planes  in  general,  further  pursuit  of  IP  was  not  undertaken. 

Branch  and  Bound 

Branch  and  bound  techniques  are  by  far  the  most  common  type 
applied  to  the  TSP  and  to  the  VRP.  This  is  true  because  relaxations  of 
the  TSP  or  the  VRP,  such  as  t Ire  assignment  problem,  arc  readily 
solvauic  and  provide  logical  branching  points.  Also,  good  lower  bounds 
on  the  optimal  completion  of  the  subproblem  can  be  readily  computed. 
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Tight  bounds  and  logical  branching  strategies  are  not  readily 
available  for  the  PUDP.  The  precedence  requirement  is  the  primary 
complicating  factor.  The  sequence  in  which  stops  are  visited  is 
irrelevant  in  the  TSP,  but  critical  with  the  PUDP.  Consequently,  the 
relaxation  solution  may  produce  infeasible  subtours.  Therefore,  one 
must  not  only  branch  to  eliminate  subtours,  but  also  to  achieve 
feasibility.  Given  both  subtours  and  infeasibility,  the  choice  of  a 
branching  rule  is  not  readily  apparent,  liven  if  such  a  rule  were 
readily  available,  the  need  to  compute  tight  bounds  remains  as  a 
necessary  criterion  for  a  good  branch  and  bound  algorithm.  Such 
bounds  are  not  readily  available. 

The  PIJUP  is  a  restriction  of  the  TSP.  Therefore,  the  optimal 
solution  to  the  ISP  is  a  lower  bound  on  the  PUDP.  lor  noncontrived 
problems,  we  will  see  that  the  TSP  provides  a  very  poor  lower  bound. 
Consequently,  to  efficiently  solve  the  PUUP  by  branch  and  bound,  a 
relaxation  of  the  PUUP  is  needed  that  produces  bounds  better  than  the 
optimal  TSP  solution  over  the  same  network,  further,  tie  bound  must 
be  easily  computed  in  polynomial  time  (.such  a  relaxation  is  referred 
to  as  a  "polynomial  time  relaxation") .  The  following  lemma  shows  that 
such  a  relaxation  probably  does  not  exist. 

Lemma  2 .  Unless  P  =  NP,  there  docs  not  exist  a  polynomial  time 
relaxation  of  the  PUUP  witli  value  R  which  is  greater  than  or  equal  to 
the  optimal  ISP  relaxation  solution  value  of  Z*  for  all  PUUP  instances. 
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Proof.  Let  P  be  the  value  of  the  optimal  Plll'P  solution  and 


suppose  Z*  <  R  <  P"  for  every  instance  of  the  Plll'P.  i.et  A  be  the 

polynomial  time  algorithm  that  produces  R. 

Consider  the  following  instance  of  the  Plll'P:  one  vehicle,  no 

capacity,  operational,  time  window,  or  quality  of  service  constraint 

other  than  for  stop  n  and  a  requirement  that  all  origins  be  visited 

before  time  T  and  all  destinations  after  time  T.  Origin  n  and  its 

corresponding  destination  dn  are  defined  such  that  c  -  i.  =  T  and 

n  n 


=  =  T  ♦  1,  c 


n,j 


.  for  all  j,  c 


1  ,  n 


1  , 


for  all  i  and 


cfi  7n  =  0.  In  effect,  stop  n  and  2n  are  identical  to  the  depot. 

As  defined,  the  PUl'P  reduces  to  two  PSP's.  One  over  the  n-1 

origins  and  one  over  the  n-1  destinations.  I.et  c..  =  0  for  all  i, 

jrD,  =  »  for  all  if.O  and  jrl>.  Apply  A  to  the  resulting  i'UDP. 

Clearly  P*  is  the  value  of  the  optimal  'ISP  solution  over  the  depot 

and  n-1  origins,  which  is  also  the  value  of  the  Plll'P  solution.  Hence, 

★  —  * 

Z  =  R  =  P  from  the  hypothesis.  But  any  TSP  could  easily  be  trans¬ 
formed  to  such  a  PUDP  instance.  Ihcrefore,  A  can  be  used  to  solve  any 
TSP.  Since  A  is  a  polynomial  time  algorithm,  this  would  imply  I'  =  NP. 
Consequently,  we  rejected  Z*  <_  R  <  P*  and  accept  the  provisions  of  the 
lemma. 


empirical  evidence  suggests  that  obtaining  a  bound  better  than  the 
TSP  bound  is  difficult  foi  any  instance  of  the  PUl'P.  No  workable 
bounding  scheme  was  uncovered  that  was  not  also  a  lower  bound  for  the 
TSP.  The  LP  solution  obtained  by  relaxing  the  integrality  requirement 


could  he  such  a  bound,  but  is  not  easily  computable  due  primarily  to  the 
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previously  noted  number  of  decision  variables  required.  Given  such 
loose  bounds,  any  branching  tree  would  clearly  become  enormous. 

Attempted  hand  solution  of  a  five  stop  (two  origin/destination  pair) 
problem  demonstrated  the  futility  of  the  branch  and  bound  technique. 
Therefore,  it  is  not  deemed  practical  for  optimally  solving  the  general 
PUDP. 

Lagrangcan  Technique 

One  of  the  most  powerful  techniques  for  solving  the  symmetric  T.S1' 
was  i —  1-tree  approach  of  Held  and  Karp  (22,23).  The  precedence 
requirement  of  the  PUIH1  can  not  be  accomodated  by  the  1-tree  approach. 
The  construction  of  the  1-tree  requires  the  construction  of  a  minimum 
spanning  tree.  The  greedy  algorithms  used  for  the  solution  of  the 
minimum  spanning  tree  problem  requires  an  independence  relationship  that 
is  immediately  contradicted  by  the  "origin  before  destination" 
constraint . 

The  resuit  is  similar  to  that  cited  above.  Whereas  the  ISP  is 
sequence  independent,  the  PIJDP  is  sequence  dependent.  This  dependence 
is  the  primary  factor  precluding  or  limiting  the  solution  of  the  PtIDP 
by  the  1-trec,  branch  and  bound,  and  integer  programming  approaches. 

Benders  Decomposition 

Ihe  Pisher-.Iaikamur  algorithm  is  based  on  iterating  between  a 
generalized  assignment  problem  (GAP)  and  a  travelling  salesman  sub- 
problem.  The  objective  of  the  GAP  is  to  optimally  assign  customers  to 
vehicles  such  that  vehicle  capacity  is  not  exceeded. 


t 


With  the  PUDP,  vehicle  capacity  may  not  even  he  a  factor.  liven  if 
it  were,  the  capacity  limitation  affects  the  problem  differently.  It 
is  dependent  upon  the  specific  routes  followed.  Remaining  unused 
vehicle  capacity  decreases  as  an  origin  is  visited  and  increases  as  a 
destination  is  visited.  Thus,  for  a  given  combination  of  customers,  one 
or  more  sequences  might  cause  the  capacity  limitation  to  be  violated, 
while  another  would  be  feasible.  This  can  be  contrasted  to  the  VRP 
where  one  normally  considers  the  vehicle  fully  loaded  when  it  departs 
the  depot.  In  this  varying  capacity  environment,  the  GAP  is  .v  t 
clearly  defined.  Assigning  customers  to  the  correct  vehicle  is  the 
crux  of  the  multiple  vehicle  PIJDP.  However,  no  logical  formulation  of 
the  customer  assignments  could  be  found  so  that  the  Renders  decomposi¬ 
tion  concept  could  be  applied  to  obtain  an  exact  solution.  Chapter  VI 
specifically  addresses  the  multiple  vehicle  PUDP. 

II.  DYNAMIC  PROGRAMMING  AND  PRIit  hi'I.NCi. 

CONSTRA I  NTS 

Dynamic  programming  has  been  used  less  on  the  ISP  than  some  other 
techniques.  The  computer  storage  requirements,  as  noted  in  Chapter  If, 
are  the  primary  difficulty  in  using  DP  to  solve  the  TSP.  Yet  dynamic 
programming  can  be  adapted  to  handle  the  precedence  relationship  of 
the  PUDP.  further,  this  can  be  accomplished  with  a  significant 
reduction  in  storage  requirements. 

The  reason  for  the  reduction  is  that  infeasible  states  are  not 
generated,  for  a  la  node  TSP,  nearly  25,000  individual  storage  location 


are  required  for  the  data  generated  by  the  recursion  equations.  For 
the  PUDP  of  equal  size  (six  customers  plus  the  depot)  and  only 
precedence  constraints  binding,  less  than  3,000  locations  are  needed. 

The  next  sections  explain  how  this  reduction  in  storage  requirements  is 
possible.  Further  reductions  are  possible  for  some  more  heavily  con- 
strainted  instances  of  the  problem.  These  instances  arc  discussed  later 
in  the  chapter. 

Induction  Schemes 

In  (2.17),  the  basic  dynamic  programming  recursion  for  the 
travelling  salesman  problem  was  given.  The  expression  for  the  shortest 
partial  tour  from  node  1  to  node  j  that  passed  through 

Sk  =  {i  i,  i2, .  .  . ,  ik  l }  (i.l) 


was 


fk(jlsk> 


=  m  m  i  mum  ( ,  ,  (  i  s.  -  i  1  +  c . 

1  k  - 1  m  k  m  i,i 

in  m  • 


In  this  form,  the  recursion  can  be  thought  of  as  forward  induction  in 
that  the  salesman  is  proceeding  forward  from  the  starting  point  to  some 
intermediate  point  in  h i s  tour.  Often,  dynamic  programming  is  thought 
of  as  backward  induction. 

In  backward  induction,  the  ri  -ursinn  -•ivU  to  Find  the  shortest 
partial  tour  from  node  i  that  passes  through 


Sk  ~  f  1  N  +  2-k  '  'n+1  -k . V 


(1.3 


and  returns  to  the  original  origin,  node  1.  The  subscr i pt  k  indicates 
the  total  number  of  nodes,  including  node  i  tiiat  are  included  in  the 
partial  tour  not  including  the  depot.  As  developed  in  Chapter  III, 


N  is  the  total  number  of  distinct  nodes.  Mathematically,  the  recursive 
expressions  are  identical.  'I he  only  difference  is  in  the  definition  of 
S.  .  Backward  induction  is  used  in  the  algorithms  deve'oped  to  solve 
the  J'UDP.  Therefore ,  all  explanations  to  follow  assume  the  backward 
oricntat  ion. 

State  Space  Generation 

Only  feasible  state  vectors  need  to  bo  generated  for  the  PHI '1' . 

Stop  j  of  expression  11. 2)  will  he  referred  to  as  the  lead.  Conse¬ 
quently,  if  j  is  the  lead  and  j  is  one  <>t  the  n  origins,  it  must  be 
true  that 

(j+nl  sk  |  i 

if  the  vector  is  to  be  feasible.  Sinn l.u  1 /.  if  i  i s  a  d. st  mat  u-n, 
it  must  be  true  that 

< J-nj  t  Sk  ,  j 

Also,  if  stop  i  is  an  origin,  then 
i  t  Sk  =>  ( i +n)  ‘  Sk 
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State  Vector  kepresentat  ion 

A  binary  vector  representing  the  lead  unioned  with  an  appropriate 
S,  is  generated  and  stored  for  each  f  ^  value.  these  vectors  are  used 
to  identify  each  Sj.  and  to  identify  the  optimal  solution  sequence. 

Each  vector  takes  advantage  of  t ho  way  in  which  decimal  numbers  are 
stored  by  a  digital  computer.  If  i  represents  any  one  of  the  N 
different  stop  numbers,  then  the  binary  representation  of  the  decimal 
number  2l  1  consists  of  one  1  in  the  appropriate  location  and  0’s  in 
all  other  positions  of  the  word.  Suppose  JUS  =  <'2,4,6}  and  the  computer 
system  employs  an  eight  bit  word,  figure  4  depicts  the  binary  repre¬ 
sentation  of  the  decimal  number  12,  which  is  obtained  by 

22'1  ♦  2*" 1  *  26' 1  =  2  +  S  *  52  =  42  .  (4.7) 

Consequent ly ,  the  decimal  number  42  corresponds  to  the  presence  of  stops 
number  2,  4  and  6  in  the  state  vector.  It  is  easily  seen  that  all 
possible  combinations  of  stops  have  unique  binary  representations. 

lor  problems  where  N  is  greater  than  tin-  length.,  !,  of  an  integer 
word  on  the  computer  system  in  use,  the  representation  is  separated 
into  two  (or  more)  vectors  such  that  each  word  represents  part  of  the 
total  state  vector.  for  example,  suppose  \  -  45  and  1  =  52.  Then  for 
given  state  * ,  two  words,  RIl'l  1*1  and  Kl.1'2  1*1,  could  he  used  to 
indicate  the  state  condition  for  f ^ ( * ) .  Stop  numbers  1  through  51  could 
he  tracked  by  Ul.l’l  (*l  and  52  through  45  by  kll’2  l‘  t.  Jhe  concept  is 
similarly  extended  to  Kll’5,  etc,  as  necessary  for  large’’  problems. 


Storage  Requirements 


As  noted,  the  number  of  storage  locations  required  to  solve  the 
PUDP  is  less  than  that  for  the  TSP.  The  feasibility  criteria,  expressed 
by  (4.4)  through  (4.6),  decrease  the  number  of  combinations  that  must 
be  considered  for  each  lead.  This  difference  can  best  be  visualized  by 
computing  the  storage  locations  required  for  one  step  of  a  sample  TSP 
and  for  the  related  PUDP. 

lixampl c  requirements.  Let  N  =  IS  and  k  =  7.  Suppose  the  lead  is 
stop  number  2,  which  corresponds  to  an  origin  in  the  PUDP  with  the 
corresponding  destination  being  stop  number  9.  l-'or  the  TSP,  with  the 
depot  fixed  as  stop  number  1,  there  arc  f1^)  =  1,716  different  S7 
values  to  be  computed,  for  the  PUDP,  stop  number  9  must  be  included 
in  all  S?  values.  Thus,  there  arc  five  free  positions  available,  which 
can  be  filled  with 

•  2  origin/destination  pairs  and  1  destination, 

•  1  origin/destination  pair  and  3  destinations,  or 

•  5  destinations. 

The  number  of  combinations  for  each  configuration  is 

•  f^)  rj)  =  IS  •  4  =  60, 

•  ((’l  (*)  =  6  ■  10  =  60, 

•  J  =  6, 

respectively,  for  a  total  of  only  126  values  to  he  computed.  This 
represents  more  than  a  tenfold  decrease  in  storage  requirements.  Were 


the  lead  position  a  destination,  141  storage  locations  would  be 
required.  For  a  problem  with  N  =  29,  over  a  hundredfold  decrease  can 
be  achieved.  However,  the  growth  of  storage  locations  remains 
exponential . 

Total  requirements.  An  upper  limit  on  the  total  number  of  storage 
locations  required  of  2n  •  3n_1  +  1  is  given  by  Psaraftis  (24),  where 
n  is  the  number  of  origin/destination  pairs.  His  development  of  the 
states  is  different,  but  requires  essentially  the  same  number  of  total 
states  for  the  PUDP  with  only  precedence  constraints  binding.  The 
inclusion  of  time  window  and/or  quality  of  service  constraints  for  some 
problem  instances  reduces  storage  requirements  even  further,  in 
contrast  to  Psaraftis'  development,  as  will  be  seen  later  in  this 
chapter. 

Lxamplc  Problem 

An  example  problem  will  be  used  to  explain  exactly  how  t lie  dynamic 
programming  algorithm  developed  solves  t lie  PUDP  with  only  precedence 
constraints  binding.  A  problem  with  three  customers  or  origin/desti¬ 
nation  pairs  will  be  assumed.  The  depot  is  node  1,  while  the  origin/ 
destination  pairs  arc  2/5,  3/b,  and  4/7,  respectively.  Table  1  depicts 
the  cost  matrix  for  the  hypothetical  problem.  Inconsistent  cells, 
such  as  a  destination  to  its  correspond ing  origin,  are  blacked  out  to 
so  indicate.  However,  any  cost  could  be  placed  in  one  of  these  cells, 
as  the  algorithm  does  not  allow  for  such  i nfeas i h i 1 i t i es .  Table  2 
presents  all  of  the  f ^ ( j ] )  values,  in  the  order  generated,  as  well 
as  the  binary  identification  vectors. 
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tabll:  2 

RliCURSIVli  DATA  TOR  l.NAMI'U  DYNAMIC 
PROGRAMMING  PROBLLM 


fk 

COUNT!.  R 

BINARY 

Rl 

PR1.SLYI  AT  ION 

VALUII 

LOAD 

VliCTOR 

8 

-7 

-6- 

5 

4- 

3- 

T 

-I 

8 

2 

2  1  5 

0 

0 

0 

1 

0 

0 

1 

0 

10 

3 

3  1  6 

0 

0 

1 

0 

0 

1 

0 

<) 

1 

4 

4  1  7 

0 

1 

0 

0 

1 

0 

0 

(1 

9 

5 

5  1 6 

0 

0 

1 

1 

0 

0 

0 

0 

9 

5 

5  1  7 

0 

1 

0 

1 

0 

0 

0 

0 

7 

6 

615 

0 

0 

1 

1 

0 

0 

0 

0 

1 

6 

617 

0 

1 

1 

0 

0 

0 

0 

n 

8 

7 

715 

0 

1 

0 

1 

0 

0 

0 

(i 

2 

7 

7  1 6 

0 

1 

1 

0 

0 

0 

0 

(1 

7  2  2 1 5(, 

11  2  2157 

9  5  3105 

9  3  3 1 67 

8  4  4175 

2  4  4 1 76 

10  5  5 1 36 

9  5  5147 

8  5  5167 

8  6  61 25 

9  6  6147 

8  6  6  15" 

11  7  7 1 25 

21  7  7136 

7  7  7156 


0  0  1  1  0  0  1  U 

0  10  1  0  0  1  II 

0  0  110  1  0  0 
(1110  0  1  0  0 
0  10  110  0  0 
0  110  10  0  0 
0  10  11  0  0  0 
0  10  110  0  0 
0  111  0  0  0  n 
0  0  1  ]  (I  0  1  0 

0  1  10  10  O  (I 

0  111  0  0  0  o 
0  1  0  0  0  0  1  0 

0  1  10  0  1  O  0 

0  1110  0  0  0 


14  2  2 1 536 

17  2  21 54  7 

8  2  21 567 

12  3  31625 

I  3  3  3  1 64  7 

8  3  31657 

II  4  41725 

15  4  41736 

4  1756 


0  0  1  10  1  10 
0  1  0  1  1  0  1  () 
0  1110  0  1(1 
0  0  110  110 
0  1  10  1  10  0 
0  1110  10  0 
0  10  110  10 
oiioi  ioo 
oil  11000 
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TABU:  2  (continued) 


fk 

VALUE 

LEAD 

COUNTER 

VECTOR 

B I  NARY 
8-7- 

REPRESENTATION 
6-5-4- 3-2- 1 

k  =  4 

9 

5 

51367 

0 

1 

1 

1 

0 

1 

0 

0 

10 

5 

51476 

0 

1 

1 

1 

1 

0 

0 

0 

1 1 

6 

6  1257 

0 

1 

1 

1 

0 

0 

1 

0 

16 

6 

61475 

0 

1 

1 

1 

1 

0 

0 

0 

8 

7 

7l2ro 

0 

1 

1 

1 

0 

0 

1 

0 

18 

7 

7  |  365 

0 

1 

1 

1 

0 

1 

0 

0 

k=5 

13 

2 

215367 

0 

1 

1 

1 

0 

1 

1 

0 

16 

2 

215476 

0 

1 

1 

1 

1 

0 

1 

0 

13 

3 

316257 

0 

1 

1 

1 

0 

1 

1 

0 

10 

5 

3  1 64  75 

0 

1 

1 

1 

1 

1 

0 

0 

8 

4 

4 |7256 

0 

1 

1 

1 

1 

0 

1 

0 

12 

4 

417365 

0 

1 

1 

1 

1 

1 

0 

0 

13 

5 

513647 

0 

I 

1 

I 

1 

1 

0 

0 

17 

6 

612547 

0 

1 

1 

1 

1 

0 

1 

0 

17 

7 

712536 

0 

1 

1 

1 

0 

1 

1 

0 

k  =  b 

15 

2 

2153647 

0 

1 

1 

1 

1 

1 

1 

0 

19 

3 

3162547 

0 

1 

1 

1 

1 

1 

1 

0 

17 

4 

4  172536 

0 

1 

1 

1 

1 

1 

1 

0 

k=7 

17 

1 

1  1255647 

0 

1 

1 

1 

1 

1 

1 

1 

Optimal  Solutions: 


1-2-5-5-4-7-6-1 
1  -4-7-2-3-5-6-1 


and 


The  algorithm  used  actually  generates  a  vector  of  length  k  that  is 
incremented  to  obtain  each  feasible  state,  like  a  counter,  bach  element 
of  the  k-vector  is  then  used  to  compute  the  binary  representation,  which 
is  stored  for  each  f^  value.  If  the  lead  position  is  an  origin,  the 
next  position  in  the  k-vector  is  automatically  set  to  that  origin's 
corresponding  destination.  Consequently,  only  k-2  positions  would 
require  incrementing,  further,  any  time  one  of  the  elements  incremented 
is  an  origin,  the  next  element  is  similarly  set  to  its  correspond i ng 
destination.  Thus,  for  k  =  5  the  first  counter  vector  is  given  in 
Table  2  as 

2  |  5367. 

Origin  2  as  the  lead  automatically  sets  the  next  position  to  2's 
destination,  5.  Origin  3  automatically  sets  the  position  on  its  right 
to  its  destination,  6.  bach  time  one  element  of  the  k-vector  is 
incremented,  all  elements  to  the  right  of  it  are  initialized  to  their 
lowest  feasible  node  number.  Then  the  vector  is  incremented  like  a 
counter,  with  the  right  most  clement  incremented  first.  Once  completed, 
each  element  of  the  k-vcctor  is  then  used  to  compute  the  binary  represen 
tation,  which  is  stored  for  each  f ^  value.  The  binary  representation 
conveys  only  the  identity  of  the  nodes,  not  the  sequence  in  which  thee 
were  generated. 

f 2 ( j I i M j •  The  algorithm  begins  by  computing  f ? (j | i  ) .  Recalling 
that  the  algorithm  uses  backward  recursion,  there  are  only  two  configtira 
tion  patterns  possible: 


*  an  origin/destination  pair  with  return  to  the  depot,  and 

•  two  destinations  with  return  to  the  depot. 

In  the  latter  case,  any  destination  could  serve  as  the  lead,  while  in 
the  former,  only  an  origin  could.  The  computation  in  either  case  is 
straightforward  and  will  not  be  elaborated  on  further. 

f 3 ( j I S 3 ) •  Beginning  with  j,  or  the  lead,  equal  to  the  first  origin, 
node  2,  wo  generate  all  possible  S3  combinations.  Because  node  2  is  an 
origin,  node  5  must  be  included,  leaving  only  one  free  position  in  S3. 
Clearly,  this  node  must  be  a  destination.  Beginning  with  the  lowest 
numbered  destination  node,  not  identical  to  any  destination  already 
fixed,  each  state  vector  is  incrementally  generated,  expression  (4.2J 
and  the  previously  generated  data  arc  then  used  to  compute  each  f3 
value.  Once  all  S3  combinations  arc  generated  for  node  2,  the  lead  is 
incremented  and  the  procedure  repeated. 

When  the  lead  is  incremented  to  node  number  5,  the  first  destination 
node,  there  are  two  free  positions  in  S-,  that  can  be  filled  In’: 

•  an  origin/destination  pair,  or 

*  the  remaining  two  destinations. 

In  the  former  case,  an  origin  requires  its  corresponding  destination. 
Thus,  when  the  origin  position  in  the  state  vector  is  being  incremented, 
there  is  no  free  position  left  in  S3.  lor  the  latter,  one  free  position 
rema ins. 

f 4 ( j I Sh ) .  The  procedure  for  f„  values  is  the  logical  extension  of 
that  for  f‘3  values.  With  an  origin  as  a  lea.l,  the  two  free  positions 
in  Si.  can  he  filled  by  either: 


(>3 

•  an  origin/destination  pair,  or 

•  two  destinations. 

Similarly,  when  the  lead  position  is  a  destination,  an  origin/destination 
pair  with  one  destination  is  tlte  only  possible  combination.  A 
destination  followed  by  three  destinations  is  impossible  since  there  are 
only  three  total  destinations  in  tlte  problem. 

Remaining  fk(j|Sk)'s.  The  same  process  is  repeated  for  f *  and  ff 
values.  Tor  ft,  values,  only  origins  are  legitimate  lead  values, 
followed  by  all  other  node  numbers.  The  optimal  solution  value  is 
obtained  by  solving 

f7(l  2,3,...,  7)  =  min  [feOjSc-i^  *  c.  )  .  1 4 .  S ) 

m  m , 

Because  only  origins  can  be  in  the  lead  position  in  f, ,  one  is  assured 
that  i  tO  and  that  the  route  begins  legitimately.  One  can  then  identify 
an  optimal  solution  with  value  1".  One  such  solution  is  1-2-.'-.'-  I  ~  r.]. 

Obtaining  the  solution.  Using  the  f^  values,  the  lead  node  lot 
that  value,  the  binary  representation  and  the  original  cost  data,  tin- 
solution  can  be  constructed  one  node  at  a  time.  lor  example,  the  f  • 
value  of  17  minus  the  cost  from  the  depot  to  node  2,  which  is  two, 
equals  the  f&  value  of  15  associated  with  node  2  as  a  lead.  Similarly, 
this  value  of  15  minus  the  cost  from  node  2  to  3  of  5  equals  the  f, 


value  of  10  associated  with  node  3  as  the  lead.  This  assignment  is 
legitimate  since  the  identification  vectors  indicate  that  no  previously 
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scheduled  stops  (in  this  case  node  2)  were  used  to  obtain  the  value. 
The  process  continues  until  all  stops  are  sequenced.  The  tour  is  then 
completed  by  indicating  the  return  to  the  depot. 

Computational  Results 

With  only  the  precedence  constraints  binding,  the  exponential 
growth  of  the  required  storage  locations  limits  the  sice  of  problem 
that  can  be  practically  solved.  Table  3  summarizes  results  for  four 
small  problems.  All  computer  runs  were  made  on  the  University  of 
Tennessee  Computing  Center's  IBM  370/3031  computer  system.  Cost  matrix 
data  were  generated  randomly.  Both  x  and  y  graph  coordinates  between 
0  and  99  were  generated  for  each  point.  Interpoint  distances  were  then 
computed  using  the  rectilinear  metric,  i.c. 

c(i,j)  =  1^  -  Xj  |  +  l>'i  -  >j  I  (4.9) 

where  x.  is  the  x  coordinate  of  the  ith  point,  etc.  The  data  generation 
technique  was  the  same  for  all  problems  examined  throughout  this 
dissertation.  All  of  the  algorithms  developed  in  tiiis  and  the  next  two 
chapters  are  independent  of  the  manner  in  which  the  cost  matrix  is 
obtained.  The  rectilinear  metric  was  selected  only  bceuase  it 
facilitated  hand  computations  during  the  early  stages  of  algorithm 
development.  Any  metric  could  be  used  in  place  of  the  rectilinear  one. 
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as  well  as  the  solution  procedure. 
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The  discussion  in  the  last  section  dealt  only  with  binding 
precedence  constraints.  In  some  real  world  applications  of  the  pickup 
and  delivery  problem,  many  other  constraints,  as  suggested  earlier, 
may  be  present.  In  Chapter  III,  two  different  approaches  for  formulating 
the  fully  constrained  version  of  the  PUMP  were  presented.  One  employed 
the  use  of  clock  time,  while  the  other  used  the  stop  number  of  a 
customer  in  the  optimal  sequence.  This  latter  approach  is  ideally 
suited  to  solution  by  dynamic  programming. 

Lead  _1  *os  it  ion  Leas  i  bi  lit  v 

Position  identification.  The  first  origin  visited  after  the 
vehicle  leaves  the  depot  is  considered  to  he  the  first  stop  number. 
Consequently,  the  last  destination  the  vehicle  visits  before  returning 
to  the  depot  is  designated  stop  number  N - 1 ,  with  all  other  origins  and 
destinations  being  sequenced  somewhere  in  between.  l.aiTier  in  this 
chapter  it  was  noted  that  the  subscript  k  on  the  indicated  the 
number  of  nodes  or  customers  in  the  partial  tour  from  node  i  to  the 
depot.  Let  p  represent  the  stop  number  associated  with  k.  Then  the 
stop  number  can  he  found  by 

"k  "  N- k  ,  k  -  d ,  a . \-l 


I  1  .  1  0  I 


Leas ibi 1 i t  v  sot . 


(1" 

The  dynamic  prog rami'! i ay  algorithm  previously 
described  proceeded  b>'  first  fixing  the  lead  customer,  and  then 
generating  all  combinations  of'  subsequent  stops.  However,  when  con¬ 
straints  other  than  precedence  requirements  are  present,  it  is  likely 
that  not  all  customers  can  feasibly  b<  assigned  to  a  given  "p  lead 
position.  Therefore,  there  is  no  need  to  generate  all  subsequent 
customer  combinations,  nor  to  store  any  data  relative  to  infeasible 
combinations.  Let 


.1 


k 


{ j  !  j  is  a  feasible  lead  for  ” [  : 
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State  vectors  will  be  generated  onlv  for  if. I 

k 

Time  windows.  In  order  for  the  complement  of  d  to  not  be  er.pty 

for  all  7‘p,  at  least  one  of  the  customer  requirements  must  have  a  not 

earlier  than  (NLTJ  requirement  or  a  not  later  than  ( N 1.  i  requirement, 

or  both .  The  more  interesting  case,  in  terms  of  solution  by  dynamic 

programming,  is  wtien  each  customer's  requirements  involve  time  windows, 

which  it  is  assumed  can  be  defined  by  stop  numbers.  When  tire  windows 

are  present  for  all  customer  requirements,  the  set  .1,  is  clearly 

k 

defined  for  all  ~ ^ ,  and  the  total  number  of  storage  locations  required 
for  the  state  and  binary  representation  vector  -  reduced.  i  quail) 
significant  is  that  main  of  flic  other  possible  con  "tr. lints  at  e  a  1  so 


definable  in  terms  of  tine  windows. 
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Handling  Other  C.'onst  ra  ints 

The  presence  of  N!:T  or  NI.'I  times  will  he  assumed  for  at  least  one 
portion  (either  the  origin  or  the  destination)  of  a  customer's  require¬ 
ment.  One  point  in  time  is  all  that  is  necessary  to  fix  delivery 
windows,  as  explained  below,  for  servicing  both  the  customer's  origin 
and  destination.  In  many  cases,  this  one  point  would  be  a  NLT  time 
for  delivery  to  the  destination.  for  ease  of  explanation,  the  dial-a- 
ridc  problem  will  be  used.  However,  any  l'DHP  for  which  the  appropriate 
assumptions  hold  could  He  substituted  equally  as  well. 

Quality  of  service.  Let  Q  represent  the  quality  of  service  for 
the  dial-n-ridc  service  problem.  Q  is  specified  by  the  provider  based 
on  his  standards.  It  represents  the  maximum  number  of  stops  before  a 
given  customer  reaches  his  destination.  Conversely,  given  a  required 
delivery  time,  represented  by  a  specific  stop  number,  that  time  minus 
Q  represents  the  earliest  time  that  a  customer  could  be  picked  up  at 
his  origin.  T lie  use  of  stop  numbers  to  represent  time  is  absolutely 
necessary.  Consequently,  the  direct  travel  time  between  two  points 
is  taken  to  be  one  stop,  regardless  of  the  actual  clock  time  required. 
Supp  isc  a  customer  must  be  at  his  destination  by  stop  1.  In  this  case, 
the  Nl.T  for  the  customer1, •  destination  is  I.  Ihe  Nl.l  time  for  the 
corresponding  origin  is  1-1  since  it  takes  a  minimum  of  one  period  to 
travel  from  a  customer's  origin  to  his  destination.  Since  Q  is  the 
maximum  difference  between  when  a  customer  must  be  at  his  destination 


and  when  a  customer  is  actually  picked  up,  the  earliest  time  that  the 


customer  could  he  picked  up  would  be  T-Q.  Thus,  T-Q  represents  the  SILT 
time  for  the  customer's  origin.  A  NUT  value  of  T-Q+l  for  the  customer' 
destination  is  obtained  since  there  is  a  minimum  of  a  one  stop  travel 
time.  Suppose  Q=5  and  T=1S  for  a  given  customer.  Then  the  Mil  and  NI.T 
times  for  that  customer's  destination  would  be  18-5  1=14  and  18.  lor 
the  customer's  origin,  these  times  would  be  18-5=15  and  18-1=17, 
respectively.  The  overall  time  window  associated  with  this  customer 
would  thus  be  stop  numbers  13  through  18,  inclusive.  In  addition, 
there  would  be  a  window  of  13  through  17  inclusive  for  his  origin  ai.d 
14-18  inclusive  for  his  destination.  Any  attempt  to  provide  servi  v  at 
stop  numbers  outside  these  windows  would  result  in  an  infeasible 
solut ion . 

It  one  wishes  to  avoid  long  waiting  periods  upon  arriving  at  his 
destination,  a  maximum  waiting  parameter,  M,  could  be  used  so  that  the 
M:T  time  for  arrival  at  the  destination  could  be  compute d  as  T-M. 
Continuing  the  numerical  example  of  the  previous  paragraph,  if  M=3,  the 
earliest  time  (stop  number)  that  the  customer  could  arrive  at  his 
destination  would  be  18-3=13.  Ihe  overall  and  origin  time  windows 
would  remain  unchanged.  However,  the  destination  time  window  would 
become  stops  13  through  18  inclusive.  Both  Q  and  M  were  used  for 
computational  experimentation,  which  will  be  presented  later  in  this 
chapt er . 

Capacity  constraints.  lor  the  dial-a-ride  problem,  the  capacity 
constraint  is  simpler  because  each  person  can  be  thought  of  as  one  unit 
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Therefore,  one  must  only  be  concerned  that  the  number  of  people  in  the 
vehicle  at  any  time  does  not  exceed  the  vehicle’s  capacity.  This  can 
be  accommodated  by  using  a  not  unrealistic  assumption  that  the  vehicle 
capacity  is  greater  than  or  equal  to  Q.  Since  Q  represents  the  maximum 
number  of  stops  before  a  customer  reaches  his  destination,  and  since 
each  customer  represents  a  unique  stop,  the  vehicle  epapeity  can  never 
be  exceeded. 

The  assertion  that  the  needed  assumption  is  not  unrealistic  can 
be  argued  in  terms  of  a  hypothetical  situation.  Suppose  our  dial-a- 
ride  service  provides  transportation  services  to  handicapped  people. 

The  van  to  be  used  can  hold  at  most  seven  passengers,  fifteen  requests 
for  service  have  been  received,  bach  request  corresponds  to  a  pickup 
at  an  origin  and  delivery  to  a  destination.  Hence,  our  driver  has  50 
specific  stops  to  make.  Given  that  the  passengers  are  handicapped  and 
require  assistance  entering  and  leaving  the  vehicle,  we  suppose  that 
the  average  time  between  any  two  stops  in  an  optimal  sequence  is 
15  minutes.  When  one  physical  location  corresponds  to  more  than  one 
stop,  the  15  minute  average  should  not  be  materially  affected  by  the 
increased  time  required  to  provide  assistance  to  more  than  one  customer. 
These  numbers  arc  arbitrary,  but  do  lead  to  a  convenient  eight  hour 
workday  for  the  driver. 

The  vehicle's  capacity  could  not  be  exceeded  unless  at  least  one 
person  were  required  to  ride  for  two  or  more  hours  to  arrive  at  his 
destination.  Note  that  such  an  eventuality  is  necessary  but  not 
necessarily  sufficient  for  vehicle  capacity  to  be  exceeded.  Two  hours 


tamm 


"i 

appears  excessive  in  terms  of  the  quality  of  service  that  any 
organization  providing  dial-a-ride  service  would  provide.  Consequently, 
the  assumption  is  judged  not  unrealistic. 


Storage  Requirements 


Reduction  based  on  lead.  By  using  the  set  .J  ,  as  defined  by 
expression  (4.11),  the  storage  requirements  for  the  state  and  repre¬ 
sentation  vectors  can  he  reduced,  for  example,  if  for  each  value  of 
7t ,  only  30°o  of  the  j  arc  included  in  J^,  then  the  total  storage 
required  is  only  SO'i  of  that  required  for  the  unconstrained  instance. 
This  figure  makes  no  assumptions  about  the  feasibility  of  each  S^. 
further  reductions  are  possible  when  this  is  considered. 


feasible  elements  in  Sg.  'lhc  sets  Sg  represent  the  customers  who 
are  served  subsequent  to  the  lead  stop  in  a  partial  tour,  for  most 
values  of  k,  not  all  stops  can  feasibly  be  included  in  Sg.  for  example, 
any  customer  j,  either  an  origin  or  a  destination,  may  have  an  earliest 
service  time  Cj  and  a  latest  service  time  t  ,  defined.  If 


j 


1  71  k 
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cannot  feasibly  be  contained  in  Sg.  I.ct. 

=  Ol*j£V 

k 


(1.1a) 


The  number  of  fg  values  is  primarily  determined  by  the  number  of 
combinations  in  Sg,  which  is  determined  by  the  number  of  customers  to 


be  taken  over  the  k-1  remaining  positions.  From  expressions  (4.12)  and 
(4.13)  it  must  be  the  case  that 


j ell  =>  j  t  S,  (4.14) 

7Tk  k 

The  reduction  from  (4.14)  in  the  number  of  customers  that  can  be 
included  in  produces  a  marked  decrease  in  storage  requirements. 
Specific  examples  are  presented  below.  However,  even  further  savings 
in  storage  and  directly  related  computational  effort  can  be  achieved. 

Required  elements  in  Sg.  While  some  customers  can  not  be  included 
in  Sg,  it  is  often  true  that  others  must  be  if  the  final  tour  is  to  be 
feasible,  regardless  of  partial  tour  feasibility.  Consider  any  customer 
),  ei tlicr  an  origin  or  a  destination,  where 

ej  >  Ttk  .  (4 .15) 


If  customer  j  is  not  included  in  Sj.,  there  can  lie  no  feasible  completion 
of  the  total  tour,  since  sucli  a  completion  would  require  j  to  be 
sequenced  before  tt^,  in  violation  of  (4.15).  It  should  be  noted  that 
the  partial  tour  by  itself  might  well  be  feasible.  bet 

V.  =  fj  |e.  >  tt,  }  (4.15.1) 

k  J 

and  define 
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(4.1b) 


Then  !:r  represents  all  customers  for  which  a  decision  can  be  made  as  to 
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whether  or  not  they  are  to  be  included  in  a  given  partial  tour.  Using 
the  symbol  {• |  to  represent  the  set  cardinality,  the  upper  limit  on  the 
number  of  t'p  storage  locations  for  a  given  value  of  k  is  given  by 


ij,kl  ■ 


k  “  1  -  I  V  iT  k  I 


(4. 1~) 


The  precedence  requirements  make  the  actual  requirements  significantly 


Sample  storage  requirements.  Consider  a  travelling  salesman 
problem,  a  constrained  and  an  unconstrained  dial-a-ridc  service  problem 
of  site  N=13.  Suppose  k=7,  which  represents  the  halfway  point  in  the 
recursions,  further  suppose  that  for  the  constrained  problem  the  set 
J rr 7  consists  of  two  origins  and  three  destinations,  the  set  U-7 
consists  of  one  origin  and  one  destination  in  the  form  of  one  pair. 

Thus,  F  consists  of  three  origins  and  four  destinations.  Table  4 
shows  the  storage  requirements  to  compute  all  feasible  f?  values  for  all 
three  problems.  The  45  actual  locations  required  for  the  constrained 
problem  is  only  60°6  of  the  limit  of  75  computed  by  expression  (4.17). 


TABU-.  4 

hXAMl’I.h  fk  STOUAC.i:  UF.cjlIIUFMFNTS  l-OR  I1IR1.P  RFl.ATFP 
PROBI.PMS  WITH  Ni=  1 5  and  k  =  7 
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As  will  be  seen  below,  this  efficient  method  of  state  generation  and 
storage  allows  for  the  solution  of  much  larger  problems  than  previously 
solved. 

Existence  of  a  feasible  Solution 

Feasibility  question.  When  all  of  the  previously  suggested  problem 
constraints  are  present,  there  is  no  guarantee  that  a  feasible  solution 
even  exists.  The  more  heavily  constrained  a  given  problem  is,  the 
greater  the  probability  that  a  feasible  solution  does  not  exist.  The 
following  lemma  examines  the  feasibility  question  for  one  instance  of 
the  PUDP  where  capacity  constraints  arc  not  binding,  or  where  the 
capacity  arguments  suggested  earlier  in  this  chapter  arc  applicable. 

Suppose  that  for  each  customer  with  origin  i  and  corresponding 
destination  j,  all  relevant  constraints  can  be  expressed  in  terms  of 
the  stop  numbers  on  the  vehicle's  route.  Further  suppose  that  each 
customer  has  a  specified  stop  number  T.  representing  the  NET  time  ( >' j  | 
for  delivery  to  the  destination.  Hie  latest  time  the  corresponding 
origin  could  feasibly  be  visited  is  lb-1  since  at  least  one  time  period 
is  required  to  travel  between  any  two  points.  Consider  a  solution  P 
generated  by  taking  each  required  stop  in  a  nondecreasing  ordering  of 
the  not  later  than  times.  For  each  customer,  there  are  two  required 
stops:  one  for  the  origin  and  one  for  the  destination.  In  such  an 
ordering,  a  customer's  destination  would  always  follow  its  corresponding 
origin  since  t.  =  ib  >  ib-1  =  P j .  Consequently,  1'  would  always  satisfy 
the  precedence  requirement  of  the  PUDP. 


4 


< 


The  construction  of  P  can  be  demonstrated  by  t lie  following  simple 
example.  Suppose  there  are  three  customers  with  NLT  times  for  their 
destinations  of  4,  5  and  8  for  (destination  1),  and  1)_,  respec¬ 
tively.  Then  the  NLT  times  for  the  corresponding  origins  arc  3,  4,  and 
7  for  Oj,  O2  and  O3.  Then  P  would  be  the  sequence: 

Depot-01(3)-02(4J-l)1(4)-D2(5)-05(7)-D3(8)-[)epot, 


where  the  numbers  in  parenthesis  are  the  NLT  values. 


Lemma  3.  Given  a  not  later  than  time  for  each  delivery,  construct 
P  as  shown  above.  Suppose  Q  (the  quality  of  service  paramentcr)  is 
defined  to  be  the  maximum  possible  number  of  stops  before  a  customer 
reaches  his  destination.  Thus,  for  a  customer  with  origin  i,  destination 
j  and  Tj ,  the  not  earlier  than  times  are  computed  to  be  Tj-Q  and  Tj-Q+1 
respectively  (discussed  earlier  on  p.  08).  Then,  there  exists  at  least 
one  feasible  solution  if  and  only  if  I’  is  feasible.  To  be  feasible  a 
point  must  be  sequenced  within  its  time  window:  T^-Q  through  T ’• - l 
inclusive  for  the  origin  associated  with  destination  j  and  T.-Q+l 
through  Tj  inclusive  for  j. 

Proof.  The  if  part  is  obvious,  lor  the  only  if  part,  suppose  I' 

is  not  feasible.  Let  P*  be  any  other  sequence.  It  will  be  shown  that 

P*  cannot  be  feasible.  I.et  h  represent  the  number  of  stops  until  the 

first  infeasible  assignment  in  P  is  encountered.  let  j  represent  that 


infeasible  assignment,  lor  example,  consider  the  following  I': 
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Depot -0 -02'D1-0^-Do-D3'Depot . 

If  the  third  stop,  which  is  D^ ,  is  infeasible,  then  h=3  and  j=U^.  These 
are  two  cases  to  be  considered,  a  point  sequenced  cither  too  early  or 
too  late. 

Case  1;  ej  >  h.  In  this  case,  the  point  sequenced,  either  an 
origin  or  a  destination,  has  been  sequenced  too  early.  In  effect,  it 
will  be  shown  that  to  move  j  later  in  the  tour  requires  placing  another 
point  k  in  spot  h  which  has  c^  >  C|v 

Let  h1  be  the  first  stop  for  which  j  could  possibly  be  feasibly 
sequenced  in  P*.  Clearly,  h'  >  h  since  ej  >  h  in  P.  Hence,  in  P*. 
h’-l  points  must  be  sequenced  before  j.  Let 

H  =  (i(ei  <  ej}  .  (4. 18) 


Because 

Ci  =  -  Q+l,  i  =  1 ,2 . N  ,  (4. IP) 

it  must  be  true  that  <  ej  if  and  only  if  f.j  f.  j .  Therefore,  it 
follows  that 

I II  i  <  h  11 .  JO.) 

where  represent  the  cardinality  of  the  set.  Were  this  not  true,  j 
would  have  occupied  a  later  position  in  P.  lor  P*  to  be  feasible,  it 
must  be  true  that 


| H |  >_  h'  -1 


(  j  . :  i  j 


i.e.,  there  are  enough  points  to  fill  the  first  h'-l  positions. 

Equation  (4.21),  in  turn,  implies 

h  >  h'-l  (4 . 22 J 


or 


h  >  h'  (1.23) 

which  is  a  contradiction.  Ilcnce,  if  e j  >  h,  there  is  no  possible 
feasible  sequence. 

Case  2;  kj  <  h.  lor  the  second  case,  the  point  has  been  sequenced 
too  late.  To  he  feasible,  it  must  be  sequenced  earlier.  However,  to 
place  j  earlier  in  the  tour,  another  point  k  which  has  ip  <_  l;  put  into 
position  h. 

Let  h'  be  the  last  possible  stop  for  which  j  could  feasibly  be 
sequenced  in  !’*,  i.e.,  h'  =  ly  Using  an  argument  which  is  the  reverse 
of  Case  1,  h'  <  h  and  h '  -  1  stops  must  be  sequenced  before  stop  h 

•ft 

(customer  j)  in  1’  .  Let 

II*  =  {i  |  <  £j  }  (4.21) 

II*  is  the  set  of  customers  that  must  be  serviced  not  later  than  Pj. 

Since  I’  is  ordered  in  nondecreasing  5j,  it  must  be  true  that 

l»l*|  >  h  ,  (4.25) 


4 


> t  a  -1  —  - 


which  in  turn  implies 


111*  I  >  li'  >  h*  -  1 

Consequent  ly ,  I’*  cannot  he  feasible  since  there  are  more  custor 
must  he  scheduled  than  there  are  available  positions.  In  othci 
if  li.  j  <  h,  there  is  no  possible  feasible  sequence. 

Therefore,  one  concludes  that  if  I’  is  not  feasible,  no  otl 
permutation  can  be  feasible,  and  the  Lemma  is  proved. 

feasibility  in  general.  If  the  conditions  of  Lemma  o  do  i 
checking  for  feasibility  is  much  more  complicated.  lor  exampl* 
varies  per  customer  requirement  it  may  be  the  case  that  !’,  as 
developed  using  the  procedure  of  Lemma  5,  is  infeasible,  but  t 
slight  modification  producer  a  feasible  tour.  for  example,  co 
the  data  in  Table  5.  An  ordering  P  would  lie 

i-o1-o2-i)1-o^-iJ2-n3- 1 

which  is  infeasible  since  0^ ,  the  first  customer’s  origin,  is 
first,  thus,  violating  ci  =  2.  But 

i-o, -o.  -i>.-o_-i)-,-n,-i 

is  feasible.  Since  the  dynamic  programming  algorithm  only  gei 
feasible  state  vectors,  infeas ibi 1 i ty  is  recognized  when  1 f^ 
some  value  of  k.  However,  for  the  tour  construct  ion  heuristi 
discussed  in  the  next  chapter,  a  simple  technique  for  insurin 
feasibility  is  most  important. 


TABU.  5 

LXAMI'LL  DATA  I  OK  WHICH  T  HOI  S  NOT  VIi:l.H  A 
l  l.ASIBU:  SOLOMON 


NL'I 

CUSTOMLK  ULSTINATI ON/OKI  GIN  g 


1  4 

2  3 

3  (> 


Coinput  at  ion.i  1  l.xperi erpce 

Algorithm  descrijit  ion  .  I  ho  a  1  gor  i  t  hr;  used  to  solve  the 
vehicle  dial-a-ride  service  problem  is  an  extension  of  that 
only  precedence  constraints  were  binding.  'Hie  extension  me i 
porates  checks  to  determine  which  customer  roqui remen ts  sirs 
and  which  must  he  considered  for  each  value  of  k  and  choice 
lead  position  as  explained  above. 

Types  of  problems  run.  Any  number  of  problems  could  h; 
solved.  However,  because  the  literature  is  void  of  sample  ] 
which  heuristics  could  he  tested,  emphasis  was  placed  on  no. 
problems.  Medium  sited  is  defined  to  be  13  customers,  whirl 
31  specific  stops.  Using  the  previously  mentioned  approx ir; 
minutes  between  stops,  a  convenient  eight  hour  workday  is  o! 

Cons t ra ini  de f i n i t i on .  Both  a  quality  of  service  para: 
a  maximum  waiting  paramenter,  M,  are  used  to  dvtoinnie  trie 
varying  the  tightness  of  the  constraints,  Ibis  impact  been 
significant  in  studying  the  heuristic  solution  procedures  a 
seen  in  the  next  chapter.  Table  o  presents  the  tine  window 
for  one  problem  with  N  =  3 1 ,  Q  =  11,  and  M  -  <  . 

Results.  fable  ~  presents  typical  results.  1  he  s! or, 
represents  the  number  of  feasible  state  vectors  tlia;  .,re  g«. 
Actual  storage  required  for  X  =  31  is  twice  the  figure  pro- 
additional  storage  is  for  the  representation  vector.  lor  ‘ 
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TABLE  7 

SAMPLE  COMPUTATIONAL  REQUIREMENT  FOR  EXACT 
SOLUTION  BY  DYNAMIC  PROGRAMMING 


N  Q  M  STORAGE  REQUIRED  RUN  TIME 

(SEC) 


31 

31 

31 

31 

31 

91 


5 

5 

7 

11 

11 

7 


3 

5 

4 

6 

11 


174 
227 
4  S3 
784 
867 
1,328 


1.280 
1.433 
1.706 
3.057 
3.623 
1  .631 
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total  storage  required  is  four  times  the  figure  presented.  The 
additional  storage  is  for  the  three  representation  vectors,  one  for 
stops  1-31,  another  for  stops  32-62  and  the  last  for  stops  63-91.  The 
run  time  is  the  average  for  10  continuous  trials. 

Analysis  of  results.  No  certain  explanation  can  be  offered  as  to 
why  the  91  stop  problem  was  solved  more  quickly  than  was  the  similar  31 
stop  problem.  Perhaps  the  answer  relates  to  the  number  of  other  jobs 
in  the  time  sharing  system  at  the  two  different  times  these  problems 
were  solved.  Notwithstanding,  91  stop  problems  were  solved  in  less  than 
two  seconds,  requiring  allocation  of  less  than  150k  bytes  of  total 
computer  storage  for  compilation,  execution  and  data.  One  easily  con¬ 
cludes  that  any  practical  sized  single  vehicle  PUDP  for  which  the 
simplifying  assumptions  -  time  windows  and  reasonable  quality  of  service 
parameter  —  apply  can  be  optimally  solved  very  efficiently.  The  single 
vehicle  dial-a-ride  service  problem  appears  to  be  an  excellent  practical 
problem  for  which  dynamic  programming  could  be  used  to  obtain  the  exact 
routing  provided  that  computer  facilities  were  available  to  the  agency 
providing  the  service. 

IV.  USES  OF  OPTIMAL  SOLUTIONS 

Despite  these  impressive  results,  dynamic  programming  is  not 


practical  for  solving  the  multiple  vehicle  PUDP.  The  problem  is  that 
the  number  of  states  generated  cannot  be  limited  to  the  extent  they  can 
with  the  single  vehicle  case.  One  can  no  longer  require  customers  to  be 
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on  one  specific  vehicle.  They  could  just  as  easily  be  assigned  to 

another  vehicle  in  the  optimal  solution.  Consequently,  the  cardinality 

of  F  defined  by  expression  (4.16),  increases  significantly,  which 
k 

translates  into  a  tremendous  increase  in  the  number  of  f^  storage 
locations  for  a  given  value  of  k.  Results  pertaining  to  the  multiple 
vehicle  PUDP  are  contained  in  Chapter  VI. 

Since  the  PUDP  has  not  been  studied  in  the  literature,  with  the 
exception  of  the  few  small  problems  Psarftis  presents,  there  arc  no 
standard  test  problems  in  the  literature  against  which  to  compare 
various  heuristics.  The  dial-a-ride  service  problem,  optimally  solved 
by  the  dynamic  programming  algorithm,  provides  the  means  to  accurately 
and  precisely  test  the  numerous  heuristics  to  be  discussed  in  the  next 
chapter.  Having  the  exact  solution  is  extremely  valuable,  since  the 
results  indicate  that  many  of  the  most  popular  heuristics  for  the  TSP 
and  the  VRP  perform  very  poorly. 

V.  I. IMITATIONS  ON  DP  SOLUTIONS 


The  results  detailed  above  were  attainable  because  the  constraints 
could  somehow  be  expressed  as  a  function  of  stop  numbers.  This  allowed 
for  only  combinations  of  customers  in  a  partial  tour  to  be  considered. 
However,  if  the  exact  permutation  of  the  members  of  a  partial  tour  must 
be  known  to  assure  that  constraints  are  not  violated,  then  dynamic 
programming  will  not  work.  For  example,  if  quantities  to  be  moved  are 
not  unity,  it  is  easy  to  create  an  instance  where  one  permutation  of 
stops  is  feasible,  while  another  permutation  of  the  same  stops  is 


infeasible  due  to  vehicle  capacity  being  exceeded.  The  same  is  true  of 
time  window  constraints  if  stop  numbers  cannot  be  assumed  to  be  good 
enough  approximations.  Similar  instances  of  infeasibility  can  be  easily 
created  for  the  quality  of  service  requirement  and  the  operational  con¬ 
straints.  Consequently,  obtaining  optimal  solutions  appear  feasible  for 
only  a  rather  restricted  instrance  of  the  PliDP.  For  the  other  instances 
heuristic  solutions  must  suffice,  since  no  exact  solution  technique  has 
yet  been  found  that  is  effective. 


CHAPTER  V 


SINGLE  VEHICLE  HEURISTICS 

Not  all  instances  of  the  pickup  and  delivery  problem  (PUDP)  appear 
to  be  able  to  be  solved  optimally.  For  those  instances  which  cannot  be, 
hueristic  solutions  must  suffice.  This  necessity  for  reliance  on 
heuristics  was  earlier  noted  for  many  instances  of  the  travelling  sales¬ 
man  and  vehicle  routing  problems  as  well.  In  this  chapter,  therefore, 
heuristic  solutions  of  the  single  vehicle  PUDP  are  investigated. 

I.  OVERVIEW  OF  HEURISTIC  DEVELOPMENT 
Type  of  Instances  Studied 

Some  instances  of  the  PUDP  can  be  efficiently  solved  optimally  as 
shown  in  the  last  chapter.  These  instances  require  that  all  time 
related  constraints  be  expressed  in  terms  of  stop  numbers,  where  the 
stop  number  refers  to  the  number  of  points  actually  visited  in  any  par¬ 
tial  tour.  Capacity  must  not  be  a  factor,  or  the  other  constraints  must 
be  such  that  customers  do  not  remain  on  the  vehicle  long  enough  for  the 
capacity  of  the  vehicle  to  be  exceeded.  When  these  conditions  hold, 
dynamic  programming  can  be  used  to  obtain  the  exact  solution.  Very 
large  single  vehicle  dial-a-ride  service  problems  can  be  solved  when  any 
practical  quality  of  service  parameter,  Q,  is  used. 

To  say  that  one  heuristic  produces  solutions  that  are  on  average  so 
many  percentage  points  better  than  another  heuristic  is,  by  itself,  not 
too  meaningful.  If  both  heuristics  produce  results  that  arc  over  twice 
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optimal,  neither  heuristic  is  very  good.  However,  to  say  that  a 
heuristic  on  average  produces  solutions  within  so  many  percentage  points 
of  optimal  is  meaningful.  It  provides  a  precise  measure  of  how  well  a 
heuristic  performs.  For  this  reason,  the  problem  instances  studied  in 
this  chapter  are  the  same  as  those  studied  in  the  last  chapter.  This 
allows  for  precise  as  well  as  relative  evaluation  of  the  heuristics. 

Specific  Heuristics  Studied 

Each  of  the  heuristics  to  be  discussed  falls  into  one  of  the 
three  broad  classes:  tour  construction,  tour  improvement  or  composite. 
Because  the  Clarke-lVright  savings  technique  dominates  the  commercially 
available  vehicle  routing  packages,  it  was  the  logical  first  candidate 
for  study.  Other  tour  construction  algorithms  studied  include  the 
greedy  approach  and  an  insertion  technique.  The  3-optimal  approach  is 
the  primary  tour  improvement  heuristic  investigated,  while  each  of  the 
tour  construction  heuristics  is  combined  with  the  3-optimal  algorithm 
to  obtain  a  composite  solution. 

Each  heuristic  algorithm  is  first  discussed  individually,  except 
for  the  composite  algorithms  which  require  only  explanation  of  the 
individual  components.  The  algorithms  developed  are  not  straightforward 
extensions  of  those  used  on  either  the  travelling  salesman  problem  or 
the  vehicle  routing  problem.  As  was  the  case  with  many  of  the  exact 
solution  techniques,  the  precedence  relationship  is  the  complicating 
factor.  Computational  results  are  finally  presented  which  compare  the 
performances  of  each  of  the  various  heuristics.  The  results  clearly 
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show  that  tour  construction  heuristics  do  not,  in  general,  produce  good 
results.  Again,  the  precedence  relationship  is  the  primary  culprit. 

II.  CLARKE-WRIGHT  HEURISTIC 


Starting  Point  Procedures 

The  basic  Clarke-Wright  savings  model  was  developed  in  Chapter  II. 
For  the  pickup  and  delivery  problem,  the  depot  or  dispatch  point  is 
taken  to  be  the  origin,  and  the  first  two  points  are  taken  to  be  those 
with  the  largest  savings  value  as  computed  by 


s.  . 
ij 


=  c 


li 


+  c 


c.  . 
IJ 


(5.1) 


which  can  feasibly  be  linked  together.  When  time  windows  are  present, 
as  in  a  dial-a-ride  situation,  many  points  cannot  be  linked  due  to  the 
time  requirements  involved. 


Criteria  for  linking.  Two  points  can  be  linked  only  if  their  time 
windows  touch  or  overlap.  It  is  important  to  remember  that  all  con¬ 
straints  are  being  expressed  in  terms  of  a  vehicle's  sequence  number. 
Thus,  an  example  time  window  of  4  through  7  inclusive  for  a  given  point, 
representing  a  customer's  origin  or  destination,  would  require  the 
vehicle  to  visit  that  point  on  its  4th,  5th,  6th  or  7th  stop  after 
leaving  the  depot.  Let 

SR  =  {stop  numbers  in  k’s  time  window}  , 

k  =  2,3, ... ,N.  (5.2) 

Then  points  i  and  j  can  only  be  linked  if  the  set 
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Lij  =  siusj 


(5.3) 


contains  all  of  the  stop  numbers  between  the  lowest  and  the  highest 

values.  For  example,  suppose  S.  =  {2,3,4}  and  S.  =  {6,7,8}.  Then 

t  J 

=  {2, 3, 4,6, 7, 8}.  Because  stop  5  t  L. ^ ,  points  i  and  j  cannot  be 
linked.  If  i  were  sequenced  at  its  latest  possible  time  4,  j  would  have 
a  sequence  time  of  5  which  is  one  stop  earlier  than  its  earliest  time. 
Similarly,  sequencing  j  at  its  earliest  time  of  6  would  require  i  to  be 
sequenced  at  5,  one  stop  beyond  its  latest  time.  Since  the  Clarke- 
Wright  heuristic  does  not  allow  for  the  breaking  of  a  link  once  formed, 
the  linking  of  i  to  j  would  preclude  a  feasible  tour  from  ever  being 
constructed.  Further,  if 


S.  n  S.  =  0 
i  J 


(5.4) 


all  elements  of  S^  are  less  than  those  of  ,  then  i  must  be  linked 


to  j  in  a  forward  manner,  i.e.. 


1  -  i  -  j  -  1 


in  the  initial  partial  tour.  The  same  is  true  if  the  roles  of  i  and  j 
are  reversed.  Suppose  S.  =  {2,3,4}  and  S.  =  {5,6,7}.  Then  D  Sj  =  0. 
Clearly  any  tour  having  i  and  j  feasibly  linked  would  require  i  in  the 
4th  position  and  j  in  the  5th.  Any  other  linking  would  violate  one  of 
the  time  windows. 

For  other  cases,  one  has  a  choice.  Suppose  =  {2, 3, 4, 5, 6, 7}  and 
S.  =  {4, 5, 6, 7, 8, 9}.  If  i  precedes  j,  then  point  i  can  feasibly  occupy 


position  3,  4,  5,  6,  7  or  8  in  a  complete  tour.  However,  letting  j 
precede  i  limits  i  to  occupy  only  positions  5,  6,  or  7.  Hence,  i 
preceding  j  allows  greater  flexibility  in  subsequent  construction  of 
the  solution.  Therefore,  if  £j  <_  £j ,  then  the  algorithm  has  i  preceding 
j  in  the  initial  partial  tour. 

Initial  feasibility.  Although  the  above  conditions  are  necessary 
for  feasibility,  they  are  not  sufficient  to  insure  that  a  feasible  tour 
can  be  constructed.  If  linking  i  and  j  precludes  the  subsequent  con¬ 
struction  of  a  feasible  tour,  then  i  and  j  cannot  be  linked.  Consider 
a  simple  example,  with  partial  data  as  given  in  Table  8.  Linking  i  and 
j  appears  feasible  with  i  in  the  stop  5  position  and  j  in  the  stop  6 
position.  However,  with  i  and  j  in  these  positions,  point  b  cannot  be 
feasibly  sequenced.  Yet  the  sequenced 

1-x-x-x-i-b-j-x-. . .-1  (5.6) 

is  feasible  with  respect  to  i,  j  and  b.  The  x's  represent  arbitrary 
other  points. 

Therefore,  in  selecting  the  first  two  points  to  be  linked,  it  is 
necessary  to  insure  that  a  feasible  final  tour  can  be  constructed.  When 
the  constraints  are  relatively  tight,  many  infeasible  linkings  are 
likely  to  have  higher  savings  values  than  those  of  the  points  that  can 
be  feasibly  linked.  To  ignore  feasibility  during  the  construction 
process  would  almost  always  result  in  an  infeasible  solution.  Lvcn  if 
one  had  a  90%  chance  of  selecting  a  feasible  linking  at  each  stage,  for 


TABLE  8 

EXAMPLE  DATA  TO  DEMONSTRATE  THE  PROBLEMS 
ASSOCIATED  WITH  LINKING  FEASIBILITY 


POINT 


NET 


NLT 


a  problem  with  15  customers  (30  origin/destination  points)  the 
probability  of  not  selecting  an  infeasible  linking  (without  insuring 
feasibility)  is  less  than  5%.  Further,  insuring  a  feasible  linking  is 
not  sufficient  to  guarantee  a  feasible  final  tour  as  discussed  above. 
Consequently,  the  probability  of  a  feasible  final  tour  (without  insuring 
feasibility)  is  near  zero. 

The  checking  process  is  complicated  and  must  be  accomplished  at 
each  step  before  a  new  point  is  added  to  the  solution.  Assuring  a 
feasible  final  tour  is  necessary  if  the  heuristic  is  to  have  practical 
value.  Infeasible  solutions  are  considered  worthless.  The  detailed 
explanation  of  this  procedure  will  be  deferred  until  later  in  this 
section. 

Adding  Additional  Points 

Recall  that  the  Clarke-Wright  algorithm  adds  new  points,  one  at  a 
time,  to  an  existing  partial  tour.  The  new  point  added  is  always 
connected  to  the  origin  and  one  end  of  the  partial  tour.  Once  two 
points  are  linked,  they  remain  linked.  For  ease  of  explanation,  at  any 
step  in  the  process,  those  points  which  arc  already  in  the  partial  tour 
will  be  referred  to  as  in  the  assigned  set,  A.  The  order  in  the 
assigned  set  will  be  that  of  the  partial  tour.  The  point  being  con¬ 
sidered  for  addition  to  the  partial  tour  is  the  candidate.  The  candi¬ 
date  can  be  inserted  either  in  front  of  or  behind  the  assigned  set, 


provided  such  an  insertion  can  lead  to  a  feasible  final  tour.  Therefore, 
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at  most  two  feasible  candidates  are  generated  at  each  step  of  the 
iteration,  a  front  end  candidate  and  a  rear  end  candidate.  It  is 
possible  that  one,  but  not  both,  of  these  candidates  will  not  exist, 
especially  near  the  end  of  the  construction.  The  partial  tour  is 
expanded  by  adding  to  the  tour  the  candidate  which  has  the  largest 
savings.  In  this  sense,  the  heuristic  is  applied  exactly  in  the  same 
manner  for  the  PUDP  as  it  is  for  the  TSP  or  the  Y'RP. 

Front  end  candidate.  There  are  three  possible  ways  that  new  point 
i  could  be  the  front  end  candidate,  disregarding  feasibility: 

•  i  is  an  unassigned  origin,  whose  corresponding  destination 
is  already  sequenced. 

•  i  is  an  unassigned  origin  as  is  its  corresponding  destination, 
but  the  assigned  set  is  small  enough  that  the  destination  can 
be  sequenced  later. 

•  i  is  an  unassigned  destination  whose  origin  has  not  yet  been 
sequenced . 

Checking  feasibility  is  expensive  in  terms  of  computational  effort. 
However,  simple  checks  can  sometimes  eliminate  a  possible  candidate 
which  passes  the  above  criteria  and  has  a  greater  savings  value  than  the 
present  incumbant.  These  checks  arc  based  on  the  earliest  time  that 
point  i  can  be  served  (c ^ )  and  the  latest  time  that  it  can  be  served 
(&i).  Let  i  be  the  candidate  and  j  the  first  element  of  the  ordered  set 
A,  which  implies  that  i  would  directly  precede  j.  If  cither 

*i  <  cj  -  1  (5.71 


< 


or 


5)4 


oA  1  (5.8) 

hold,  then  i  is  not  feasible  in  the  proposed  partial  tour,  and  can  be 
eliminated  from  further  consideration.  Otherwise,  the  more  complicated 
feasibility  check  is  used. 

Expressions  (5.7)  and  (5.8)  can  be  illustrated  with  two  simple 
examples.  First,  suppose  =  3  and  ej  =  5  so  that  (5.7)  holds.  If  i 
is  placed  in  position  3,  j  must  be  in  position  4,  which  violates  ej  =  5. 
If  j  holds  position  5,  i  must  be  in  position  4,  which  violates  £j_  =  3. 
Second,  suppose  e^  =  £j  =  4.  Then  i  in  position  4  requires  j  in  5, 
violating  =  4,  while  j  in  position  4  requires  i  in  3,  violating 
ei  =  4. 

Rear  end  candidate.  There  are  three  possible  ways  that  point  i 
could  be  a  rear  end  candidate: 

•  i  is  an  unassigned  destination  whose  corresponding  origin  is 
already  sequenced. 

•  i  is  an  unassigned  destination  whose  corresponding  origin 
has  not  yet  been  sequenced,  but  the  assigned  set  is  small 
enough  that  the  origin  could  be  sequenced  before  the  first  stop. 

•  i  is  an  unassigned  origin  whose  destination  has  not  yet  been 
sequenced . 

Letting  i  again  be  the  candidate  and  letting  j  be  the  last  ordered 
element  of  A,  which  imples  a  j-i  linking,  i  can  be  ignored  if  either 


« 


l3 


+  1  <  ej 
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(5.9) 


or 


ej  >  h 


(5.10) 


hold. 

Again,  two  examples  will  illustrate  why  these  expressions  are  true. 
First,  suppose  Hj  =  6  and  ej  =  8  so  that  (5.9)  holds.  If  j  occupies  its 
last  possible  position,  position  6,  point  i  must  occupy  position  7.  If 
i  is  in  position  8,  j  must  be  in  position  7.  In  both  cases,  one  of  the 
points  violates  a  time  window.  For  expression  (5.10),  suppose 
ej  =  =  9.  Since  j  precedes  i,  there  is  no  way  both  windows  can  be 

satisfied.  The  situation  becomes  worse  as  Cj  increases,  i.e.,  j  cannot 
be  served  until  later  times. 

Candidate  selection.  Suppose  that  two  feasible  candidates  have 
been  identified.  The  one  with  the  greatest  savings  is  added  to  the 
partial  tour.  If  both  candidates  have  the  same  savings  and  the  front 
end  candidate  is  an  origin,  it  is  sequenced.  Otherwise,  the  rear  end 
candidate  is  added  to  the  assigned  set.  It  is  not  necessarily  possible 
to  add  both  the  front  end  and  the  rear  end  candidate  to  a  partial  tour 
in  the  event  of  a  tic  in  savings.  Final  tour  feasibility  can  be 
destroyed,  as  was  discovered  during  initial  implementation  of  the 
algorithm.  In  some  cases  the  front  end  candidate  and  the  rear  end  candi¬ 
date  are  the  same  point. 


if 


4 


Constructing  a  Feasible  Tour 

When  only  the  precedence  constraints  are  binding,  it  is  a 
relatively  simple  matter  to  construct  a  feasible  tour  using  the  simple 
checks.  However,  when  time  windows  are  present,  guaranteeing  that  a 
partial  tour  can  be  extended  into  a  legitimate  final  solution  is  much 
more  complicated. 

Ordering  of  points.  Lemma  3  of  Chapter  IV,  establishes  necessary 
and  sufficient  conditions  for  the  existence  of  a  feasible  solution  for 
one  type  of  the  PUDP.  Although  the  lemma  is  not  valid  for  all  PUDP 
instances,  the  nondecreasing  ordering  of  the  NLT  times  is  a  useful  and 
reasonable  heuristic  for  testing  feasibility.  As  before,  let  P  denote 
the  sequence  of  points  ordered  by  nondecreasing  NLT  times.  Therefore, 
in  P  each  customer  point,  origin  or  destination,  has  a  position  number, 
which  would  correspond  to  the  stop  number  if  P  were  a  proposed  tour. 

One  assumes  P  is  feasible  so  that  at  least  one  feasible  tour  exists.  If 
P  is  infeasible,  the  algorithm  terminates. 

Proposed  tour  Y.  To  test  feasibility,  a  complete  tour  will  be 
constructed  using  the  already  constructed  partial  tour,  the  candidate 
point,  and  the  ordering  P.  The  complete  tour  to  be  so  constructed  is 
termed  Y.  The  ordering  P  becomes  a  starting  point  from  which  the 
proposed  tour  Y  is  created.  Y  is  initially  a  null  vector.  If  the  final 
Y  is  feasible,  the  candidate  can  be  feasibly  linked  to  the  partial  tour. 
If  y  is  not  feasible,  the  proposed  condidatc  is  discarded;  there  is  no 


assurance  that  a  feasible  tour  can  be  constructed. 
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Concept  of  Y.  Y  is  constructed  one  point  at  a  time  in  reverse 
order.  Thus,  the  last  position  is  the  first  filled,  the  next  to  last 
position  is  the  second  filled,  etc.  The  ordering  P  is  used  whenever 
possible.  P  is  also  examined  in  reverse  order.  Thus,  the  last  element 
of  P  not  already  in  Y  is  first  considered  to  fill  the  next  position.  It 
may  be  that  other  points  or  the  entire  partial  tour  will  be  added  at 
this  point.  The  decision  is  based  on  assigning  each  point  as  late  as 
feasibly  possible,  which  in  turn  increases  the  overall  likelihood  of  a 
feasible  complete  tour. 

Latest  position  in  A.  Let  L  represent  the  position  number  of  the 
point  in  A  which  has  the  highest  position  number  in  P.  It  is  not 
necessary  that  the  point  in  question  be  the  last  ordered  element  in  A. 
More  often  than  not,  it  will  not  be.  Suppose  there  are  10  points  in  P 
and  L  =  7.  This  means  that  those  points  in  positions  8,  9  and  10  of  P 
are  not  yet  in  the  Clarke-W'right  generated  set  A.  Therefore,  all  points 
in  P  with  position  numbers  after  L  can  be  assigned  identical  positions 
in  Y.  The  set  A  may  or  may  not  be  the  next  points  added  to  Y.  If  A 
contains  origins  whose  destinations  are  not  in  A  nor  yet  in  Y,  these 
must  be  inserted  before  A  is,  to  preclude  violation  of  the  precedence 
constraints.  An  example  test  will  be  used  to  illustrate  the  procedure. 
For  ease  of  explanation,  the  customer  points  will  be  assigned  letters. 
Table  9  gives  hypothetical  time  windows  for  a  problem  with  N  =  11. 

Figure  5  depicts  the  associated  sequence  P,  a  given  assigned  set  A, 
and  the  initial  assignment  of  stops  to  based  on  L.  The  fact  that 


TABLE  9 

HYPOTHETICAL  TIME  WINDOW  USED  TO  DEMONSTRATE 
TOUR  CONSTRUCTION  FEASIBILITY 


ORIGINS 

DESTINATIONS 

POINT  # 

NET 

NLT 

POINT  ff 

NET 

NET 

a 

1 

4 

f 

2 

5 

b 

1 

5 

g 

i 

(> 

p  = 

a  1 
I 

b 

f  1 

g 

c 

h 

d 

_ _ 

e 

i 

j 

(a)  P: 

123456789  10 

the  nondecreasing  ordering  by  NLT  times 

A  = 

r- 

Lh_j 

b 

0>) 

Given  assigned  set  A  including  candidate 

Y  = 

d 

c 

i 

j 

123456789  10 

(c)  Initial  construction  of  V 


Figure  5.  Initial  construction  of  the  proposed 
tour  Y  using  P. 
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L=6  is  based  on  the  observation  that  point  h  occupies  position  6  in 
P,  whereas  the  other  points,  d,  e,  i  and  j,  which  correspond  to  stops 
7,  8,  9  and  10  in  P,  occupy  identical  positions  in  Y. 

Destination  check.  If  the  assigned  set,  A,  contains  origins  whose 
corresponding  destinations  are  not  in  A,  these  destinations  must  follow 
A.  Further,  if  these  destinations  have  position  numbers  in  P  less  than 
L,  they  are  not  included  in  the  already  sequenced  points  in  Y.  Such 
stops  are  inserted,  in  nondecreasing  order  of  NLT  times,  before  the 
points  in  A  are.  In  the  example,  destination  g,  corresponding  to  origin 
b,  which  is  included  in  A,  has  a  position  in  P  of  4.  Consequently,  g 
is  not  yet  included  in  Y,  but  must  follow  all  elements  of  A. 

Figure  6(a)  shows  stop  g  sequenced. 

Flushing  out  Y.  The  remaining  process  places  A  directly  preceding 
the  first  sequenced  point  in  Y,  and  fills  out  by  taking  the  as  yet 
unsequenced  stops  in  the  same  order  as  in  P.  Figure  6(b)  depicts  the 
final  proposed  tour,  which  is  then  checked  for  feasibility.  The  example 
data  results  in  a  proposed  tour  Y  which  is  feasible.  Therefore,  the 
candidate,  either  c  or  b,  can  be  added  to  the  partial  tour,  and  will  be 
it  it  subsequently  represents  the  greatest  savings. 

Clarkc-Wright  Solution 


This  process  is  continued  until  the  final  tour  is  constructed. 
Final  tour  feasibility  is  assured  because  a  candidate  is  added  to  any 
partial  tour  only  if  a  feasible  final  tour  has  actually  been  verified. 
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The  final  Clarke-Wright  solution  was  also  used  as  an  input  to  the 
3-optimal  algorithm,  thus  forming  a  composite  heuristic.  Both  the 
3-optimal  algorithm  and  results  of  the  computational  experiments  are 
presented  later  in  this  chapter. 

III.  ROUT.1;  INSERTION  HEURISTIC 
Description  of  Concept 

The  route  insertion  heuri  x  combines  the  philosophy  of  the 
insertion  heuristics  with  the  savings  concept  of  the  Clarke-Wright 
heuristics.  One  of  the  expl.-’i:  (tions  offered  for  why  the  Clarke-Wr i ght 
heuristic  works  so  well  on  the  travelling  salesman  problem  (TSP)  and 
vehicle  routing  problem  (VRP)  is  that  the  more  distant  customers  are 
considered  early  for  inclusion  in  the  tour.  This  procedure  precludes 
expensive,  last  minute  adjustments  to  sequence  these  customers.  The 
route  insertion  procedure  will  attempt  to  incorporate  this  feature  of 
the  Clarke-Wright  model . 

The  insertion  procedures  for  the  TSP  and  the  VRP  call  for  the 
selection  of  the  next  single  point  to  be  added  to  a  partial  tour  based 
on  some  selection  criteria.  Because  of  the  precedence  relationship  of 
the  PUDP,  origin/destination  pairs,  rather  than  a  single  point,  arc 
selected  for  insertion  into  the  partial  tour.  In  this  manner,  the  route 
insertion  heuristic  specifically  addresses  the  precedence  requirement. 
Further,  when  an  origin  is  sequenced,  the  vehicle  must  subsequent  1 y 
travel  to  the  corresponding  destination.  Therefore,  considering  the 
locations  of  both  the  origin  and  the  destination,  relative  to  a  partial 
tour,  seems  both  logical  and  advantageous. 


I  .■+ 
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Selection  Criteria 

Pairs  are  considered  for  insertion  in  nonincreasing  order  of 
individual  service  cost.  Individual  service  means  a  pair  is  serviced 
directly  from  the  depot  and  the  vehicle  returns  to  the  depot  upon  com¬ 
pletion.  For  origin  i  and  corresponding  destination  j,  the  individual 
service  cost  would  be 


li 


+  c . 


ij 


Cjl 


(5.11) 


The  pair  with  the  highest  individual  service  cost  starts  the  process. 

The  highest  cost  pair  is  taken  first  so  that  the  most  difficult  or 
costly  customers  to  serve  are  considered  first,  just  as  they  were  in  the 
basic  Clarke-Wright  model.  Thus,  if  pair  (i,j)  has  the  greatest  cost, 
the  initial  partial  tour  would  be 

1  -  i  -  j  -  1  (5.12) 


where,  as  before,  the  depot  is  designated  as  the  first  point.  The  pair 
with  the  greatest  remaining  cost  is  the  next  pair  to  be  inserted  into 
the  partial  tour,  and  so  on. 


Insertion  Criteria 

The  next  pair  is  inserted  into  a  partial  tour  so  as  to  minimize 
the  total  increase  in  cost.  Two  cost  formulas  apply  depending  on 
whether  the  new  destination  is  inserted  directly  after  its  origin,  or 
at  a  point  later  in  the  partial  tour.  Figure  7  (a)  depicts  a  partial 
tour.  Pair  (i,j)  is  to  be  inserted  into  this  partial  tour.  Figure  7  (b) 
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Figure 


1-a-b-c-d-l 


(a)  Partial  tour 


i  -  j 
/  \ 

1-a  b-c-d-1 


(b)  One  insertion  of  pair  (i,j) 


1-a 


i  j 

\  '  \ 
b  -  c  d  - 


1 


(c)  Another  insertion  of  pair  (i,j) 


7.  Example  insertion  patterns  for  the  pair 
insertion  heuristic. 


< 
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depicts  one  possible  arrangement  in  which  j  directly  follows  i.  The 
increase  in  cost  is 


c  . 
ai 


+  c.  . 
ij 


+  c 


jb 


'ab 


(5.13) 


Figure  7  (c)  depicts  the  more  common  occurrence  where  there  are 
intervening  stops  between  i  and  j.  For  this  configuration  the  increase 
in  cost  is 


c  .  +  c.,  -  c  ,  +  c  .  +  c.  ,  -  c  ,  .  (5.14 

ai  lb  ab  cj  jd  cd 

In  either  case,  the  insertion  of  i  and  j  must  be  such  that  a  feasible 
final  tour  can  be  constructed.  Insuring  that  a  feasible  final  tour  can 
be  constructed  from  a  partial  tour  is  absolutely  necessary,  as  noted 
above  in  the  discussion  of  the  Clarkc-Kr ight  heuristic.  A  somewhat 
similar  procedure  is  used  to  verify  the  existence  of  such  a  feasible 
tour . 


Feasibility  of  Insertion 

Rough  checks.  If  point  y  is  inserted  directly  after  point  x,  point 
x  will  precede  y  in  all  subsequent  tours,  although  x  will  not  directly 
precede  y  if  other  stops  are  inserted  between  them.  Let  (i,j)  be  the 
pair  to  be  inserted  and  let  a-b  be  the  two  points  between  which  i  is 
to  be  considered  for  insertion.  Then  i  cannot  be  feasibly  inserted 
between  a  and  b  if  either 

c  >  £.  ,  (5.15) 

a  —  i 


I 


4 


e .  >  £, 
l  —  b 


(5.16) 


Botli  represent  configurations  that  cannot  possibly  be  feasible. 
Similarly,  let  c-d  be  the  two  points  between  which  j  is  to  be  considered 
for  insertion.  Then  j  cannot  be  feasibly  inserted  if  either  of  the 
following  hold: 


c  >  £ . 
c  -  J 


( 5 .  n 


e  .  >  1 . 
J  - 


(  5 .  1 S  ) 


More  complicated  check.  When  the  rough  check  indicates  that  a 
given  insertion  may  be  feasible,  a  somewhat  more  complicated  check  can 
be  used.  This  check  insures  relative  feasibility  among  the  element^ 
of  the  partial  tour  which  includes  the  pair  currently  being  inserted. 
Let  k  represent  the  difference  between  the  position  arbitrary  point  b 
occupies  in  the  partial  tour  and  the  position  point  a  occupies.  If 


c  +  k  >  £, 
a  b 


(5. 19) 


for  any  a  and  b  in  the  partial  tour  where  a  precedes  b  then  t lac  partial 
tour  cannot  be  feasibly  extended. 


i 
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Detailed  Feasibility  Check 

When  a  proposed  insertion  passes  all  of  the  above  checks  it  is 
still  possible  that  a  feasible  tour  cannot  be  constructed.  Therefore, 
a  complete  tour  is  generated  using  the  ordering  P  discussed  earlier.  If 
this  tour  is  feasible,  the  proposed  insertion  positions  replace,  or 
become,  the  incumbent  positions.  In  either  case,  the  next  set  of 
positions  are  then  checked  and  the  process  repeated  until  the  feasible 
positions  that  minimize  either  (5.13)  or  (5.14)  are  found. 

Difference  in  concept.  The  primary  difference  between  checking 
feasibility  with  route  insertion  and  with  the  Clarke-Wright  approach  is 
that  with  route  insertion  only  the  order  within  the  partial  tour  A  must 
be  preserved.  With  Clarke-Wright ,  the  entire  set  A  was  sequenced  as  a 
block . 


Construction  of  Y.  A  complete  tour  V  is  constructed  and  then 
checked  to  see  if  it  is  feasible.  V  is  constructed  one  point  at  a  time 
in  reverse  order  using  P  and  the  ordered  set  A.  Iilements  from  set  A  are 
sequenced  in  Y  when  they  are  the  last  clement  and  normally  appear  in  the 
reverse  sequence  of  P  or  when  the  last  element  of  A  must  be  sequenced 
to  preclude  violation  of  a  "not  earlier  than"  requirement. 

Modification  of  e j .  When  points  appear  in  a  partial  tour,  the 
earliest  time  one  stop  can  actually  be  visited  may  be  altered  by  a  pre¬ 
vious  point,  for  example,  consider  the  partial  tour  a-b-c-d,  and 
suppose  ea  is  greater  than  any  of  the  other  times.  Recall  that  all  time 


windows  are  defined  by  the  stop  numbers  of  the  vehicles  route.  Thus, 
direct  travel  between  two  points  requires  one  stop.  I'hen  the  earliest 
that  b,  c  and  d  can  be  visited  i.s  ea+l,  ea  +  2,  and  ea  +  5  respectively. 

Were  any  of  these  points  visited  earlier,  ca  would  be  violated,  for 
example,  if  a  tour  were  being  constructed  in  reverse  sequence,  stop  d 
must  be  sequenced  by  the  ea+3  position  (appear  in  the  tour  at  e..+3  or 

later)  or  ea  will  be  violated.  The  concept  of  adjusting  e;'s  is  used 

in  constructing  Y. 

fxamplc  Y.  Consider  the  example  P  given  in  figure  5(a),  on  p.  up, 
and  time  windows  defined  by  the  data  in  Table  S,  on  p .  91.  Suppose 
that  the  partial  tour,  made  up  of  two  pairs,  is 

1-b-d-g-i-l  .  (5.20) 

Table  10  gives  the  revised  not-earl i cr-than  times  for  each  of  the 
points  in  the  partial  tour.  Note  the  times  have  changed  for  points  g 
and  i.  figure  8  gives  the  step  by  step  construction  of  Y.  Steps  1 
insert  the  last  three  elements  of  P.  At  step  4,  d  is  not  added  to  Y 
because  g,  not  d,  is  on  the  end  of  the  not -yet-sequenced  points  in  A. 

Therefore,  h,  the  point  preceding  d  in  P,  is  inserted  at  this  step. 

At  step  5  and  (>,  g  and  d  are  inserted  because  of  their  adjusted  not- 
earl i er-than  times.  The  remaining  points  are  inserted  based  on  their 
position  in  P.  The  final  Y  is  seen  to  he  feasible.  A  forward  con¬ 
struction  using  adjusted  not  - 1 atcr-than  times  could  also  have  been  used. 


TABU:  10 

KfiVl Sl;D  NOT  TARUTK  THAN  TIMI.S  I  OK  Till.  POINT 
IN  A  TOR  CXAMPLi:  TROBLTM 


POINT  n 


NTT 
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IV.  GREEDY  HEURISTIC 

The  "greedy"  or  "nearest  neighbor"  heuristic  is  significantly 
simpler  to  implement  than  either  the  Clarke-Wright  or  the  route 
insertion  algorithms.  One  begins  at  the  depot  and  always  proceeds  to 
the  closest  point  which  can  lead  to  a  feasible  completion. 

Candidate  Set 

At  each  step,  a  candidate  set  is  identified.  Initially,  the 
candidate  set  is  made  up  of  all  origins  that  can  be  visited  during  the 
first  time  period,  i.e.,  all 

j  e  0  3  ej  <  2  .  (5.21) 

The  closest  origin  is  added  to  the  tour.  Its  destination,  if  feasible, 
and  all 

j  e  0  3  e-  <  3  (5.22) 

are  added  to  the  candidate  set  for  the  next  stop.  So  long  as  a  feasible 
complete  tour  can  be  attained,  this  process  continues  until  all  customers 
are  scheduled.  The  candidate  set  at  each  point  consists  of  unsequenced 
origins  and  of  unsequenced  destinations  whose  corresponding  origins  are 
sequenced,  both  of  whose  time  windows  overlap  the  next  stop  number. 
Arbitrary  selection  of  the  nearest  point  among  this  set  will  not 
necessarily  lead  to  a  feasible  tour. 
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Feasibility  Check 

The  nearest  neighbor  is  selected  from  the  candidate  set  and  a  trial 
complete  tour  Y  is  constructed.  If  Y  is  feasible,  the  point  selected 
is  added  to  the  tour.  If  Y  is  not  feasible,  the  point  is  removed  from 
the  candidate  set  and  the  process  repeated  until  a  feasible  point  is 
found.  As  before,  the  sequence  P  is  used,  this  time  in  a  forward 
manner.  The  k  stops  selected  by  the  greedy  procedure  become  the  first 
k  stops  in  Y.  All  remaining  stops  are  then  sequenced  using  the  same 
relative  order  as  these  stops  have  in  P. 

Random  Feasible  Tours 

The  candidate  set  and  feasibility  check  used  for  the  greedy 
heuristic  were  also  used  with  a  random  selection  of  the  elements  from 
the  candidate  set.  These  tours  were  used  as  a  random  input  to  the 
3-optimal  heuristic.  The  greedy  as  well  as  the  randomly  generated  tour 
results  will  be  discussed  after  the  tour  improvement  heuristics  are 
examined . 

V.  TOUR  IMPROVEMENT  HEURISTICS 


A-Optimal ity 

The  A-optimal  heuristic  of  Lin  and  Kerninghan,  discussed  in 
Chapter  II,  is  one  of  the  most  powerful  heuristics  available  for  use  on 
the  travelling  salesman  problem  and  the  vehicle  routing  problem.  One  of 
the  keys  to  the  computational  efficiency  of  the  algorithm  is  that 
profitable  reconnection  patterns  are  pursued  only  if  it  is  possible  to 


4 


113 


actually  exchange  the  A  arcs  and  obtain  a  feasible  solution.  Insuring 
feasibility  is  relatively  easy  for  the  ISP,  but  much  more  complicated 
for  the  PUDP. 

Consider  the  example  depicted  in  f  igure  '.Hal  .  There  is  no 
possible  way  to  remove  two  arcs  and  reconnect  them  in  a  new  pattern  such 
that  a  new  tour  satisfies  the  precedence  requirements.  Vet  Figure  9(b) 
shows  a  feasible  3-optima  1  reconnection.  The  A-optimal  approach,  which 
only  continues  if  the  arcs  already  removed  can  be  feasibly  rejoined, 
would  have  stopped  at  >  =  2  and  never  have  investigated  the  configura¬ 
tion  shown  in  Figure  9(b). 

One  should  note  that  for  a  directed  TSP  tour,  no  2-arc  exchanges 
are  possible.  A  directed  TSP  is  a  TSP  on  a  directed  network,  i.e., 
arcs  only  go  in  one  direction  between  nodes.  The  2-arc  reconnection 
pattern,  shown  earlier  in  Figure  9,  requires  one  of  the  aths  remaining 
after  the  arcs  are  removed  to  be  traversed  in  the  reverse  direction. 
Since  reverse  arcs  do  not  exist  m  the  directed  network,  no  2-arc 
exchanges  are  possible.  'I wo -arc  exchanges  are  quite  liielv  in  the  I'llhP. 
Traversing  paths  in  a  reverse  direction  is  not  prohibited  m  the  Fillip 
unless  doing  so  violates  a  precedence  constraint  or  one  ot  the  time 
window  constraints.  Traversing  paths  in  the  reverse  direction  is 
discussed  in  greater  detail  later  in  this  chapter  as  part  of  the 
discussion  of  an  r-optimal  algorithm.  Ihe  implication  of  this 
discussion  is  that  there  is  no  way  of  knowing  a  priori  wlutlur  or  not 
the  removal  of  more  arcs  will  eventually  produce  feasibility.  Ihe 
situation  depicted  in  Figure  9(a)  is  a  fairly  common  occurrence. 
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When  the  problem  constraints  are  relatively  tight,  many 
reconnections  will  appear  profitable  from  a  cost  savings  point  of  view. 
Therefore,  the  alternative  of  continuing  until  no  more  savings  are 
possible  and  then  testing  feasibility  would  be  both  inefficient  and 
ineffective.  No  other  alternative  was  found  that  allowed  for  modifica¬ 
tion  of  the  A-optimality  concept  so  that  it  could  be  applied  to  the 
PUDP.  This  is  yet  another  example  of  how  the  precedence  constraints  of 
the  PUDP  complicate  or  relegate  ineffective  solution  techniques  that  are 
very  powerful  when  applied  to  the  TSP  or  the  VRP. 

r-Optimality 

The  concept  underlying  the  r-optimal  heuristic  was  presented  in 
Chapter  II.  A  3-optimal  algorithm  was  developed  and  used  to  solve  the 
PUDP.  The  3-optimal  heuristic  was  used  on  randomly  generated  feasible 
tours  as  well  as  in  conjunction  with  solutions  generated  by  the  tour 
construction  heuristics. 

Order  of  complexity.  Tor  the  TSP,  the  r-optimal  heuristic  is  of 
order  Nr,  where  N  is  the  number  of  problem  nodes,  for  the  PUDP  with  a 
fixed  quality  of  service  parameter  Q,  the  algorithm  can  effectively  be 
used  in  order  NQr"^ .  Recall  that  Q  is  the  maximum  difference  between 
pickup  and  latest  delivery  for  a  customer.  The  reason  for  this  is  that 
it  is  highly  improbable  that  a  feasible  reconnection  pattern  exists 
when  the  difference  between  the  first  arc  removed  and  the  last  arc- 
removed  is  greater  than  or  equal  to  Q  arcs.  In  fact,  in  all  of  the 
problems  run,  not  once  was  a  better  solution  obtained  by  considering 
more  than  a  span  of  Q  arcs. 


4 


no 


Figure  10  depicts  a  span  of  Q  arcs  using  a  hypothetical  tour.  In 
this  figure,  letters  represent  the  points  visited  and  numbers  the 
connecting  arcs.  Suppose  Q=5  and  arc  3  is  the  first  arc  removed. 

Then  the  remaining  two  arcs  can  be  selected  from  among  arcs  4  through  7. 
The  question  of  feasibility  provides  an  insight  as  to  why  such  a 
reduction  in  computational  effort  is  possible.  Notwithstanding, 

Figure  11  shows  a  contrived  example  where  a  span  of  2Q-1  arcs  results 
in  a  feasible  reconnection  pattern.  In  the  figure,  Q=4;  the  points 
represent  cither  all  origins  or  all  destinations;  the  squares  the 
beginning  and  ending  points  for  the  six  arcs  in  the  span;  the  numbers 
inside  circles  the  position  number  for  a  point  in  the  tour;  and  the 
paired  numbers  the  applicable  time  windows.  Each  point  represents 
either  an  origin  or  a  destination.  Consequently,  the  time  windows 
contain  Q  stops,  not  the  Q+l  stops  which  define  the  time  window  for  the 
pair.  This  is  true  because  Q  is  effectively  the  difference  between 
the  latest  delivery  and  the  earliest  pickup  for  a  given  customer. 

Feasibility  and  reverse  tours.  Figure  12  is  a  copy  of  Figure  3, 
p.  28,  previously  discussed  in  Chapter  II.  Note  that  only  in 
reconnection  pattern  (2)  are  all  paths  (dark  lines  remaining  after  the 
arcs  are  removed)  traversed  in  the  same  forward  direction  as  is  the 
original,  feasible  tour.  In  all  others,  at  least  one  sequence  of  stops 
is  traversed  in  the  reverse  order.  In  patterns  (1)  and  (5),  both  paths 
are  traversed  in  reverse  order.  Because  of  the  precedence  requirements 
of  the  PUDP,  reverse  paths  place  very  stringent  requirements  on  the 


(a)  Initial  tour  with  3  arcs  removed 
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Figure  10.  example  of  a  2Q- 1  feasible 
reconnection  pattern  for 


120 


points  that  can  be  legitimately  included  in  these  paths.  Any  set  of 
points  to  be  traversed  in  reverse  order  cannot  feasibly  contain  an 
origin/destination  pair.  Also,  interchanging  the  relative  order  that 
the  disconnected  paths  are  to  be  traversed  restricts  the  permissablc 
location  of  a  corresponding  destination.  For  example,  in  patterns  (4) 
and  (5),  any  origin  among  the  points  between  4  and  S  must  have  its 
destination  at  or  beyond  point  6. 

Feasibility  and  timing.  Because  of  the  assumed  time  windows 
associated  with  each  customer  point,  many  reconnections  result  in  one 
or  more  time  windows  being  violated.  For  example,  if  any  of  the  stops 
from  points  2  to  3  are  already  sequenced  at  their  not-latcr-than  times, 
patterns  (2),  (3),  (4)  and  (5)  will  all  result  in  that  stop  being 
visited  at  a  later  time.  For  noncontrived  problems,  many  solution 
points  are  sequenced  at  or  near  one  end  or  the  other  of  their 
respective  time  windows.  Consequently,  slight  shifts,  imposed  by  the 
various  reconnection  patterns,  can  result  in  a  violation  of  the  time 
window.  For  this  reason,  a  span  of  Q  stops  was  sufficient  to  obtain 
the  best  5-optimal  tour  with  a  very  high  probability. 

Checking  feasibility.  A  proposed  tour  V  is  constructed  and  then 
checked  for  feasibility.  Because  any  tour  presented  to  the  3-optima  1 
algorithm  is  required  to  be  feasible,  it  is  only  necessary  to  check 
that  portion  of  Y  between  points  1  and  b  to  determine  if  the  entire  tour 
if  feasible.  The  check  insures  that  the  origin  precedes  its 
corresponding  destination  and  that  all  time  window  constraints  are 


satisfied.  When  a  feasible  improvement  is  found,  this  solution 
replaces  the  incumbent  and  the  algorithm  is  repeated.  When  no  improve¬ 
ment  can  be  found  for  an  entire  tour,  the  current  solution  is  the 
3-optimal  solution.  Despite  the  restrictions  imposed  on  the  reconnec¬ 
tion  patterns  by  the  precedence  and  time  window  constraints,  the 
5-optimal  algorithm  was  clearly  superior  to  all  other  heuristics 
tested . 


VI.  COMPUTATIONAL  ki.SULTS 

All  of  the  heuristics  discussed  above  were  tested  using  randomly 
generated  data.  In  addition,  the  ordering  !'  was  used  as  an  initial 
feasible  solution  for  the  5-optimal  algorithm.  Therefore,  the  3-optimal 
approach  was  used  both  as  a  tour  improvement  heuristic  and  in  conjunction 
with  one  of  the  tour  construction  heuristics  as  a  composite  heuristic. 

Heuristics  Verses  Optimali t y 

Data .  When  time  windows  are  present,  as  in  the  dial-a-ride  service 
problem,  optimal  solutions  arc  readily  available  which  can  be  used  to 
compare  various  heuristic's  effectiveness.  Table  (i,  on  p.  81,  contains 
the  not- latcr-than  times  used  for  defining  the  time  windows.  Table  11 
presents  the  solution  values  for  10  problem  instances.  The  instances 
are  different  in  that  the  location  of  each  point,  either  an  origin  point 
or  a  destination  point,  is  randomly  generated.  Consequently  this 
results  in  entirely  different  cost  matrices.  Costs  arc  again  computed 


using  the  rectilinear  metric,  ihc  data  in  Table  11  assume  a  quality  of 
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service,  Q,  of  five  stops  and  no  restriction  on  maximum  waiting  time, 

M.  Later,  varying  the  quality  of  service  parameters  is  examined. 

i'our  construction  comparison.  The  data  indicate  that  none  of  the 
tour  construction  heuristics  perform  particularly  well.  While  the 
average  performance  on  the  ten  problems  is  in  the  neighborhood  of  12T 
over  optimal,  there  are  several  cases  where  solutions  of  20T  or  more 
above  optimal  are  obtained.  Clearly,  there  is  nothing  to  indicate  that 
any  of  the  tour  constuction  heuristics  is  superior  to  another.  These 
results  are  interesting  in  that  the  Clarke-Wright  heuristic,  which  has 
shown  such  considerable  success  on  the  related  TSP  and  VRP,  performs  so 
poorly  on  the  more  heavily  constrained  l’UPP.  Also  noteworthy  is  that 
the  fairly  simple  greedy  procedure  produces  results  that  arc  on  an 
average  as  good  as  the  other  more  complicated  tour  construction 
heuristics.  An  explanation,  which  again  related  to  the  precedence 
relationship,  is  offered  beginning  on  page  129. 

Tour  improvement  and  composite  heuristics.  The  3-optimal  heuristic 
applied  to  the  ordering  P  unquestionably  produces  better  results  than 
any  of  the  tour  construction  heuristics.  In  fact,  these  results  are 
equal  to  those  of  the  composite  heuristics.  On  average,  the  3-optimal 
procedure  on  any  given  feasible  tour  produces  results  within  3"  of  the 
optimal  solution,  although  some  of  the  3-optimal  solutions  are  more  than 
10°o  above  optimal.  When  the  best  of  five  3-optimal  solutions  is  used 
the  results  arc  on  average  within  1°,  of  optimal.  Ihc  five  initial 
feasible  tours  are  randomly  generated.  These  results  are  considered 
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excellent,  but  the  computational  effort  is  five  times  as  great  as 
finding  the  3-optimal  solution  from  the  ordering  P.  There  does  not 
appear  to  be  any  rationale  for  using  a  composite  heuristic.  The  expense 
of  finding  the  initial  feasible  tour  by  one  of  the  tour  construction 
heuristics  is  too  great,  and  offers  no  apparent  advantage  over  the 
simpler  techniques. 

Varying  Service  Parameters 

The  results  discussed  above  assumed  a  quality  of  service  of  five 
stops  and  no  special  restriction  on  early  arrivals  (Q=M=5) .  lhe 
question  next  to  be  addressed  is  how  altering  these  parameters  affects 
the  efficiency  of  the  heuristics.  Tables  12  thru  15  present  solution 
data  for:  (1)  Q=M=11;  (2)  Q=ll,  M=6;  (3)  Q=7,  M=4 ;  and,  (4)  Q=5,  M=3 
respectively.  The  cost  data  and  not- later-than  times  for  the  destina¬ 
tions  are  identical  to  that  used  to  obtain  the  data  discussed  above. 
Therefore,  all  solution  differences  are  directly  attributable  to 
changes  in  the  service  parameters. 

Comments  on  results.  As  the  constraints  become  less  binding,  the 
tour  construction  heuristics  perform  even  more  poorly.  Solutions  of 
more  than  50”o  above  optimal  were  obtained.  The  3-optimal  heuristic 
also  showed  some  deterioration  in  effectiveness  as  the  constraints  were 
relaxed.  However,  the  average  results  were  still  within  5°o  of  the 
optimal.  As  before,  there  is  no  advantage  to  using  any  of  the  tour 
construction  heuristics  in  a  composite  mode.  The  3-optimal  algorithm 
applied  to  any  feasible  tour  performs  as  well  as  any  of  the  composites. 
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Large  Problems 

It  was  anticipated  that  the  basic  results  discussed  above  would 
carry  over  to  large  scale  problems.  Table  16  verifies  this  fact  for 
three  problems  with  N=91,  Q=7  and  M=4 .  The  best  of  the  five  3-optimal 
solutions  again  produced  excellent  results  that  were  within  2 *  of  the 
optimal  value. 

Analysis  of  Results 

The  precedence  requirement  is  deemed  to  be  the  primary  culprit 
which  renders  the  tour  construction  heuristics  ineffective.  Because 
each  of  these  heuristics  constructs  tours  which  maintain  either  actual 
or  relative  order,  the  order  becomes  locked  in.  Not  only  is  the  point 
being  sequenced  locked  into  a  given  relative  order,  its  corresponding 
pair-mate  is  also.  Only  in  the  pair  selection  heuristic  is  the  cost 
of  sequencing  the  pair-mate  considered  at  all  when  attempting  to 
sequence  a  given  point.  However,  based  on  the  empirical  data,  even 
explicitly  considering  such  costs  is  not  necessarily  effective. 

Figure  13  depicts  an  example  of  how  the  precedence  requirement 
might  lead  to  problems  with  the  Clarke-Wright  or  greedy  heuristics.  In 
panel  (a)  one  supposes  that  a  partial  tour,  depicted  by  the  dark  line, 
has  been  constructed  and  origin  i  is  determined  to  be  the  next  point  to 
be  added  to  the  tour.  Because  of  the  location  of  i's  corresponding 
destination,  j,  the  tour  must  at  sometime  retrace  itself  in  order  to 
visit  stop  j,  as  shown  in  panel  (b).  Stop  j  would  not  sequenced  until 
it  was  absolutely  necessary  to  do  so.  Panel  (c)  presents  another 
possible  partial  tour  which  is  clearly  much  cheaper. 


I 
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For  the  pair  insertion  heuristic.  Figure  14  presents  an  example  of 
the  possible  problems.  Panel  (a)  depicts  a  partial  tour.  Pair  number 
3  is  the  next  pair  to  be  added.  Panel  (b)  shows  the  insertion  of  0-^ 
and  Dj  into  the  partial  tour  in  the  prescribed  manner.  Panel  (c)  shows 
another  partial  tour,  which  is  obviously  superior  to  that  depicted  in 
panel  (b) .  Because  the  relative  order  of  Dj  preceding  Cb  in  panel  (a), 
the  configuration  in  panel  (c)  could  never  be  obtained,  since  the 
relative  order  there  has  Cb  preceding  Pj . 

Because  of  the  problems  created  by  the  precedence  requirement  in 
attempting  to  construct  a  tour  point  by  point  or  pair  by  pair,  it  is 
considered  doubtful  that  any  tour  construction  heuristic  could  be 
developed  that  is  superior  to  the  3-optimal  heuristic.  Further,  the 
3-optimal  heuristic  is  more  efficient  in  terms  of  computational  effort. 
Therefore,  any  further  attempt  to  develop  tour  construction  heuristics 
for  the  PUDP  is  not  recommended.  In  Chapter  VII  subject  areas  that 
appear  to  be  much  more  promising  for  further  research  are  suggested. 


ligure  14.  An  example  demonstrating  pa i r 
insertion  problems. 
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MUl.il  I’LL  VLII1CI.1.  l’ROBLl.M 

In  this  chapter,  the  algorithms  used  in  the  previous  chapters  to 
solve  the  single  vehicle  pickup  and  delivery  problem  are  extended  to 
the  multiple  vehicle  case.  The  multiple  vehicle  problem  is  more 
difficult  to  solve  than  is  the  single  vehicle  problem.  This  is  true 
both  for  exact  and  heuristic  solutions  to  the  problem.  Lxact  solution 
to  the  problem,  using  the  dynamic  programming  algorithm  detailed  in 
Chapter  IV,  is  first  discussed.  Then  some  of  the  heuristic  solutions 
detailed  in  Chapter  V  are  examined. 

I.  LXACT  SOLUTION  10  Till.  MUl.il  I'LL 
VLIIICLI.  I’ROIU.UM 

The  discussion  contained  in  Chapter  IV  suggested  that  dynamic 
programming  was  the  only  exact  solution  algorithm  that  appeared 
promising  for  optimal  solution  of  the  single  vehicle  I'UUl’.  Ihe  inherent 
complexity  of  the  problem  was  the  primary  reason  that  other  techniques 
were  determined  to  be  ineffective.  Because  the  multiple  vehicle  I’UUl’ 
is  even  more  complex  than  the  single  vehicle  version,  these  same 
arguments  appear  to  apply  even  more  strongly.  Therefore,  only  dynamic 
programming  was  considered  for  extension  to  the  multiple  vehicle  case. 
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Dynamic  Programming  extended  to  the 
Multiple  Vehicle  Problem 

Very  large  single  vehicle  PUUP's  were  solved  using  the  dynamic 
programming  algorithm  developed  in  Chapter  IV.  The  key  factor  that 
allowed  for  such  success  was  that  all  problem  constraints  were  assumed 
to  be  expressable  in  terms  of  stop  numbers.  The  state  space  could, 
therefore,  be  drastically  reduced  by  only  generating  feasible  combina¬ 
tions.  A  binary  representation  vector  was  used  to  identify  the  specific 
customer  points  that  comprised  each  of  these  vectors.  The  state  space 
cannot  be  as  drastically  reduced  in  the  multiple  vehicle  problem.  In 
fact,  given  a  reasonable  quality  of  service  parameter,  Q,  the  constrained 
problems  that  can  be  optimally  solved  are  nearly  identical  to  their 
unconstrained  couterparts.  The  reason  for  this  lies  in  the  manner  in 
which  the  state  vectors  are  generated  and  then  used  to  determined  the 
optimal  solution. 

Solution  on  an  expanded  network.  In  Chapter  II,  the  solution  of 
the  multiple  travelling  salesman  problem  on  an  expanded  network  was 
discussed.  The  same  technique  can  be  used  for  the  PUDP.  Ihc  only 
special  precaution  necessary  for  using  the  dynamic  programming  algorithm 
is  that  one  of  the  points  representing  the  depot  can  only  be  considered 
for  the  lead  position  if  the  set  S  contains  nothing  other  than  origin/ 
destination  pairs  an  1  possibly  points  representing  the  common  depot. 

This  insures  that  a  customer  is  not  picked  up  by  one  vehicle  and 
delivered  by  another. 
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Instances  of  interest.  Consider  the  situation  where  each  vehicle 


is  to  provide  service  to  an  equal  number  of  customers.  Therefore,  the 
number  of  customers  is  assumed  to  be  an  integer  multiple  of  the  number 
of  vehicles.  For  example,  if  there  were  five  vehicles,  then  there 
would  be  either  5,  10,  15,  20,  25.  .  .customers.  Given  this  scenario, 
those  values  of  k  where  the  lead  position  in  the  dynamic  programming 
state  vector  is  to  be  a  depot  point  are  well  defined.  Suppose  there 
were  25  customers  to  be  assigned  among  the  five  vehicles.  Recalling 
that  k  represents  the  total  number  of  points  considered  for  inclusion, 
the  lead  positions  would  be  the  11th,  22nd,  33 rd,  etc.  There  arc  ten 
positions  to  accomodate  five  pairs  prior  to  the  first  depot  point  (one 
route);  21  positions  to  accomodate  two  sets  of  five  pairs  and  the  first 
depot  point  prior  to  the  second  depot  point  (two  routes);  etc. 

Although  this  problem  could  be  solved  on  the  expanded  network,  as 
suggested  above,  a  simpler  method  is  available  for  this  instance  that 
requires  considerably  less  storage  than  docs  solution  on  an  expanded 
network.  Unfortunately,  even  the  simpler  method  requires  too  much 
storage  for  the  solution  of  problems  with  more  than  about  10  origin/ 
destination  pairs. 

Storage  requirements  for  expanded  network  solution.  Suppose  only 
the  precedence  constraints  arc  binding.  Let  b  represent  the  number  of 
or i gin/ dost i nation  pairs  assigned  to  a  vehicle  and  recall  that  n  is 
the  number  of  or i gin/dcst inat i on  pairs.  In  the  example  above  with  five 
vehicles,  n  =  25  and  b  =  5.  These  data  will  he  used  at  each  step  to 
illustrate  the  computations. 
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First,  consider  the  number  of  storage  locations  for  k  values 
directly  preceding  k  =  2b  +  1  or  11,  which  is  the  point  at  which  one  oi 


the  depot  points  is  first  considered  for  the  lead  position.  For 
k  =  2b  or  10  there  are 

n- 1  ' 

n  (0.1) 

b-1 

V. 

storage  locations  required  for  the  f10  values.  The  lead  position  must 
always  be  one  of  the  n  origins.  Expression  (0.1)  is  based  on  n  possible 
origins  visited  immediately  upon  departure  from  the  depot  and  the 
n-1  taken  b-1  ways  of  selecting  the  remaining  customers  to  be  serviced 
by  one  of  the  vehicles.  Thus,  the  example  data  result  in 

'24  ' 

25  =  265,050 

4  J 

required  locations.  For  k=2b-l  or  9  and  an  origin  in  the  lead 
position,  the  storage  required  is 

n-1  ( n-b+ 1 

n  .  (0.2) 

b-2  J  {  \ 

In  terms  of  the  example,  this  represents  25  choices  of  t he  lead  origin, 
24  remaining  pairs  taken  3  at  a  time,  and  one  choice  among  the  21 
remaining  destinations.  This  destination  is  not  a  match  to  one  of  the 
origins  being  considered,  but  would  match  the  origin  selected  subse¬ 
quently  at  the  k  =  10  iteration.  The  value  of  this  computation  is 


1,002,600.  If  a  destination  were  the  lead,  the  required  storage  would 
be  the  same  as  that  given  by  expression  (0.1).  In  the  example,  the 
required  storage  would  be  1,002,000  +  205,050  =  1,528,250  for  k  =  9 . 
Considering  k  =  10  in  this  same  example,  the  storage  required  for  just 
the  fi^  values  would  be  in  excess  of  08  million  storage  locations. 
Clearly,  practical  solution  to  the  multiple  vehicle  PUD!’  on  an  expanded 
network  is  limited  to  rather  small  problems. 

Alternative  Solution  Using 
Dynamic  Programming 

In  the  above  discussion,  once  all  of  the  f_.  ,  or  f , ,  values  are 

2b+l  11 

computed,  the  data  necessary  to  solve  the  problem  is  at  hand.  This  is 
so  since  the  f 2t>  + 1  values  represent  the  cost  of  all  feasible  tours  of 
the  n  customers  taken  b  at  a  time.  In  the  example,  the  f j j  values 
contain  data  on  all  combinations  of  customers  taken  5  at  a  time,  (liven 
5  vehicles  being  used,  the  optimal  solution  is  found  by  taking  the  5 
f 1 1  values  whose  sum  is  minimal  and  whose  combined  individual  state 
vectors  indicate  that  service  is  provided  to  all  customers.  In  general, 
this  can  be  accompl i shed  by  a  direct  comparison  of  all  of  the  fq,+j 
values  and  their  binary  representation  vectors.  The  FORTRAN  routine 
to  accomplish  this  requires  V  nested  DO  loops,  where  V  represents  the 
number  of  vehicles.  bach  loop  corresponds  to  one  of  the  vehicles,  and 
an  assignment  of  customers  to  that  vehicle  comes  from  the  identification 
vector.  Suppose  V  =  2  and  the  customers  set  is  { 1  ,  2 , 5,  -1 , 5 , 0  !  . 

Therefore,  three  customers  must  he  assigned  to  each  vehicle.  If  the 
identification  vector  for  the  first  loop  associated  with  a  specific 
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f2  value  identifies  the  set  {1,3,4},  while  the  one  for  the  second 
identifies  {2,5,6},  their  combination  results  in  service  to  all  of  the 
customers.  Furthermore,  if  the  sum  of  the  two  f -j  values  is  minimal  with 
respect  to  all  other  feasible  combinations,  tiie  two  optimal  routes  can 
be  solved  for  individually,  each  as  a  single  vehicle  problem,  thus 
solving  the  two  vehicle  problem. 

The  range  of  the  exterior  loop  is  over  the  number  of  f-,j +,  values. 
For  the  earlier  example  above,  there  would  be  five  loops,  each  with  an 
initial  range  of  approximately  265,650.  Let  F  =  { f 2b+ 1  vaiucs^-  At 
first  look,  it  might  appear  that  this  routine  requires  on  the  order  of 
I F | V  comparisons  or  265,650s  ±  1.323  x  1027  comparisons .  However,  the 
inner  loops  need  not  be  entered  unless  there  is  a  possibility  that  the 
optimal  solution  is  contained  within. 

Requirement  for  no  overlapping  customers.  The  representation 
vectors  identify  the  customers  serviced  by  each  one  of  the  |l:|  indi¬ 
vidual  tours.  If  the  same  origin/destination  pair  is  known  to  be  in 
any  two  of  the  loops,  then  there  is  no  way  that  a  feasible  solution  can 
exist.  The  reason  for  this  is  that  eacli  loop,  in  effect,  assigns 
exactly  b(N/V)  customers  to  a  vehicle.  If  the  same  customer  is  assigned 
to  more  than  one  vehicle,  it  must  be  true  that  at  least  one  customer  is 
not  assigned  to  any  vehicle.  Let 

Rj  =  representation  vector  of  the  ith  individual  tour.  {6.3} 

Suppose  V  >  3,  then  if 


A. 
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R.  fl  R.  t  0  (9-4J 

i  3 

for  the  first  two  loops,  the  second  loop  can  he  immediately  incremented 
since  at  least  one  origin/destination  pair  is  contained  in  both  partial 
tours.  In  terms  of  the  previous  example  where  n=25,  V=5  and  b  =  5, 
if  one  customer  pair  showed  up  in  the  intersection,  the  first  two  loops 
provide  assignments  for  only  9  customers.  The  best  that  the  three 
remaining  loops  could  do  is  provide  assignments  for  IS  more.  Conse¬ 
quently,  one  customer  is  not  assigned  to  any  vehicle,  and  the  solution 
is  thereby  infeasible.  If  the  intersection  is  null,  a  feasible  complete 
tour  can  be  found  and  the  third  loop  is  entered.  Let  Rk  be  a  vector 
for  this  third  loop.  Now,  if 

(Ri  n  Rk)  U  (Rj  n  Rk)  /  0  (6.S) 

the  third  loop  is  incremented  and  another  Rk  selected.  Otherwise  the 
fourth  loop  is  entered,  finally,  the  Vth  loop  is  entered,  a  representa¬ 
tion  vector  is  found  such  that  all  customers  are  serviced.  This  repre¬ 
sents  a  feasible  solution  and  the  cost  of  this  solution,  represented 
by  the  V  ^2b+l  va4ues>  4s  computed.  If  the  value  of  this  solution  is 
better  than  the  incumbant,  it  becomes  the  new  solution.  The  innermost 
loop  is  then  incremented  until  another  solution  is  found  or  until  the 
loop  parameter  is  exhausted.  These  basic  steps  are  continued  until  the 
lowest  cost  feasible  complete  tour  is  found. 


FORTRAN  intersection.  Graves  (20)  provides  techniques  for 

performing  logical  operations  in  FORTRAN  without  recourse  to  assembly 

level  routines.  His  approach  is  an  order  of  magnitude  quicker  than 

using  an  assembly  level  routine  since  it  eliminates  several  storage  and 

retrieval  operations  for  register  values  as  the  program  shifts  control 

from  the  main  routine  to  the  subroutine  and  back  again.  The 

LQUIVALLNCli  statement  is  the  key  as  shown  in  the  following  algorithm 

which  computes  the  binary  intersection  A  of  two  representation  vectors 

B  and  C  and  prints  out  the  decimal  equivalent: 

INTLGliR  A ,  B ,  C 

L0GJCAI.*4  LA,LB,1,C 

FqiJIVAI. i:\ci:  (A,  LA ),( B,  LB)  Cj 

I  A  =  LB.AND.LC 

PRINT  A 

Thus,  if  A  =  0,  the  intersection  is  null.  If  A  f  0,  two  or  more 
origin/destination  pairs  arc  continued  in  both  B  and  C.  Identification 
of  which  pairs  are  redundant  could  be  found  by  converting  the  decimal 
value  of  A  into  its  binary  equivalent. 

Limitation  on  Problem  Sice 

As  mentioned  above,  exact  solution  to  problems  with  more  than  about 
ten  origin/destination  pairs  is  not  practical  in  the  case  of  the 
multiple  vehicle  PUUP.  lime  windows  can  still  be  used  to  determine 
those  points  which  are  valid  condidates  for  the  lead  position  or  which 
are  eligible  to  be  used  to  complete  a  given  state  vector.  They  cannot 
be  used  to  require  points  to  be  included  in  the  completion  of  a  state 
vector.  A  point  which,  if  not  sequenced  would  preclude  anv  feasible 
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solution  in  the  single  vehicle  problem,  can  be  assigned  to  any  of  the 
V  vehicles  in  the  multiple  vehicle  case,  lor  the  alternative  solutions 
approach  discussed  above,  there  is  never  an  instance  where  a  point  must 
be  included.  The  feasible  state  space  is  also  larger  due  to  the  fact 
that  the  candidate  set  for  the  lead  position  is  approximately  V  times 
greater  than  it  is  for  t he  single  vehicle  problem.  Consequently,  the 
number  of  feasible  states  for  each  value  of  k  and  choice  of  lead  is 
considerably  greater. 

Computational  Results 

Table  17  presents  typical  results  for  various  sited  problems.  Only 
the  precedence  constraints  were  binding.  The  largest  of  these  problems 
actually  solved  involved  only  nine  origin/dcst ina t ion  pairs  and  three 
vehicles.  This  equates  to  a  21  point  multiple  travelling  salesman 
problem  on  an  expanded  graph.  Attempting  to  add  one  more  customer  per 
vehicle  resulted  in  exceeding  a  five  minute  execution  time  limit.  In 
fact,  the  limit  was  exceeded  during  the  computation  ol  one  of  the 
9,900  fc;  values.  This  means  the  problem  was  less  than  one  half  solved 
at  the  time  of  termination,  since  storage  required  up  to  and  including 
f 5  is  only  14,994  locations. 

Constrained  problems  produced  very  similar  results.  The  12 
customer,  .>  vehicle  problem  was  too  large  to  solve.  Using  a  value  of 


Q=5,  the  time  windows  for  the  remaining  problems  are  very  similar  to 
those  for  an  unconstrained  problem.  Consequently,  t he  number  of  feasible 
states  generated  is  not  significantly  reduced.  Table  18  provides  a  side 


I 


Lj  TABLE  17 

L'l  SELECTED  OPTIMAL  RESULTS  FOR  THE.  MULTIPLE  VEHICLE 

L  !  PICKUP  AND  DELIVERY  PROBLEM 


NUMBER 

OF 

VEHICLES 

NUMBER 

OF 

CUSTOMERS 

N 

STORAGE 

REQUIRED 

RUN 

TIME 

(SEC) 

2 

4 

9 

65 

3.1 

2 

6 

15 

667 

4.3 

2 

8 

17 

6,337 

74.53 

3 

6 

13 

157 

3.2 

3 

9 

19 

2,566 

22.4 
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12 

25 

39,414 
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by  side  comparison  for  the  9  customer,  3  vehicle  problem.  While  the 
unconstrained  problem  required  2,566  storage  locations,  the  constrained 
counterpart  requires  2,128  locations.  Consequently,  optimal  solution 
to  any  real  world  multiple  vehicle  problem  is  not  deemed  practical. 

The  state  space  becomes  too  large,  just  as  it  does  in  the  travelling 
salesman  problem. 

II.  IILURISTIC  SOLUTION  TO  Till:  MULTI  I'LL 
VLHICLL  PROBLLM 

The  multiple  vehicle  1’UUi’  was  shown  in  the  preceding  section  to  he 
more  difficult  to  solve  optimally.  It  is  also  more  difficult  to  obtain 
a  good  solution  by  means  of  a  heuristic,  although  a  precise  measure  of 
"goodness"  is  difficult  because  only  small  problems  can  be  solved 
optimally.  For  the  single  vehicle  problem,  the  3-optimal  heuristic 
produced  solutions  generally  within  a  few  percentage  points  of  the 
opt ima 1  solution.  As  will  be  seen  below,  the  3-optimal  technique  is 
severely  limited  because  of  the  precedence  relationship  in  the  multiple 
vehicle  problem. 

The  key  to  good  solutions  in  the  multiple  vehicle  problem  is 
determining  which  vehicle  customers  should  lie  assigned  to.  Once 
customers  are  assigned  to  a  given  vehicle,  the  problem  reduces  to  V 
single  vehicle  problems  for  which  optimal  or  very  good  solutions  arc 
readily  attainable.  Several  attempts  were  made  to  determine  a  good 
method  to  make  such  assignments.  The  results  are  somewhat  discouraging 
as  will  be  seen  later.  Consequently,  the  multiple  vehicle  problem 
appears  to  remain  relatively  unsolved. 
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Clarke-Wright  and  Pair  Insertion 

Two  of  the  tour  construction  heuristics  discussed  in  the  last 
chapter  were  the  Clarke-Wright  savings  heuristic  and  the  pair  insertion 
heuristic.  both  of  these  heuristics  were  expensive  in  terms  of  compu¬ 
tational  effort  and  performed  relatively  poorly.  Other  heuristics 
produced  much  better  results.  Therefore,  neither  of  these  approaches 
was  pursued  as  a  possible  solution  to  the  multiple  vehicle  problem. 

Orcody  Heuristic 

The  greedy  heuristic  is  also  a  tour  construction  heuristic,  but  is 
much  less  expensive  computationally.  The  heuristic  also  lias  the 
advantage  that  it  is  easily  understood  and  can  often  be  implemented 
manually.  Therefore,  the  greedy  heuristic  was  selected  to  determine  how 
well  a  tour  construction  heuristic  performed  on  the  multiple  vehicle 
problem.  The  results  were  surprisingly  good. 

Construction  concept.  There  arc  two  ways  that  tours  could  be 
constructed  using  the  nearest  neighbor  concept.  One  would  be  a 
sequential  approach  where  vehicles  arc  dispatched  one  at  a  time.  In 
this  case,  1 / V  of  the  customers  would  be  assigned  to  the  first  vehicle 
before  the  second  is  dispatched.  The  second  approach  dispatches  all 
vehicles  at  the  same  time  and  builds  routes  simultaneously.  The  latter 
approach  was  used  for  the  multiple  vehicle  I’llDI*,  since  it  produced 
individual  routes  of  nearly  equal  length.  The  former  method  tends  to 
produce  at  least  one  relatively  short  route  and  at  least  one  relatively 
long  route  which  may  be  undesirable. 


Description  of  algorithm.  Initially,  the  V  nearest  origins  to  the 
depot  which  satisfy  the  time  window  constraints  are  assigned  one  each  to 
the  V  vehicles.  Then  a  second  stop  is  determined  for  all  the  vehicles, 
followed  by  the  third,  fourth,  etc.  A  candidate  set  is  used  to 
identify  all  unsequcnced  points  that  satisfy  the  time  window  constraints. 
For  each  vehicle,  the  next  nearest  feasible  neighbor  from  this  candidate 
set  is  added  to  the  partial  tour,  provided  a  check,  made  to  insure  that 
adding  the  nearest  neighbor  will  result  in  a  feasible  tour,  is 
successful . 

Feasibility  check.  Insuring  that  points  added  allow  for  a  feasible 
set  of  individual  tours  is  essential  for  the  multiple  vehicle  problem, 
just  as  it  was  for  the  single  vehicle  problem.  The  construction  or 
flushing  out  of  the  partially  completed  tours  is  a  bit  more  complex  in 
the  multiple  vehicle  problem.  The  nondecrcasing  ordering  I'  is  used, 
just  as  in  the  single  vehicle  case.  Ilie  next  element  of  P,  not  already 
sequenced  by  the  greedy  selection,  is  placed  on  the  vehicle  with  the 
fewest  elements  unless  it  is  a  destination.  Destinations  must  be 
placed  on  the  same  vehicle  as  their  corresponding  origins.  Fach  of 
the  V  individual  tours  is  then  checked  to  verify  that  it  is  feasible. 
Unless  all  V  of  the  tours  are  feasible,  the  nearest  neighbor  being 
considered  for  addition  to  one  of  the  tours  is  not  added.  Rather  it  is 
removed  from  the  candidate  set  for  that  vehicle  and  another  candidate 
(the  nearest  remaining  feasible  neighbor)  is  selected  and  the  feasibilit 
check  repeated. 
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Optimal  individual  tours.  Once  all  V  tours  have  been  constructed, 
it  is  often  possible  to  determine  the  optimal  routing  for  each  tour. 

This  is  possible  when  all  problem  constraints  are  expressable  in  terms 
of  the  vehicles'  stop  numbers  as  discussed  in  Chapter  IV.  The  greedy 
solution  identifies  which  pairs  are  on  each  vehicle,  thus  defining  V 
single  vehicle  problems.  If  the  constraints  arc  expressed  in  a  manner 
such  that  the  dynamic  programming  algorithm  can  not  be  used,  or  if 
the  problem  is  too  large  for  the  dynamic  programming  algorithm,  the 
3-optimal  heuristic  can  be  applied  to  the  greedy  solution.  These 
techniques  were  applied  to  all  of  the  problems  studied  herein. 

Greedy  look  ahead.  The  basic  greedy  algorithm  is  myopic  in  that 
one  selects  the  nearest  neighbor  without  thought  of  where  the  tour  will 
next  go.  The  precedence  requirement  of  the  ITIDI'  necessitates  subsequent 
travel  to  the  as  yet  unscquenced  destinations.  It  was  conjectured  that 
by  considering  (minimizing)  the  distance  over  the  next  two  points  that 
a  better  solution  could  be  obtained,  especially  in  terms  of  which  pairs 
were  assigned  to  which  vehicle.  The  basic  greed)'  algorithm  is  an  order 
Z  algorithm,  where  Z  represents  the  cardinality  of  the  current  candidate 
set.  The  look  ahead  greedy  is  approximately  an  order  Z'  algorithm. 
Results  indicated  rio  appreciable  difference  between  the  greedy  and  the 
look  ahead  greedy.  Therefore,  the  look  ahead  version  was  not  tested 


further . 


The  heavier  computational  effort  could  not  he  justified. 


r-Optimal  Heuristic 

Limitations.  The  results  of  Chapter  Y  clearly  showed  the 
superiority  of  the  3-optimal  heuristic  when  compared  to  any  of  the  tour 
construction  heuristics.  However,  the  precedence  constraint  nullifies 
much  of  the  effectiveness  of  the  r-optimal  heuristic  when  applied  to 
the  multiple  vehicle  problem.  The  solution  presented  to  the  5-optimal 
algorithm  is  on  an  expanded  network.  Since  a  very  good  (if  not  optimal) 
solution  can  be  obtained  for  the  individual  tours,  interchanging  arcs 
between  tours  and  thus  creating  a  change  in  the  customers  assigned 
to  a  given  vehicle  is  what  is  desired.  The  problem  stems  from  the 
necessity  that  both  the  origin  and  destination  be  on  the  same  vehicle. 
For  the  pickup  and  delivery  problem,  this  requirement  can  only  be  met 
when  the  vehicles  are  empty ,  as  is  shown  below. 

Reconnect  ion  pattern.  The  crux  of  the  problem  can  be  readily  seen 
by  considering  a  2-optimal  example.  In  Figure  13(a),  the  two  indi¬ 
vidual  tours,  each  with  an  arc  removed,  are  shown  side  by  side.  ihe 
four  (j*)'s  all  represent  the  common  depot,  and  the  individual  tours 
begin  at  the  bottom  and  proceed  upward.  Fach  letter  represents  the 
set  of  points  between  the  depot  and  a  removed  arc.  Figure  1  a  (. b  ) 
depicts  the  only  available  reconnection  pattern.  Let  i  and  j  represent 
an  origin  and  its  correspond) ng  dost i nat ion  respectively.  For  the 
configuration  in  Figure  15(a)  to  be  initially  feasible  the  following. 


conditions  must  hold: 


0  G  O  0 


Original  tours  with  (bj  New  tours  after 

two  arcs  removed  reconnection 


I'ijpirc  15.  I'ictoral  of  2-optim.il  reconnection 
patten)  for  the  multiple  veld  do 
prob  1  ein . 


i eA  =>  j  e ( AU  B ) 


(6.6) 


it-C 

=  >  j  i-  (C1JL)) 

(0.7) 

j  E.B 

=  >  it(Al)B) 

(6.8) 

and 


j  lD  =>  ic  (CUL)J  .  (6.9) 

Clearly,  if 

it.A  and  jcB  or  if  ieC  and  jel) 

removal  of  either  of  these  arcs  could  never  result  in  a  feasible 
solution,  as  reconnection  would  result  in  the  origin  assigned  to  one 
vehicle  and  the  destination  to  another.  Only  when  the  vehicle  is  empty 
can  an  arc  feasibl)  be  removed.  In  order  to  maintain  an  equal  number  of 
customers  serviced  by  each  vehicle,  both  vehicles  must  become  empty  at 
the  same  stop  number.  The  likelihood  of  this  occurring  is  rather  small 
based  on  the  empirical  evidence.  liven  when  such  an  occurancc  does 
materialized,  the  interchange  is  seldom  profitable. 

hxtension  to  5-optimal.  Given  the  restriction  that  each  vehicle 
must  provide  service  to  the  same  number  of  customers,  there  is  only  one 
way  three  arcs  can  be  feasibly  removed.  Removal  of  two  arcs  from  one 
tour  alters  the  number  of  customers  served  by  a  given  vehicle.  There¬ 


fore,  the  arcs  must  come  from  each  of  three  separate  tours,  all  at  the 
same  vehicle  stop  number,  and  with  all  vehicles  empty  prior  to  the  break. 
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The  likelihood  of  this  occurring  is  even  less  than  with  a  two  arc 
removal.  Preliminary  results  suggested  the  approach  was  not  effective 
in  solving  the  multiple  vehicle  PUDP  and  was  not  pursued  further.  This 
is  but  another  example  of  how  the  precedence  relationship  of  the  PUUP 
severely  limits  the  effectiveness  of  heuristics  that  otherwise  work  well 
on  related  problems. 

Pair  Selection 

The  pair  selection  heuristic  identifies  the  profitabil ity  of 
assigning  each  pair  of  origin/destination  pairs  to  the  same  vehicle. 
Using  a  savings  value  computed  for  each  pair  of  pairs,  a  group  of 
customers  is  eventually  identified  for  each  of  the  Y  vehicles.  Once 
the  V  groups  are  identified,  the  problem  is  reduced  to  V  single  vehicle 
problems,  each  solved  either  optimally  or  by  the  5-optimal  heuristic. 

Concept .  The  motivation  for  the  pair  selection  heuristic  lies  in 
taking  advantage  of  the  requirement  to  move  from  an  origin  to  its 
corresponding  destination.  If  two  customers  cun  be  served  as  cheaply 
as  one,  or  nearly  so,  it  appears  advantageous  to  do  so.  figure  16  shows 
two  origin/destination  pairs  tiiat  can  both  be  serviced  for  the  same 
cost  as  serving  the  first  alone.  Therefore,  the  savings  value  is 
computed  to  determine  how  much  can  be  saved  by  serving  two  customers 
on  the  same  route.  This  value  is  termed  SAW 

SAV  formulas.  For  any  pair  of  customers  taken  by  themselves  there 


are  only  six  ways  in  which  they  could  be  serviced.  Suppose  both 
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Figure  16.  Example  of  how  two  customers  can 
be  served  at  the  cost  of  one. 
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customers  are  considered  separately,  and  the  only  cost  for  each  is  that 
from  its  origin  directly  to  its  destination.  How  much  could  be  saved 
by  having  both  customers  on  the  same  route?  The  savings  is  taken  to  be 
the  cost  of  linking  the  two  origin/destination  pairs  minus  the  cost  of 
separate  service.  The  SAV  value,  which  may  be  either  positive  or 
negative,  is  taken  to  be  the  minimum  cost  of  each  of  the  six  possible 
service  patterns.  Table  19  gives  the  six  different  reconnection  or 
service  patterns  and  the  savings  formulas  in  IOUTRAN  symbology  for  each. 
The  blank  spaces  result  from  costs  that  appear  with  both  a  positive 
and  a  negative  sign  in  the  same  formula  and  are  thus  cancelled  out. 

Selection  procedure.  The  SAV  value  that  is  the  most  negative 
(least  positive)  identifies  the  first  two  origin/dcstination  pairs  to 
be  assigned  to  a  vehicle.  The  next  smallest  SAV  value  identifies  the 
next  two  pairs.  Provided  neither  is  assigned  to  the  first  vehicle,  they 
are  assigned  to  the  second  vehicle.  If  one  had  been  assigned  to  the 
first  vehicle,  the  other  would  be  also.  This  would  leave  the  first 
vehicle  with  three  assigned  customers.  With  only  precedence  constraints 
binding,  this  process  is  continued  until  at  least  one  customer  (origin/ 
destination  pair)  is  assigned  to  each  of  the  V  vehicles. 

Subsequent  pair  assignment.  Once  all  V  vehicles  have  been  assigned 
at  least  one  customer,  the  process  becomes  slightly  more  complex.  Let 
I  and  J  represent  the  next  pair  to  be  considered  based  on  the  SAV(1,J) 
being  the  minimum  of  the  remaining  values.  If  1  lias  already  been 


assigned  but  .J  has  not,  .1  is  assigned  to  the  same  vehicle  as  I  is, 
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TABLE  19 

PAIR  SELECTION  PATTERNS  AND 
SAVINGS  FORMULAS 


PATTERN 

SAVING  FORMULA 

(1) 

I-J-ID-JD 

C  (  I , J) +C ( J ,  I D  )  +C  ( I  D,  .ID)  -  C  (  J  ,  1D.)-C(J,JD) 

(2) 

I-J-JD- ID 

C(I,.J)  +C  (.ID,  I D)  -C  ( I ,  ID) 

(3) 

I-1D-J-JD 

C(1D,J) 

(4) 

J  - 1  -  JD- ID 

C(J,  1)+C(I,JD)+C(JD, I D) -C  C I , ID)-C(J.JD) 

(5) 

J-I-ID-JD 

C(.J,I)  +C  ( I D ,  JD)  -C  (J  ,.JD ) 

{  ('  1 

J-JD-I-ID 

C (JD, I ) 

NOTE:  I,  the  first  origin 

J,  the  second  origin 
C(',*),  Cost  or  distance  from  •  to  * 
ID,  I ' s  destination 
JD,  J's  destination 


i 
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provided  the  vehicle  does  not  already  have  its  full  complement  of 
customers.  If  I  has  been  assigned  to  one  vehicle  and  J  to  another,  the 
pair  is  already  taken  care  of  by  previous  assignments.  If  neither  can 
be  matched  with  a  vehicle  at  this  time,  the  pair  is  held  until  one  of 
the  customers  is  assigned  to  one  of  t lie  V  vehicles.  Then  the  other  is 
assigned  to  the  same  vehicle  provided  that  the  first  customer  did  not 
complete  the  vehicle's  complement.  Also,  when  time  windows  are  present, 
it  is  necessary  to  insure  that  each  assignment  can  lead  to  feasible 
individual  tours. 

feasibility  of  assignments.  Insuring  feasible  tours  given  the 
presence  of  time  windows  has  been  a  critical  factor  in  all  of  the  I’UUP 
heuristics  examined.  The  pair  selection  procedure  is  no  exception. 

The  procedure  parallels  that  used  in  both  the  single  vehicle  cases 
and  the  multiple  vehicle  case  by  greedy  selection.  V  individual  tours 
are  constructed  using  the  ordering  P  and  then  each  tested  for 
feasibility.  The  customer  or  customers  being  considered  for  service 
by  a  given  vehicle  are  assigned  to  that  vehicle  only  if  all  V  of  the 
tours  are  feasible. 

Interchange  Heuristic 

Concept .  Given  V  feasible  individual  tours,  the  interchange 
heuristic  attempts  to  identify  two  customers  on  different  tours  that 
can  be  switched  so  as  to  reduce  the  combined  cost  of  the  two  respective 
tours.  This  is  accomplished  by  first  determining  for  each  customer  the 
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penalty  associated  with  being  on  the  tour  it  is  on  and  not  on  one  of 
the  other  tours,  and  then  matching  the  penalties  to  find  those  that 
appear  to  be  good  candidates  for  a  switch. 


Penalty  computation.  The  penalty  for  origin  i  and  its 
corresponding  destination  ID  on  tour  k  not  being  on  tour  k'  is  defined 
as  the  difference  between  the  cost  of  connecting  to  this  pair  in  its 
present  tour  and  the  minimum  cost  to  connect  this  pair  to  tour  k'. 

Let  IP(IDP)  and  I F ( IDF)  represent  the  stops  directly  preceding  and 
directly  following  I  (ID)  in  its  present  tour.  Also  let  .1  and  J’  with 
J  preceding  J'  represent  the  closest  points  to  I  and  ID  on  tour  k' 
respectively.  Then  the  penalty  can  be  computed  by 

min  [C(IP,  I)  +C(  1 ,  Ii:J  +C(  I  DP,  ID)  +C(  ID,  I  DP  ) ; 

C ( IP, I) +C ( 1 , ID) +C ( ID , 1DP) ] 

-2 [C ( I , J J  +C ( ID, J ' ) ]  (0.10) 

Large  penalty  values  suggest  it  may  be  better  to  have  the  pair  on  the 
other  tour. 

Switching  pairs.  tt  would  not  be  practical  to  attempt  all  of  the 
possible  switches  of  two  pair.  Therefore,  switches  of  the  pairs  com¬ 
prising  the  highest  four  or  five  penalty  values  were  attempted.  "I ho 
switch  effects  a  change  in  the  composition  of  each  vehicle's  customers. 
Therefore,  the  new  configurations  are  optimally  solved  for  both  indi¬ 
vidual  tours.  If  the  combined  tour  cost  is  less  after  the  switch,  the 
switch  becomes  permanent  and  the  process  repeated. 
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interchange  results.  The  heuristic  was  tried  on  several  problems 
of  varying  size  and  complexity.  The  heuristic  very  seldom  identified 
a  profitable  switch  when  the  initial  tours  were  the  optimal  routing 
obtained  from  either  the  greedy  or  the  pair  selection  procedures. 
Consequent ly ,  this  heuristic  was  eliminated  from  further  consideration. 

Computational  Results 

Because  it  is  much  more  difficult  to  obtain  optimal  solutions  for 
the  multiple  vehicle  I'U l)P  than  it  is  for  the  single  vehicle  version, 
only  in  the  case  of  very  small  problems  can  the  heuristics  be  compared 
to  the  optimal  solution.  Tor  larger  problems,  results  are  limited  to 
a  relative  comparison.  As  in  the  single  vehicle  case,  the  impact  of 
varying  the  service  paramenters  is  also  investigated. 

Heuristics  verses  optimal.  Table  20  presents  results  from  five 
sample  problems.  These  data  indicate  that  the  pair  selection  technique 
is  very  effective  for  problems  of  this  size.  The  largest  deviation 
from  optimal  was  less  than  S°o.  The  data  also  suggests  that  the  greedy 
solution  followed  by  optimally  sequencing  the  individual  routes  is 
not  especially  effective.  Deviations  of  greater  than  20"o  arc  noted. 
However,  as  will  be  seen  later,  this  greedy  procedure  is  still  superior 
to  a  random  assignment  of  the  customers  among  the  vehicles.  Also  noted 
below  is  the  fact  that  on  larger  problems  with  time  windows  present, 
the  greedy  selection  heuristic  demonstrates  equal  or  slightly  better 
performance  when  compared  to  the  pair  selection  procedure. 


4 


TABLE  20 

COMPARISON'  OF  MULTIPLE  VEHICLE  HEURISTIC 
TO  THE  OPTIMAL  SOLUTION  VALUES 


N 

V 

OPTIMAL 

GREEDY 

GREEDY 

SELECTION 

PAIR 

Sl.I.I.Ci  ION 

17 

2 

604 

918 

7  34 

<>04 

17 

2 

688 

74  2 

730 

704 

17 

2 

530 

746 

598 

5  5  <-> 

13 

3 

802 

1 096 

974 

802 

13 

3 

704 

926 

808 

808 

160 


Comparison  of  greedy  atul  pair  selection  on  urn ons t ra i ned  problems. 
When  there  are  no  time  window  constraints  present,  the  pair  selection 
procedure  generally  produces  better  results  than  does  the  greedy 
selection  procedure.  Hie  data  in  Table  21  verifies  this  for  ten 
problems  with  16  customers  (N'=33)  and  four  vehicles.  Only  in  one 
instance  (problem  8)  was  the  greedy  pair  solution  better  than  the  pair 
selection  solution.  Only  a  relative  comparison  is  possible  since 
\  =  33  and  V=4  is  too  large  to  be  optimally  solved  by  dynamic  programming. 
Consequently,  these  results  may  or  may  not  be  anywhere  near  optimal. 

The  questions  of  whether  these  selection  procedures  are  superior  to  a 
haphazard  selection  procedure  remains  to  be  answered. 

Selection  comparison  on  unconstrained  problems.  Table  22  presents 
ten  solutions  to  each  of  two  problems  with  only  precedence  constraints 
binding.  The  first  and  second  solutions  for  each  problem  are  the 
solutions  obtained  from  the  greedy  selection  and  pair  selection 
heuristics  respectively.  Solutions  three  through  ten  represent  a  random 
assignment  of  four  customers  to  each  of  the  four  vehicles.  in  all 
cases,  each  vehicle  is  individually  optimally  routed  using  the  dynamic 
programming  algorithm.  Unquestionably,  both  the  greedy  and  pair 
selection  procedures  produce  better  results  than  just  randomly  or 
haphazardly  assigning  customers  to  the  vehicles.  It  is  also  evident 
that  the  key  to  a  good  solution  for  the  multiple  vehicle  I’UI'P  is  in  the 
assignment  of  customers  to  the  vehicles.  This  is  especially  true  for 
those  instances  that  allow  for  an  optimal  sequencing  of  the  assigned 
pairs.  Such  an  assignment  is  complicated  since  there  are 


TABU;  22 

COMPARISON  01  SOLUTION  VALIIPS  RASPI'  ON 
I)  I  II LRPN )  SI:  I. PC  I  I  ON  OP  CIPSTOMPRS, 
N=SS,  V  -  1 


SOU  IT  I  ON 
NUMB PR 


1  1  Creed)'.) 

1  SO  4 

1  ISS 

2  1  Pa  i  r  ) 

ISO  2 

MU 

3 

180-1 

1  (-11(1 

4 

1792 

1  S9i. 

;■> 

17  SO 

I  192 

(> 

1  7  S  8 

14-0 

7 

1  7  SO 

I  S9<  ■ 

8 

1  092 

1 7-  - : 

9 

ISIS 

l  i-s : 

11) 

1  7(4- 

1  SSS 

PRO  HI.  PM  1  I 'ROB!.  PM 

solution  sounio: 


n  /  v 


ways  to  distribute  the  customers  among  the  V  vehicles.  lor  the  problems 
discussed  immediately  above  this  represents 


1,820 


d i f ferent  poss i hi  1 i t ies . 


Varying  serv ice  parameters .  The  previous  discussion  assumed  that 
only  the  precedence  constraints  were  binding.  Uhen  time  windows  are 
present  as  they  are  in  the  dial-a-ride  service  problem,  the  apparent 
superiority  ot'  the  pair  selection  heuristic  disappears.  Tables  2 a 
through  2()  present  results  tor  1(>  customers  l.\  =  .vdi,  Tour  vehicles  and 
varying  values  of  the  quality  of  service  parameters  (J  and  M.  I  he 
solutions  are  presented  for  three  solution  methods:  greed)',  greed)  pair 
selection,  and  pair  selection.  flearly,  the  individual  greed)'  tours 
should  he  subjected  to  the  optimal  or  r-eptima!  algorithms  to  improve 
the  initial  tour.  Ivhcn  the  quality  of  service  constraints  are  tight, 
as  shown  in  tables  2d  and  21,  the  greedy  selection  procedure  shows  a 
slight  superiority  over  t  pair  selection  technique.  As  the  con¬ 
straints  U-ss  binding,  the  preference  swings  toward  the  pair 

sr  1  ■  .  t  on  ii'.uii  tu.  lii'wrv  rr,  tor  larger  problems  in  the  const  ru  ined 
r  ii  :  :  •  :  >  ■ : :  •  .  *  ‘ .  ■  dat  a  d  ■  not  ■■nrivr!  i  i  t  he  i  T.etir  i  st  i  c  being  superior 
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TABU:  25 


COMPARISON  OB' 
CRLLDY  PAIR 
\ 

SOLUTION  \  A  LI  II  S  1  OR  IIIL  CRLLDY 
AND  PAIR  SPLIT!  1  ION  III.UR1SI  ICS, 
=  38,  V=  l,  Q=3,  M  =  8 

I’ROBl.TM 

CRHIIDY 

GRI -1:1  0  PAIR  P 

AIR  SL Lit!  1  ION 

NIIMBLR 

SOLUTION 

SOLUTION 

SOLDI  ION 

1 

1962 

1982  * 

200  1 

2190 

2084  * 

2138 

.*» 

1891 

1894  * 

1  99  8 

4 

1974 

1958  4 

1  962 

5 

1900 

18  18  * 

19SS 

(.) 

1  Son 

1  80S  * 

2100 

7 

20 1  2 

2  3  ^  > 

25on  * 

8 

2308 

204  8  * 

2112 

9 

1  98  1 

1  '.1  “  2 

1914  ‘ 

10 

2060 

1930  4 

22m) 

‘Best  Solution 


160 


TABI.H  25 

COMPARISON  01-  SOLUTION  VAI.1ILS  I  OR  Till:  GRLLDY, 
GRLLDY  PAIR  AND  PAIR  SLI.LGT  ION  IIUIRISTICS, 

N =35,  V=4,  Q=7,  M=4 


I’ROBLLM 

NIIMBL.R 

c;ri:i;dv 

SOLUTION 

ORIiLDY  PAIR 
SOLUTION 

PAIR  SIT. LOTION 
SOLUTION 

1 

1782 

1716 

1580  * 

1 566 

1548  * 

1038 

.s 

1988 

1  7  10 

1518  * 

•1 

1  764 

1650 

1484  * 

s 

1  668 

1624 

1502  * 

(> 

1  588 

I  r>/>8 

1  568 

7 

1668 

1620  * 

1~24 

8 

1644 

1  526 

1480  * 

9 

1880 

1  60  1 

1562  4 

10 

1772 

1654 

1624  + 

*Bcst  Solution 
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TABU-!  2(> 

COMPARISON  OP  SOLUTION  YAI  IM  S  TOR  Till:  GRPPDY, 
GRPPDY  PAIR  AND  PAIR  SPPPCTIOX  HI.UR  ISTICS, 
N=33,  V=4 ,  Q=1 1  ,  M=(> 


P  ROB  I.  PM 
MJMBPR 

GRPPDY 

SOLUTION 

GRPPDY  PAIR 

SOLUTION 

PAIR  SI!  PUT  ION 
SOLUTION 

1 

2062 

168  1 

16!  K  * 

9 

1518 

1-16  2  * 

1  560 

7> 

1602 

14"0 

1368 

4 

iy  32 

158-1 

1  18  3  * 

5 

1876 

1558  * 

15.84 

6 

1596 

1504 

1  444  * 

7 

1770 

1728 

1580  * 

8 

1  802 

1  322  * 

1458 

0 

1754 

1 624  * 

1  <-36 

10 

1740 

15  "6 

1  1 86 

*Bcst  Solution 


Performance  on  larger  problems.  Tables  27  through  2 1  present 


results  for  30  customers  and  either  three  or  five  vehicles.  The 
results  parallel  those  discussed  above.  Neither  the  pair  selection  nor 
the  greedy  selection  produces  consistently  superior  results.  When  the 
constraints  are  relatively  tight,  the  greedy  selection  procedure  appears 
to  produce  better  results. 

Discussion  of  results.  The  fact  that  a  greedy-like  procedure  often 
produces  the  best  results  is  uncomforting.  Greedy  heuristics  seldom 
produce  the  best  solutions  on  other  difficult  combinatorial  problems. 

The  fact  that  these  solutions  cannot  be  compared  to  the  optimal  is  also 
distressing.  There  is  no  way  of  telling  the  true  effectiveness  of 
these  heuristics  on  other  than  very  small  problems.  The  pair  selection 
heuristic,  as  developed,  and  the  greedy  selection  heuristics  may  well 
be  excellent  ones.  However,  it  seems  likely  that  a  better  technique 
for  determining  which  customers  to  assign  to  what  vehicle  exists. 

What  that  technique  is  or  what  the  basis  is  for  its  development  are  as 
yet  unanswered  questions.  Such  questions  represent  one  area  for  further 
investigation.  This  and  other  areas  that  appear  to  offer  research 
opportunities  that  expand  on  this  work  arc  discussed  in  the  next 
chapter . 
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TABLi:  27 

COMPARISON  OP  SOLUTION  YAI.UPS  TOR  1111:  GREEDY, 
GREliDY  PAIR  AND  PAIR  SELECTION  HEURISTICS, 
N=61 ,  V = 5 ,  Q=5 ,  M=5 


PROBLEM 

NUMBER 

GREEDY 

SOLUTION 

GREEDY  PAIR 
SOLUTION 

PAIR  SELECTION 
SOLUTION 

1 

3696 

3640  * 

3932 

2 

3SS8 

3396  * 

3824 

3 

3780 

353-1 

3656 

4 

3234 

3130  * 

3238 

5 

5354 

3134 

3316 

6 

34  1 0 

3374  * 

354  1 

7 

3036 

3052  * 

3242 

8 

3505 

34(18  * 

34  28 

9 

3662 

3576 

3468  * 

10 

3652 

5622 

3940 

*Best  Solution 


TABI.i:  28 

COMPARISON'  OP  SOLUTION  TAMILS  TOR  Till.  GRLLDY, 
GRLLDY  PAIR  AND  PAIR  SLLLC'P'N  IILURISTICS  , 
N=01,  V = 3 ,  Q=7 ,  M  . 


PROBLLM 

NIJMBLR 

GRLLDY 

SOLDI  ION 

GRLLDY  PAIR 

SOLUTION 

PAIR  SLLLC I I ON 
SOLUTION 

1 

294  0 

2824  * 

28  38 

2 

3230 

2912  * 

2994 

3 

287b 

2730 

2832 

4 

3088 

2870 

26  28  * 

5 

2902 

2606  * 

T  —  -  T 

6 

2872 

2668  * 

2814 

7 

2876 

2646  * 

207b 

8 

2708 

25-2  * 

2'6S 

9 

3  5 1 4 

ji)  22 

2-12  * 

10 

3214 

2930 

2906 

*Bcst  Solution 
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tabu:  si 


COMPARISON  OP  SOLUTION  VAl.lll  S  I'OR 
GRLLPV  PAIR  AN' D  PAIR  SLLLCT I  ON 
N=hl,  V = 5 ,  Q=7,  M=4 


Tiih  (;rp:i:dv, 
iiliiristics, 


PRO BP PM 
NIIMBPR 


GRPHUY 

SOLUTION 


ORI.P.UV  PAIR 
SOLUTION 


PAIR  SPI.I.C  I  ION' 
SOUP!' ION 


CI.APTLK  VII 


SUC.GLS'i'Ci)  ARliAS  IOR  I-IJK  1 1 li.K  Rf.SI.ARCi, 

The  work  documented  in  the  preceding  chapters  represents  one  of 
the  first  attempts  to  define  and  solve  the  pickup  and  delivers  problem. 
As  such,  it  represents  only  a  first  step  in  that  direction.  Mans 
aspects  of  the  problem  remain  to  he  investigated.  Therefore,  this 
chapter  outlines  several  research  opportunities  that  have  suggested 
themselves  during  the  course  of  this  effort. 

I.  SUMMARY  Of  ASPhCTS  STUIMfU 

In  Chapters  I  and  III  the  l’llUP  was  defined  and  discussed  in  its 
most  general  form.  Subsequently,  a  more  restricted  form  of  the  problem 
was  developed  and  explored.  Specifically,  Chapters  IV  through  Y 1  dealt 
with  instances  of  the  Pilin’  where  the  constraints  were  expressablo  in 
terms  of  the  stop  or  sequence  number  of  the  vehicle  serving  a  particu¬ 
lar  customer.  Capacity  was  assumed  to  not  he  a  binding  constraint. 

This  allowed  for  the  exact  solution  of  very  large  single  vehicle 
problems  as  well  as  modest  sized  multiple  vehicle  ones.  The  heuristics 
discussed  considered  the  same  problem  instances  in  order  to  have  a 
precise  measure  of  how  well  a  given  heuristic  performed.  Relaxation 
of  these  restrictions  oiler  the  first  area  of  opport un i t >  . 
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II.  CONSIDLRATIUN  OI:  IIIL  GF.NFRAI.  TROHhh.M 


The  most  general  version  of  the  I’lll'P  is  constrained  by  vehicle 
capacity,  time  windows,  quality  of  service,  operational  considerations 
and  of  course  the  precedence  requirement.  When  these  constraints 
cannot  be  expressed  as  discussed  above,  the  heuristic  techniques 
developed  herein  will  not  work  without  modifications.  'I tie  dynamic 
programming  algorithm  will  not  work  at  all  in  most  cases.  Consequently, 
exact  solutions  to  the  general  problem  appear  doubtful.  Therefore, 
comparison  among  heuristics  probably  will  have  to  be  made  on  a  relative 
basis. 

111.  qlJhSITO.N  OF  1  III;  FIX  I  S  IT. NCI:  OF  F  IAS  1 11 1  F  1  IT 

Given  the  more  general  problem,  one  question  that  must  he  addressed 
is  that  of  the  existence  of  a  feasible  route  or  set  of  routes.  Assuming 
that  an  r-optimal  heuristic  is  to  he  used,  the  existence  of  an  initial 
feasible  solution  is  critical.  Finding  such  a  solution  nay  well  he  an 
extremely  difficult  task  for  those  problem  instances  with  relathely 
tight  constraints.  Use  of  any  tour  construction  heuristic  implicitely 
requires  a  guarantee  that  a  given  partial  tour  can  he  extended  to  a 
feasible,  complete  route.  Such  a  determination  was  often  complicated 
for  the  problems  discussed  in  previous  chapters.  It  would  appear  that 
it  would  be  even  more  complicated  in  the  more  general  case.  Notwith¬ 
standing,  this  represents  the  logical  next  step  in  studying  the  1‘III'P. 
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IV.  MUI.T  1 PIT.  \UiU:i.I  IX  I  INS  IONS 

The  multiple  vehicle  I’UIU’  is  much  more  complex  and  difficult  to 
solve  than  is  the  single  vehicle  PUDP.  It  also  offers  the  greatest 
potential  for  research.  lhe  results  of  Chapter  V 1  suggest  that  better 
solutions  techniques  may  exist  for  the  multiple  vehicle  problem. 


Vehicle  Assignment 

Given  a  good  algorithm  for  solving  the  single  vehicle  1’Uhi’,  the 
key  to  a  good  solution  to  the  multiple  vehicle  problem  lies  in  the 
distribution  of  the  customers  among  the  vehicles.  Although  the 
techniques  discussed  is  Chapter  XI  did  not  perform  as  well  as  hoped,  the 
basic  concept  still  appears  fruitful.  That  is,  obtain  a  good  initial 
assignment  and  then  improve  the  solution  hv  switching  customers  am. nog 
the  vehicles.  lhe  problem  is  in  finding  more  <  ffective  techniques 
for  doing  so. 


Number  P or  V e h  i  cj  e 

The  results  for  the  multiple  vehicle  problem  assumed  that  the 
same  number  of  customers  would  he  assigned  to  each  vehicle.  such  an 
assignment  may  not  lie  the  most  efficient.  Consequently,  relaxing  this 
restriction  offers  anothc  area  for  consideration. 


V.  POSSIUII.I  IV  Of  SI.M  !  PI  Phil's 


In  all  the  c.isos  eon  s  i  d  e  red  ,  there  were  a  !  w ,  i  \  s  taiou  gh  c ;  e- !  .  >  ":i  t 


v;  i\i  11  t 
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that  the  vehicle  was  alwj\  in  use. 


O’  a  I  .It  1  . 1 1  I  ell  w  ,i 


where  the  vehicle  might  be  idle  awaiting  the  earliest  t.me  the 


next  customer  could  be  picked  up.  1. specially  in  the  dial-a-ride 
service  problem,  such  n  scenario  is  not  unrealistic  and  should  be 
considered.  Other  possibilities  tor  slave  periods  might  include  lunch 
breaks  or  coffee  breaks  for  the  drivers. 

VI.  Ol’l.RATIOXAI.  COX'S  I  KA  I  XI  S 

1  he  possibility  of  limitations  on  the  total  time  or  distance  that 
a  vehicle  could  be  operated  is  a  possible  constraint  tint  was  not 
addressed  in  Chapter  VI  .  Such  operational  constraints  could  uppl..  to 
the  single  vehicle  problem,  but  with  optimal  or  near  optimal  selutm:.- 
possible,  imposing  them  could  result  in  no  feasible  solution.  Ihc  :•  o i c 
likely  scenario  for  op  itional  constraints  would  be  in  the  multiple 
vehicle  problem.  It'  and  when  more  effective  assignment  and  i  mprov  erne  lit 
procedures  are  developed  for  the  multiple  \  eh  i  c  1 1  problem,  cousider.it  ion 
of  operational  constraints  would  represent  a  jo-  Pule  suhscqiien t 
ex  tens  ion . 

Vi  i  .  I'OSS  I  HI  I  Ml  I'RISI  It's  !  OR  III!  s  i  pop; 

VI  H n '1.1  S  I'lipr 

Any  number  ot  additional  heuristic  pi  oci-du  o  ■  -  could  be  d  •  re  i  opt  J 
and  appl’ed  to  the  single  vehicle  I’lll'l’.  However,  two  lu  ir  i  - ;  :c-  i,  ; ; 

to  otter  good  potential  for  all  instances  ,  q  tin  i  bbb.  One  i -•  ui 
extension  ot  the  S-optimnl  heuristic,  while  the  ■  - 1 1  v  ■.  r  a  ,  tin  ,,pt  .■  ,  | 


so  1 ut  i on  to  a  r<  i  o  -  1  prof  I  cm . 
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Four  and  5 -Optimal 

For  the  TSP  and  the  VRP,  an  r  value  of  3  in  the  r-optimal  heuristic 
was  usually  the  largest  value  used.  The  reason  is  becuase  of  the 
exponential  increase  in  computations  as  the  value  of  r  increases. 
Consequently,  values  of  r  greater  than  3  were  deemed  practical.  The 
number  of  reconnection  patterns  also  increases  exponentially  as  the 
value  of  r  increases.  The  same  phenomenon  is  true  with  the  PUDP. 
However,  due  to  the  precedence  constraints  of  the  problem,  it  may  not 
be  impractical  to  use  larger  values  of  r.  This  might  be  accomplished 
by  only  considering  those  reconnection  patterns  that  offer  the  highest 
probability  of  obtaining  a  feasible  reconnection.  The  range  of  stop 
numbers  included  in  each  loop  might  also  be  limited  so  as  to  include 
only  those  stops  which  are  most  likely  to  yield  a  feasible  reconnection 
pattern.  For  the  3-optimal  solution,  this  range  was  taken  to  be  Q 
stops.  Given  the  demonstrated  superiority  of  the  3-optimal  solution 
over  the  other  heuristics  tested,  a  4  or  5-optimal  heuristic  appears  to 
be  a  most  fruitful  area  to  investigate. 

Optimal  Solution  to  Related  Problem 

For  the  general  problem  where  dynamic  programming  cannot  be  used 
to  obtain  the  exact  solution,  it  may  be  possible  to  use  dynamic 
programming  to  obtain  a  good  solution.  This  could  be  accomplished  by 
approximating  the  general  constraints  by  stop  numbers  and  then  solving 
the  related  problem  by  the  dynamic  programming  algorithm.  If  the 
resulting  solution  is  feasible  to  the  original  problem,  one  would 
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hopefully  have  a  good  solution.  If  the  solution  were  not  feasible,  some 
interchange  routine  would  be  needed  to  find  a  feasible  solution  with  a 
minimum  of  additional  cost.  Intuitively,  this  procedure  should  perform 
well  and  could  be  used  on  either  the  single  vehicle  problem  or  as  the 
routing  portion  of  a  two  step  solution  of  the  multiple  vehicle  PUDP. 

VIII.  PUDP  POTENTIAL 

The  PUDP  is  a  new  problem  which  is  only  just  beginning  to  attract 
research  attention.  Given  both  the  complexity  of  the  problem,  which 
makes  obtaining  solutions  difficult,  and  the  practicality  of  the 
applications,  the  PUDP  should  appeal  both  to  the  theoretician  and  to 
the  practitioner  for  some  time  to  come.  The  work  documented  herein 
as  well  as  the  few  related  articles  represent  only  the  tip  of  the 
iceberg.  A  great  deal  remains  to  be  accomplished  before  the  PUDP  can 
be  considered  fully  solved. 
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